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FOREWORD 



This document is the final report for Contract NAS8-28614, "Design 
Criteria for Low Profile Flange Calculations, " needed for the establishment 
of a Low Profile Flange Standard, Low Profile Flanges are characterized 
by featuring a small width but large height. Testing of Low Profile Flanges 
showed their superiority in performance, weight, and envelope volume in 
comparison to commonly used flanges for space application. 

The work was initiated by the Layout and Assembly Engineering Branch, 
EngineeringDivision, Astronautics Laboratory of the NASA, Marshall Space 
Flight Center, in a joint effort with the Lockheed Missiles and Space Company, 
Inc,, Huntsville Research and Engineering Center. 

The primary objective of this effort was to evaluate the existing design 
procedure shown in the publication "Application of Low Profile Flange Design 
for Space Vehicles, " and other flange design literature to establish a standard 
for Low Profile Flange calculations. 


The period of performance of this study was from May. 18, 1972, to 
March 22, 1973. 



NASA 

Marshall Space Flight Center 
Huntsville, Alabama 
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Section 1 
INTRODUCTION 

The purpose of this study was to develop an analytical method and a 
design procedure to design flanged separable pipe connectors based on the 
previously established algorism for calculating low profile flanges. These 
flanges demonstrated their superiority with respect to leak- tightness and 
weight savings in comparison with other flanges used for space application. 

When the low profile flange was first considered for space vehicle and 
launch application no design procedure was established and conventional flange 
design methods were used for the basic analysis. To remedy the situation 
Prasthofer (Ref. 1) devised a simple but effective design procedure considering 
the strength of the flange ring cross section as the design criteria. 

It has been shown by Schwaigerer (Ref. 2) and through experiments by 
Buhner et al. (Ref. 3) and Haenle (Ref. 4), that there is a major contribution 
of the adjacent tube wall to the strength of a flanged connection. If one plots 
the flange roll angle x versus the applied moment one finds a gradual decline 
of the slope of the curve (Fig. 1-1). This points to the existence of a plastic 
hinge at the most highly stressed section of the tube. The location of the 
plastic hinge is close to the neck of the flange, depending, in part, on the 
variation of the wall thickness of the tube in the area of the neck (Fig. 1-2). 
Buhner et al. (Ref. 3) present a large number of data relating to the perform- 
ance of flanges after the formation of the plastic hinge. The comparison for 
flanges with identical cross-section (Fig. 1-3) reveals that the next best choice 
to a conventional design with conical hub is one with a fillet (c). This com- 
parison is not realistic, though, since design (b) can be replaced by a much 
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C onf i gu r ati on 

(a) (b) (c) (d) 



0.2% permanent set at ® 

❖ 

1 deg permanent roll X 


Fig. 1-3 - Comparison of Structural Performance Based 
on Identical Flange Cross Section (Ref. 3) 
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narrower one of type (c), thus reducing the applied moment and therefore 
not requiring the larger moment capacity available with type (b). This is 
one of the advantages of the low profile flange which leads to the attendant 
weight saving, weight being proportional to the cross-sectional area and the 
centroidal radius. It also should be remembered that the strength of a flange 
increases approximately linearly with the flange width, b, but quadratically 
with the flange height, h, while the stiffness (resistance against roll) increases 
even cubically with h. This explains the better performance of the low profile 
flange over conventional wide profile flanges. 

Thus, the advantage of the low profile flange is seen as being twofold: 
first, to reduce the lever arm, e, between the gasket and bolt circle, and 
second, to have the material of the flange where it is most effective, i.e., 
have the height, h, larger than the width, b (see Fig. 1-2). 

Most of the available data on flange performance in the plastic range 
has been devoted to designs in steel at moderate temperatures. Steel has, 
however, a distinct yield point in its stress- strain diagram as compared to 
aluminum or titanium. The development of a plastic hinge for the latter 
materials at different temperatures would be a most interesting subject for 
further experimental investigations since the flange design method in this 
report is partially based on the assumption of a plastic hinge. 

The plastic design method has been made part of the German flange design 
code DIN 2505 (Ref. 5), whereas American practice is based on an elasticity 
approach (Refs. 6 and 7). The use of the plastic design method is valid when 
the material is capable of undergoing large strains -without fracture. The 
plastic method assumes a ductile failure. If a brittle failure is the predominant 
mode such as for certain high strength steels then the elastic method is more 
suitable. 
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1.1 DESIGN CRITERIA 

The condition for sufficient strength of a structural component requires 

that 


max 


K 

(F.S.) 


( 1 . 1 ) 


where a at the maximum equivalent stress, K at the reference strength 

'max 

of the material, and (F,S.) is the factor of safety (with or without subscripts)* 
In this paragraph these three quantities are briefly reviewed* 

(a) Maximum Equivalent Stress : The computation of the maximum equivalent 

stress, a , to be compared with the uniaxial material strength is based 

max 

on the type of expected failure. For a failure associated with plastic deforma- 
tion (yielding) or fatigue, the hypothesis of the limit of the elastic distribution 
energy by Huber (Ref. 8) and vonMises (Ref. 9) is used. The equivalent stress 

is 


_ yr ^ CT 1 ' ct 2^ 2 + + ( a 3 " CT i) 2 (1.2) 

where cr-,. are the principal stresses. For components subjected to 
high tensile stresses, i.e., if ct^X), the equivalent stress is 

a e = CT 1 


This failure mode is fracture. 
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A third equivalent stress occasionally considered is the one defined by 
Tresca {Ref. 10) and is used for shear failures, 

(7 = <7 - a . = 2t 

e max min max 


In the development of the flange design procedure the Huber and von Mises 
hypothesis is used. 


(b) Material Strength : The material strength K to be used in Eq. (1.1) 

depends on the type of failure envisioned and must correspond to the type of 
equivalent stress computed. The two most frequently encountered 


max 

types of uniaxial stress -strain diagrams are shown on Fig. 1-4. Diagram 

(a) has a distinct yield point with the tensile yield strength F^, . The ultimate 

tensile strength is F^. Diagram (b), on the other hand, has a gradual change 

in slope requiring the definition of a yield strength from permanent strain 

considerations. Typically, the yield strength is F^. = F where F ^ is the 

0.2% stress at permanent set. If F^ is much larger than F ^ the definition 

of the yield strength may be based on F Q 5 or F 1 This is the case with 

highly ductile materials. For the subsequent use in a design formula the 

stress -strain diagram is replaced by an idealized diagram as shown on 

Fig. 1-5 which could be called elastic — ideally plastic. This diagram must 

specify, however, a limit of its validity by giving a maximum allowable strain, 

c . A component which has been designed according to a plastic design 
mix 

method, such as the one proposed in this report for flanges, needs to be 
checked for strains under the design conditions, i.e., the ultimate load. This 
load condition will be discussed later in more detail. 


When temperature effects are to be considered, the appropriate strength 
values at the design temperature must be used. Similarly, fatigue strength 
and creep rupture strength can be the dominant strength values to be considered. 
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Fig. 1-5 - Idealized Stress -Strain Diagram 
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{c) Safety Factors : The proper use of the safety factors is important in a 
complex system such as a bolted connection, but it also leaves room for 
different design philosophies. For a pipe system three pressure levels are 
considered: operating pressure, proof pressure and burst pressure. Usually 
the proof pressure is times the operating pressure, and the burst pressure 
two times the operating pressure. These factors are implied safety factors 
against uncertainties in the prediction of the operating pressure due to pres- 
sure surges at valve closure or vehicle vibrations accompanied by pressure 
oscillations. The design pressure is mostly chosen as to be proof pressure, 
that is, the structure is to be able to withstand the proof pressure without 
damage. That condition occurs at least once in the lifetime of the structure. 

If an elastic design method is employed and only stress peaks are checked, a 
small safety factor of say 1.2 against yielding at critical points is sufficient. 
For the plastic design method, however, instead of a safety factor an ultimate 
factor of at least 1.5 is used by which the load is multiplied. This magnified 
load is called the ultimate load. For example, the maximum applied moment 
on the flange due to the proof pressure condition is nip and the ultimate 
moment that is to be carried is 


m Fu = (F - S * } m F 


(1.5) 


A structural capacity has to be provided for m p u * This capacity can be 
expressed as 


m r = Z„ F. 
Fu F ty 


( 1 . 6 ) 


where F^ is the tensile yield strength of the material and Zp is the combined 
section modulus of the flange and the adjacent tube after a plastic hinge has 
formed in the neck. 
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The flange is then still in an elastic state of stress with only the extreme 
fibers yielded. Schwaigerer suggests in Ref. 2 that the flange cross section, 
too, should be considered being in a fully plastic state of stress. This condi- 
tion is, however, somewhat exaggerated. 

Other safety factors are needed to cover uncertainties such as in the 
computation of the stresses (using average values) and uncertainties in the 
material properties, i.e., the values of K chosen for the different materials. 
Possibly the material properties of the flange and the bolts are much more 
accurately known than those of the gasket, requiring a higher gasket safety 
factor. 

The total safety factor may be defined as the product of the individual 
safety factors 


(F.S.) =<F.S.) r (F.S.) 2 . (F.S.) 3 . 


(F.S.) 


n 


(1.7) 


In this study the values for the safety factors have been chosen more or less arbi- 
trarily. Also, some design formulas such as the ones for the tube wall thickness 
contain implied safety factors. These are explained where they occur in Section 2. 

1.2 PAST EXPERIENCE WITH LOW PROFILE FLANGES 

The initial idea for the low profile flange concept came from Boon and 
Lok (Ref. 11) which was taken up later by Prasthofer (Ref. 12) for the design 
of launch vehicle pipe connections. Qualification testing reported by 
(Ref. 13) and experimental stress analysis of one photoelastic model configura- 
tion by Kubitza and Hearne (Ref. 14) showed the soundness of this flange 
concept. Design procedures for a similar type of flange have been established 
by Trainer et al. (Ref. 15), by Aerojet General (Ref. 16) and Pratt and Whitney 
(Ref. 17) although not in a usable form. The latter procedures are tailored to 
specific seal configurations and are therefore not generally applicable. In 
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addition most of the existing methods require an excessive amount of computa- 
tions if carried out manually. An automated design study reported by Rathbun 
{Ref. 18) is set up to produce wide profile flanges with conical hubs, being 
undesirable in the context of the low profile concept. 

Previous design methods were not definite on the minimum spacing re- 
quirements for the bolts. Bolt spacing is the driving parameter for the flange 
width. Minimum bolt spacing assures a low profile flange. The spacing re- 
quirements are discussed in more detail in Section 2. It should be noted here 
that the bolts should be as small as possible and be located as closely to the 
tube wall as can be accomplished within the constraint of wrench clearance 
requirements. This can be accomplished by providing countersunk spot faces 
for the bolts in Configuration (c) on Fig. 1-3 or by a machined groove as in 
Configuration (d). The latter sacrifices some ultimate moment capacity. 
Configuration (c) has been used exclusively, so far, in all past low profile 
flange designs. The introduction of stress peaks around the spot faces has 
been impressively demonstrated in Ref. 14. It is therefore suggested to supply 
contoured washers that would eliminate countersinking, possibly even with 
spherical glide surfaces to accommodate rotations of the bolt head or nut with 
respect to the flange. This brief summary may suffice to characterize past 
experience. 

Seen in the light of the long history of flange design and analysis methods, 
beginning with Westphal's classical paper (Ref. 19) of 1897, the methods 
presented in this report constitute the logical extension of current ideas. 
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Section 2 

DESIGN PROCEDURE 

In this section the individual steips of the design procedure are derived. 
Some of the steps have several alternatives and the most suitable ones are 
selected. A summary of the design procedure to serve as a guideline for 
manual computations and as an outline for the design section of the computer 
program is given in Appendix A. 

2.1 TUBE DESIGN 


The computation of the required wall thickness for a cylindrical tube 
under internal pressure is based on the stresses. For a thick-walled tube 
these are 


a 


r 


(b/r) 2 + 1 

P 

(b/a) 2 - 1 


( 2 . 1 ) 


a 


1 


(b/a ) 2 - 1 


= (b/r f ± 1 

u (p , 9 ^ 

* (b/a) 2 - 1 


( 2 . 2 ) 


(2.3) 


The coordinate system is defined on Fig. 2-1, The equivalent stress is 

.2 


a = 

e 


V2 


=- V <CT r - V * + ( V - + - «r> 2 ' ■^7 J / £i r P ■ < 2 - 4 > 

(b/a) - 1 
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Fig. 2-1 - Coordinate System for a Thick-Walled Tube 
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with a maximum at r = a, 


- P . 

e max (b/a) - 1 


(2.5) 


If cr^ is larger than given by Eq. (2.2), say a ^ > a then 


o = 

e 


V 3 (b/r ) 4 + (1 - <^/a x ) 2 


(b/a)' 


P . 


( 2 . 6 ) 


The average equivalent stress is generally 

= (b/a) 2 + 1 £ 
av (b/a) - 1 


(2.7) 


When the equivalent stress a according to Eq. (2.4) is equal to the yield 
strength F t which is assumed to be a constant for mathematical simplicity, 
then 



The equivalent stress in this case (Ref. 20, p. 106) is 



( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 
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The fully elastic and fully plastic states of stress are shown on Fig. 2-2. The 
reversal of the stresses when going from the elastic to the plastic state of 
stress is obvious. 


In Ref. 2, p. 29, it has been shown that up to a ratio 


b 

a 


1.2 


( 2 . 12 ) 


the average equivalent stress a can be used for the design of a tube since 

e av 

it is almost equal to the maximum equivalent stress o g and to the equivalent 
stress of a fully plastic state, i.e., max 


where 


and 


-2 — ~ —2— ^ -2- if -<1.2 

e °e o a 

max av e 



p _ 2t 

a 2a + t 

e 

av 


(2.13) 


(2.14) 


(2.15) 


( 2 . 16 ) 


The latter equation (2.16) was derived from Eq. (2.7) with t = b-a, where t is 
the wall thickness. When the strength design criterion Eq. (1.1), is applied 
the relation becomes 


2t K 

p = 2a + t (F.S.) 


(2.17) 
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which can be solved for the wall thickness, giving 

t E5 (2.18) 

* " £ • 

(F.S.) 2 

This formula is the basis for most pressure vessel design codes, for example, 
the American "ASME Pressure Vessel Code" (Ref. 21) or the German 
"Dampfkessel-Bestimmungen" (Ref. 22). 


In order to accommodate wall thickness tolerances At and a factor 4* 
weakening by welds the wall thickness t in Eq. (2.12) is replaced by 


t = (t - At) 4* 


(2.19) 


which gives the formula for the thickness as 


t = 


E* r 

7 .11— . . E\ 

\(F.S.) 2 


+ At 


* 


( 2 . 20 ) 


This is the formula used in Ref. 22. 


In Ref. 21 this equation has been modified by taking 


t = 


1.1 P a 


K 


(F.S.) 


>4 P 


( 2 . 21 ) 


while a formula used by Pratt & Whitney (Ref. 17) is 


t = 


P_a 


~ K ~ 

(F.S. I 


+ 2 At 




( 2 . 22 ) 
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The simplest formula, based only on circumferential stress, is 


t 



(F.S.) 


(2.23) 


The weakening factor 4 1 has been used in the order of 

0.70 < iJj < 1.00 (2.24) 

For the weakening by a weld the stress component perpendicular to the weld is 
most important. If a weld is at an angle y with the cylinder axis then the stress 
perpendicular (normal) to the weld is given by 

a = a sin ^y + a. cos^y (2.25) 

n x (p 

Therefore in Eqs . (2.19) and (2.20) the weakening factor is generally 

(2.26) 

l + cos^y 

The factor <|/ has to be determined by test. 

Considerations other than internal pressure, such as creep, vibrations, 
bending and shear, may influence tube design. These are briefly reviewed in 
the following paragraphs. 

Creep : To simplify the derivation, only steady state creep is considered. 

This problem was studied in depth by several investigators (Refs. 23 through 27). 
Let the material law be given by 


e 


e 



where c e is the equivalent strain defined by 


(2.27) 
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= V (e r ‘ e c/ + (e ^ " a x )2 * ^ 


and the creep constant B is 


(2.28) 


b = pe 


a T 


(2.29) 


where a and P are independent of temperature, T is the absolute temperature 
and G is the base of the natural logarithm. The stresses in this case, corre- 
sponding to Eqs. (2.1) through (2.3) are given by 


a 


r 



(2.30) 



(2.31) 


(2.32) 


It can be seen that for n = 1 the elastic case of Eqs. (2.1) through (2.3) is 
obtained. Using the foregoing relations the accumulation of strains can be 
computed for the lifetime of the tube, thus serving as a design criterion for 
the selection of the tube thickness. 


From Eqs. (2.30) through (2.32) the maximum equivalent stress at the 
inside of the tube is, similar to Eq. (2.5), 
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(2.33) 


To design a tube for a given lifetime until creep rupture occurs, the ultimate 
equivalent strength is computed by 

o = (F.S.) a (2.34) 

U e 

max max 


and from a plot of the equivalent stress versus the creep parameter, P, the 

value of P for a of Eq. (2.34) is obtained, 

max 

The creep parameter, if for example, Larson and Miller's (Ref. 28) 
formulation is used, is defined by 

P = C,T [log (t rupt ) + c 2 j (2-35) 

which can be solved for t^ . giving 

‘rupt = antilog (^T - C z) • (2 ’ 36) 

This excursion into creep analysis methods may suffice. 

Vibration : The oscillations involving propellant feedlines, engines and 

longitudinal structural modes of a launch vehicle are described in a paper by 
Ryan et al. (Ref. 29). To cope with the problem from a design point of view 
the following approach may be taken 

t = (A BlS + c (2.37) 

K F Ap 
(F.S.) " 2 
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where Kj. is the appropriate £atigue strength of the tube material for the stress 
ratio R = 0 assuming a maximum internal pressure of 


Ap 


max 


mm 


2.38) 


The maximum pressure p would be determined from a vibration analysis 

max 

such as the one cited. 


Bending and Shear ; The presence of a bending moment and a shear force 
introduce stresses into the tube wall which may control the design. The ex- 
pressions for the stresses in terms of a bending moment Mj and a shear force 
are 

M 1 

o = 5 - cos0 (2.39) 

X n R £ t 


and 


'1 


xip 


JTRt 


siaq> 


(2.40) 


where if) is the circumferential coordinate. From Eq. (2.39) an equivalent 
axial force of 


M 

n = , (2.41) 

X • TTR^ 

being the maximum, should be used in addition to the axial stresses resulting 
from internal pressure. \ 

2.2 BOLT SIZE 

The essential idea of the low profile flange concept is to have a large 
number of small diameter bolts. Therefore, the design calculation is started 
by selecting a bolt diameter and then trying to accommodate the number of 
bolts required to keep a leakproof connection. To find a basis for selecting the 
bolt diameter it is assumed that the wall thickness computed previously can 
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support an internal pressure of 

P 

which is obtained from the simplified formula given by Eq. (2.23), K^, being 
the tube strength. On the other hand the entire pressure has also to be carried 
by the bolts, i.e., the bolt force 


K_ 


a (F.S.) r 


(2.42) 


TD 2 

P B = >ra p 

has to be equal to 

„ Z ” r B ’ A l K B 
B " s 4 (F.S.) b 


(2.43) 


(2.44) 


where is the bolt circle radius, s is the spacing of the bolts, d B is the 
nominal bolt diameter and K B the bolt strength. The bolt strength is usually 
chosen to be the ultimate tensile strength, together with the appropriate safety 
factor. When the pressure value from Eq. (2.42) is substituted into Eq. (2.43), 



1 r a 


K, 


a (F.S.) t 


(2.45) 


and Eqs. (2.43) and (2.44) set equal, then 


7T a 


t K T 


a (F.S.). 


2 7T r fi 7Td B K b 


4 (F.S.) B * 


or, after rearranging, 


t 


K r 


TC d 


a ^:>T 


B r B 


K 


B 


(F.S.K 


(2.46) 


(2.47) 
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An estimate for the ratio d^/t is arrived at by assuming the following ratios 


k b / k t ^ 

(F.S.) b /(F.S.) t ~ * 


(2.48) 



(2.49) 



(2.50) 


so that 



2 

7T 


(2.5) 


_1 1 _ 

(1.1) 1.5 


1.0 . 


(2.51) 


This is only an initial estimate to get the design calculations started. 
The final bolt diameter will be determined when checked against all the other 
requirements to be discussed later. 


When the initial selection of the bolt diameter has been made the require- 
ments for wrench clearance and adequate spacing of the bolts from each other 
and from the edge of the flange have to be considered. On Figs. 2-3 through 2-5, 


and Tables 2-1 through 2-3, which were 
values for 



^2 = 


taken from Ref. 18, non-dimensional 


s / d B 

(2.52) 

e l/ d B 

(2.53) 

e 2 /,d B 

(2.54) 
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Dimensionles s 
Parameter (Minimum) 



Figure 2-3 - Design Parameters for Open-End Wrenching 
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Dimensionless 
Parameter {Minimum) 



Fig. 2-4 - Design Parameter for Socket Wrenching 
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Table 2-1 

BOLT TABLE FOR OPEN-END WRENCHING (REF. 18) 


Size 

d B (in.) 

% 

H 

^2 

A oB (in 2 ) 

1 

0.2500 

3.00 

2.00 

1.50 

0.03182 

2 

0.3125 

2.60 

1.80 

1.40 

0.05243 

3 

0.3750 

2.67 

1.67 

1.33 

0.07749 

4 

0.4375 

2.57 

1.57 

1.29 

0.10631 

5 

0.5000 

2.50 

1.62 

1.24 

0.14190 

6 

0.5625 

2.45 

1.56 

1.22 

0.18194 


0.62 50 

2.40 

1.50 

1.20 

0.22600 


0.7500 

2.33 

1.49 

1.08 

0.33446 


0.8750 

2.35 

1.43 

1.07 

0.46173 


1.0000 

2.25 

1.37 

1.06 

0.60574 


1.1250 

2.22 

1.33 

1.00 

0.76327 


1.2500 

2.25 

1.40 


0.92905 


1.3750 

2.23 

1.36 


1.15488 


1.5000 

2.17 

1.33 

1.00 

1.40525 


Legend: 
Size = 
d B = 
^0 = 

*1 = 

^2 = 

A OB = 
d hole 


size number of the bolt 
nominal diameter of the bolt 


‘B 


V spacing parameter (dimensionless) 


B 


B 

stress area of one bolt 


dg + 0. 005 in. 
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Table 2-2 

BOLT TABLE FOR SOCKET WRENCHING (REF. 18) 




H 

^1 

EM 

1 

0.2500 

2.76 

■ 

1.40 

2 

0.3125 

2.53 

1 

1.28 

3 

0.3750 

2.37 

1.33 

1.20 

4 

0.4375 

2.26 

1.25 

1.14 

5 

0.5000 

2.18 

1.20 

1.10 

6 

0.5625 

2.20 

1.22 

1.11 

7 

0.62 50 

2.22 

1.25 

1.12 

8 

0.7500 

2.12 

1.17 

1.07 

9 

0.8750 

2.28 

1.31 

1.14 

10 

1.0000 

2.19 

1.25 

1.10 

11 

1.1250 

2.14 

1.22 

1.07 

12 

1.2500 

2.09 

1.18 

1.04 

13 

1.3750 

2.00 

1.16 

1.02 

14 

1.5000 

2.02 

1.13 

1.00 


Legend: 
Size - 



A oB = 
d hole " 


size number of the bolt 
nominal diameter of the bolt 


B 


-r— } spacing parameter (dimensionless) 


B 


B 


stress area of one bolt (see Table 2-1) 


d D + 0.005 in. 

-D 
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Table 2-3 

BOLT TABLE FOR INTERNAL WRENCHING (REF. 18) 


Size 

d B (in.) 

^0 

*1 

^2 

1 

0.2500 

1.92 

1.16 

0.96 

2 

0.3125 

1.86 

1.09 

0.93 

3 

0.3750 

1.79 

1.04 

0.91 

4 

0.4375 

1.80 

1.03 

0.91 

5 

0.5000 

1.78 

1.00 

0.90 

6 

0.5625 

1.76 

0.98 

0.89 

7 

0.6250 

1.75 

0.96 

0.88 

8 

0.7500 

1.68 

0.91 

0.84 

9 

0.8750 

1.69 

0.90 

0.85 

10 

1.0000 

1.67 

0.89 

0.84 

11 

1.1250 

1.86 

0.96 

0.92 

12 

1.2500 

1.67 

0.87 

0.83 

13 

1.3 7 50 

1.80 

0.93 

0.89 

14 

1.5000 

1.65 

0.85 

0.82 


Legend: 

Size 

d B 


^0 

" d B 

’ll 

e l 
" d B 

^2 

e 2 
' d B 


= size number of the bolt 
= nominal diameter of the bolt 
S 


= -j— \ spacing parameter (dimensionless) 


oB 


= stress area of one bolt (see Table 2-1) 


hole 


dg + 0.005 in. 
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are given. The distances e^ and refer to a flange with machined spot- 

faces and are shown on Fig. 2-6. Also, the tables contain the stress area 

A of each bolt size. These were taken from Ref. 30 for the ISO-inch coarse 
oB 

thread series for l/4 ^dg _< 3/2 inch. The corresponding recommended 
metric series is given in Tables 2-4 through 2-6, where 6.3^dg^36.0 mm 
(0.2480 < d_ < 1.4173 inch). Where 14 different sizes were used for the indi- 

JD — 

cated diameter range in the ISO-inch series, only nine different sizes are 
given for the metric series. This series, however, is tentative and subject 
to further studies by the Industrial Fasteners Institute in Cleveland, Ohio. 

The bolt tables given are not to be taken as definite data. They were 
merely used for the numerical example problems of this project. The design 
procedure and the corresponding program are configured to allow additional 
bolt data to be incorporated such as data for 8 or 12 point heads. 


The diameter of the bolt hole is taken as d^ o ^ e = dg + 0.005 inch (+0.1 
mm). These clearances have been assumed to be able to compute numerical 
examples and are not to be taken as definite design data. 


The spot face diameter is assumed as ^ S p 0 £ = 2e^, where e^ is given 

in the bolt tables. A fillet radius of r , = 0.062 inch (1.5 mm) is provided. 

spot 


When a machined groove is selected both distances are e^ = e^ = rj ^ dg 
as shown on Fig. 2-7. 


The selection of the fillet radius on Fig. 2-6 and the groove radius on 
Fig. 2-7 is somewhat arbitrary. While the machined groove is intended to 
reduce stress concentrations due to notch effects at the neck, it cannot reduce 
the high stresses in the cylinder portion. The fillet on Fig. 2-6 is intended to 
do this. A basis for the size of the fillet radius can be found by considering the 
wavelength L of the stress pattern along the shell meridian. This stress pat- 
tern alternates sinusoidally with exponentially decreasing amplitudes. The 
ratio of two successive amplitudes, considering only the edge disturbance 
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Table 2-4 

METRIC BOLT TABLE FOR OPEN-END WRENCHING (REF. 18) 


i . 
size 

^ d B 
(mm) 

d B 
(in. ) 

^0 

^1 

^2 

a °b 

(mm 6 ) 

A „ 
(in. Z ) 

i 

6.300 

0.2480 

3.00 

2.00 

1.50 

22.276 

0.035 

2 

8.000 

0.3150 

2.85 

1.90 

1.43 

36.126 

0.055 

3 

10.000 

0.3937 

2.70 

1.80 

1.85 

57.261 

0.089 

4 

12.500 

0.4921 

2.59 

1.68 

1.28 

91.524 

0.142 

5 

16.000 

0.6299 

2.45 

1.57 

1.18 

155.070 

0.240 

6 

20.000 

0.7874 

2.35 

1.48 

1.12 

242.297 

0.375 

7 

25.000 

0.9843 

2.28 

1.40 

1.07 

382.801 

0.593 

8 

30.000 

1.1811 

2.23 

1.35 

1.02 

555.296 

0.861 

9 

36.000 



1.4173 

2.19 

1.33 

1.01 

809.423 

1.255 


Table 2-5 

METRIC BOLT TABLE FOR SOCKET WRENCHING (REF. 18) 


i . 
size 

d B 

(mm) 

^0 

^1 

^2 

i 

6.300 

2.80 

1.60 

1.40 

2 

8.00 

2.69 

1.53 

1.36 

3 

10.000 

2.57 

1.48 

1.29 

4 

12.500 

2.46 

1.41 

1.23 

5 

16.00 

2.34 

1.33 

1.16 

6 

20.000 

2.25 

1.27 

1.11 

7 

25.000 

2.15 

1.22 

1.07 

8 

30.000 

2.08 

1.17 

1.03 

9 

36.000 

2.02 

1,14 

1.01 
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Table 2-6 

METRIC BOLT TABLE FOR INTERNAL WRENCHING (REF. 18} 


i . 
size 

d B 

(mm) 

^0 

^1 

^2 

i 

6.300 


1.16 

0.96 

2 

8.00 


1.12 

0.94 

3 

10.000 

1.85 

1.08 

0.92 

4 

12.500 

1.81 

1.03 

0.90 

5 

16.000 

1.77 

0.99 

0.87 

6 

20.000 

1.73 

0.94 

0.85 

7 

25.000 

1.70 

0.91 

0.84 

8 

30.000 

1.68 

0.88 

0.83 

9 

36.000 

1.66 

0.86 

0.82 
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t 



Fig. 2-6 - Low Profile Flange with Machined Spot Faces 
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introduced by the flange, is 


where 



(2.55) 


(2.56) 


The radius of the cylinder middle surface, r , is 

r = a + t/2 . (2.57) 

o 

From the logarithmic decrement 

pL/r Q = Zn (2.58) 

the wavelength is 

L = 27 T r Q /p (2.59) 

The fillet radius should cover approximately one-eighth of this wavelength 
in order to reduce the shell stresses at the neck. Equations (2.55) through 
(2.58) are illustrated on Fig. 2-8, To simplify the design procedure, approxi- 
mate fillet and groove radii are listed in Tables 2-7 and 2-8, respectively. 


2.3 BOLT CIRCLE RADIUS AND FLANGE WIDTH 

The magnitude of the bolt circle radius and the flange width are deter- 
mined by the space required on the upper surface as shown on Figs. 2-6 and 
2-7. In the case of machined spot faces a minimum distance c^ from the tube 
wall is maintained to accommodate the tool for making the spot face. The 
formulas for the bolt circle radius can be written for the machined spot face, 


2-24 


LOCKHEED • HUNTSVILLE RESEARCH 4 ENGINEERING CENTER 



LMSC-HREC TR D306492 


Stress Amplitude 

* 





LMSC-HREC TR D306492 


Table 2-7 

FILLET RADIUS FOR FLANGE WITH MACHINED SPOT FACES 


t 



(in. ) 

(mm) 

(in. ) 

J.1 

(mm) 

> 0.2 

5 

0.3750 

10 

< 0.2 

5 

0.3125 

8 

< 0.15 

3 

0.2500 

6 

< 0.10 

2.5 

0.1875 

1 

4 

< 0.05 

1 

0.1250 

3 


Legend: 

= tube thickness 
= fillet radius 



Table 2-8 

GROOVE RADIUS FOR FLANGE WITH MACHINED GROOVE 


t 


_ _ r fil 

(in.) 

(mm) 

(in. ) 

(mm) 

>0.2 

5 

0.1250 

3 

< 0.2 

5 

0.1000 

2.5 

< 0.15 

3 

0.0750 

2. 

< 0.10 

2.5 

0.0500 

1 . 

< 0.05 

1 

0.0250 

0.5 


Legend: 


t = tube thickness 
r,.. = groove radius 
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Fig. 2-6, as 


r B = r. + t + + e x (2.60) 

where = 0.0625 inch or 1.5 mm was assumed, and for a machined groove. 
Fig. 2-7, as 

r B = r j + ^ + 2r^ + . (2.6 1) 


The flange width is then 


b = r B + e 2 " r i * (2.62) 

The inside radius of the tube, a, is denoted by r. in this and the following 
sections . 

2.4 GASKET 


In selecting the gasket and computing the required contact force, two 
phases must be considered. The first phase is the initial precompression 
phase for which a total flange force of 

P .G ^ = 27rr G S G^ (2.63) 


is required, where is the corresponding line load per unit length of the 

centerline of the contact surface. The radius of this centerline is r . The 
second phase is the operational phase in which a certain minimum contact 
force is to be maintained to have zero leakage. This contact force is written 
as 

p a ' = 2 * r G s g ' <2.m> 
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The total required flange force Pp is the sum of the force caused by internal 
pressure, P , 

P = 7 T r?, p (2.65) 

p G ^ 

( 2 ) 

and the contact force P' , 


p(2) 

F 


2 

rrr G P 


+ 27 r r G 


c(2) 


( 2 . 66 ) 


It is important to understand what ^ and are for various gaskets 

and how they are related to the interface leakage rate. Starting out from a 
macroscopic view, i.e,, looking at the whole flange assembly, the relations 
between the flange force and the internal pressure at which a given leakage 
occurs are shown on Fig. 2-9. In the low pressure range a nonlinear relation 
exists between the initial precompression force P^ ^ and the corresponding 
pressure. When the force P^ has been reached, as given by Eq. (2.63), this 
relation becomes linear as the pressure increases. As the pressure is reduced 
from above p^ the relation remains linear all the way to p =0. Under renewed 
pressurization the relation remains linear. Thus P^. ^ has been established as 
the minimum load for precompression of the material. The initial precompres- 
sion force P^ in terms of P^ and P^ is approximately 

p(0) _ a p U) + (i _ a ) (2.67) 

The coefficient a should be selected to match the test data. 


To characterize a gasket material two numbers are needed. First, a 
number characterizing P^ ^ > snd second, a number characterizing the slope 
of the straight line. 
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For flat gaskets the line loads are 


!<»> . 


b iff a G 


( 2 . 68 ) 


and 


s' 2 ' = b' 2 »k p 

G eff p v 


(2.69) 


For other than flat gaskets the line load is given directly while 


5 ( 2 ) _ 

>G “ 


K p 

P 


(2.70) 


The quantity o_ has the units of stress and k is dimensionless . The quantity 
G P (i ) /2) 

has the dimensions of a length. The effective widths and b^' depend 

on the shape of the interface, such as tongue and groove, etc., and the type 

of gasket, such as serrated, etc. 


Data for conventional applications can be found in Ref. 6. For cryogenic 
or storable propellants in liquid or gaseous form these data are not readily 
available and will have to be obtained from testing or be established from data 
not currently in this format. 


One source of information is the comprehensive study by Bauer et al. 
(Ref. 31), which contains data on interface leakage as related to material 
hardness, contact surface topography (surface finish) and contact stress, 
expressed as 



m 

a 


(2.71) 


3 

where h is the "conductance parameter", K g is a constant and a is the 
contact stress. The exponent m depends on both the material hardness and 
the surface topography. This relation of Eq. (2.71) is obtained from graphs 
of the form shown on Fig. 2-10. There are four regimes identified. The 
fourth one indicates the hysteresis effect shown more clearly on Fig. 2-11. 
To use the graphs, the anticipated leak rate, either gas (volume) or liquid 
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Fig. 2-10 - Relation Between Conductance Parametei 
and Interface Contact Stress {Ref. 31) 
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(Leak- rate of the 
Connection) 



(Interface Total 
Force of the 
Connection) 


Fig. 2-11 - Hysteresis Effect (Ref. 31) 
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(weight) leak rate, is the basis for computing the required conductance 
parameter. For laminar, isothermal, compressible (gas) flow the volume 
leak rate is 


Q 


, 2 2 , 

w(p 2 -Pi) 3 

24/xp o L 


(2.72) 


where w is the width and L is the length of the leak path, ju. is the viscosity 
of the medium and p^ the standard atmospheric pressure. The inlet pressure 
is p^ and the exit pressure is pj. Other effects such as inertia, transition 
flow with molecular correction and adiabatic frictionless flow are presented 
in Ref. 31. These are, however, less important than the one given by Eq. 
(2.72). The volume leak rate Q of a gas can be converted into a mass leak 
rate by the relation 

W = rY Q (2.73) 


where R is the gas constant, T is the absolute temperature, p is the pressure 
and m is the molecular weight. For laminar, incompressible (liquid) flow 
the mass leak rate is 


W 


Pw(p 2 -pj) 

12 p L 


(2.74) 


where p is the mass density of the medium. The width and the length of the 
leak path for a gasket are, respectively 


and 


w ~ Zn i-q 
L ^ b G • 


(2.75) 


Usually the definition of zero-leak is defined in terms of volume 
per unit time of helium at standard temperature and pressure. This leak 
rate can be converted into an equivalent liquid leak rate by using coversion 
graphs. One such graph is described in a report by Weiner (Ref. 32), which 
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represents the Poiseuille equation for gas and liquid flow, as given by Eqs. 
(2.72) and (2.74). The procedure is illustrated on Fig. 2-12. 

The design procedure for finding the width of a flat gasket is based 
on the condition that 


This leads to 


p(l) _ P (2) _ p , p (2) 
G " F " p G 


Zn b iff r G ct G = * r G P + 2,r b i« r G k p P 


which can be solved for after substituting 


b eff r i b G ’ 
b eff = 7 2 b G ’ 


(2.76) 


(2.77) 


(2.78) 

(2.79) 


resulting in 


. r G P 

G " 2 <*1 C G ~ 7 Z k P P) 


(2.80) 


The safety factor should be attached to P^ , taking it as (F.S.) = 1 for the proof 
pressure condition and (F.S.) = 1.5 for the operating pressure condition. 


As a numerical example consider ALLPAX flat gaskets of thickness 1 /8 ; 

l/l6 and 1/32 inch, having a yield strength of = 10.0 ksi. The minimum 

stresses for precompression based on experience in conventional applications 

are o„ = 1.6; 3.7; 6.5 ksi and the slope of the straight line as shown on Fig. 

G 

2-9 is k = 2.0; 2.75; 3.5. If the gaskets are precompressed to yield stress 

P 

then the following gasket widths are obtained, as given in Table 2-9. 
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Table 2-9 

COMPARISON OF GASKET WIDTHS 


Pressure 
P (psi) 

Gasket Thickness 

h G 

1/8 

1/16 

1/32 

100 

0.0102 r 

G 

0.0103 r _ 

0.0104 r Q 

1000 

0.125 r G 

0.138 r G 

0.154 r G 


Since it is not necessary to precompress to the yield stress an alternate 
procedure would be to start with the available gasket width after the design 
has progressed to this point. It is 

b = r - h ^ le - r - 2Ar (2.81) 

avail 

A tolerance of Ar is provided in this formula. Then from Eq. (2.77) the re- 
quired flange force under operating condition is 

p ^ 2) = ir + 2* b Q y 2 r Q k p (F.S.) (2.82) 


where (F.S.) is the appropriate factor of safety. The condition of Eq. (2.76), 
making the initial flange force equal to the one under operating condition, 
gives the required initial contact stress as 

o(2) 


G 


2;r b. 


avail 


^1 r G 


(2.83) 


If is less than the minimum precompression stress required to seat the 

gasket, the initial flange load should be increased to achieve this minimum 

precompression stress. For example, if for a 1/10 inch ALLPAX gasket cr^, 

as computed with Eq. (2.83) is less than 3.7 ksi, the initial flange load should 

be P^ = 2jr b^, ^1 r G ksi). 

° avail 
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The safety factor in Eq. (2.82) will compensate for reduction of gasket 
stress due to the elastic deformation of the connection. The analysis of this 
deformation is described in Section 3. 

2.5 PRESSURE ENERGIZED SEAL 

The application of a pressure energized seal in both a cantilever flange 
and a flange with metal-to-metal contact is illustrated on Figs. 2-13 and 2-14, 
respectively. 

The basic difference of the two flange configurations can be seen in the 
accompanying calculations of bolt forces. The bolt force required for the 
cantilever flange is simply 


B = < F - S -» P p 

(2.84) 

P p ='6 

(2.85) 


where 


For the metal-to-metal flange the pivot point (X) is outside the bolt circle, 
while previously it was in line with the middle surface of the tube wall. Taking 
the bending moment about point (X) , i.e., assuming a situation where differ- 
ential axial motion exists at the seal-to-flange interface, the required bolt 
force is approximately 


P B = (F.S.)P p (l +”) - (2.86) 

In any case the required bolt force is by the factor of (1 + e/e^) higher than 
the corresponding cantilever flange. 


The size of the seal gland depends on the type and size of seal to be 

used. These dimensions, h and b , are supplied by the seal manufacturers' 

s s 

catalogs. The height of the recess, h^, for the cantilever flange is to be 
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Fig. 2-14 - Metal -to -Metal Flange with Pressure 
Energized Seal 
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determined, from the roll angle and corresponding differential axial displace- 
ment at the outer edge of the flange cross section. If the recess is not high 
enough, this will result in the same situation as for the metal-to-metal flange 
and is therefore undesirable. 


The width b_, carrying the same label for convenience of notation as 
G 

the width of a flat gasket, is assumed as being equal to the tube wall thickness, 
t. The distance of the seal gland from the bolt holes was assumed as being 
the same as the width of the seal gland itself. These dimensional relations 
cannot be readily defined and will have to be determined by a developmental 
test program. 

2.6 BOLT FORCE, NUMBER OF BOLTS AND BOLT SPACING 


The required bolt force is determined by the gasket initial stress and 
minimum stress during operation, or in the case of a pressure energized 
seal, by the force required to prevent separation near the seal-flange inter- 
face. The maximum force of the ones determined by different criteria is 
used to compute the number of bolts required. Two design criteria are used. 
Under proof pressure the bolts should not yield and under burst pressure they 
should not break. These two criteria can be formulated as 


n 


B1 



(2.87) 


where P_ is the bolt force under proof pressure, including the safety factor, 
B 

and 


n 


B2 


^ (burst) 

f' b > a b 

tu oB 


( 2 . 88 ) 


where P R is the bolt force at burst pressure. The minimum number 

(burst) 

of bolts was assumed as six. This will give an even stress distribution for 
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flanges with low internal pressures and small inner diameter, for which 
less than six bolts would be computed according to Eqs. (2.87) and (2.88). 


The bolt spacing is simply 


2jt r. 


s = 


n 


(2.89) 


B 


where is the maximum of the numbers of bolts, or n j$2’ com P u t e( * 

previously. The spacing should not increase beyond a certain level, which 
has been arbitrarily fixed at s = 8 dg, where dg is the nominal bolt diameter. 
This maximum spacing depends on the thickness of the flange, too, since 
bending out of the plane of the flange would introduce a reduction in interface 
stress in the space between the bolt holes. For the low profile flanges, how- 
ever, this situation is not critical since the aspect ratio of h/b for the flange 
cross section is usually greater than 3/4, mostly being around 1. Therefore 
it acts quite differently from a flat plate assumed in previous flange design 
methods. A minimum spacing is provided by the value of s = r) Q dg, where 
n o is tabulated for various types of bolt heads as a function of nominal 
diameter . 


2.7 FLANGE HEIGHT 


The computation of the flange height is based on the capacity to carry 
the ultimate applied moment 


(F.S.) Pg e 
m Fu = 2 tt r 

o 


(2.90) 


where Pg is the maximum bolt force considered for the design and e the 
internal lever arm between the bolt circle and the gasket circle, 


e 


(2.91) 
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The radius r is the one of the middle surface of the tube wall, 
o 

r = r. + t/2 (2.92) 

The width of the flange cross section has already been determined, either 
from considerations to accommodate the boltheads or to accommodate the 
gasket. The effect of the bolt holes, however, has to be taken into account. 

A simple rule has been suggested by Schwaigerer (Ref. 2) based on experience, 
by computing an effective width, b, from the bolt hole diameter d^ole an< ^ 
the bolt spacing, s, 

5 ■ » - <2 ' 93) 

Previous design methods have suggested to subtract the entire hole diameter. 
This would be unduly conservative as proven by tests (Ref. 3), 

The computation of the flange height assumes a linear stress distribution 
in the flange and the development of a plastic hinge in the neck. This procedure 
of designing a statically indeterminate structure by introducing plastic hinges 
to reduce redundancies was first used for steel frames (Ref. 33) and resulted 
in more efficient designs. The state of stress in the neck of the flange is 
three-dimensional, however, and the method used in frame design is, there- 
fore, not rigorously applicable. 

The derivation of the concept of a plastic section modulus, Z^,, for the 
tube wall is described in detail in Section 3. The result is 

t 2 . t 2 

Z T * h— I 2 - 94 ' 

where and are coefficients determined by the state of stress in the 
flange neck. Since this state of stress is unknown at this point they are 
assumed to be 
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= .8, £ 2 = .18 (2.95) 

In Section 3 the computation of and C 2 for a given flange under given loading 
conditions is shown in detail. 


The elastic section modulus of the flange cross section is given by 


b h Z 

’f " 6 r 


{2.96) 


The design formula for the flange is now derived by requiring 

tt(T) 


m _ f (F) 

m Fu - F ty 


J3L 


S F + (F) ~T 

_ F * 


(2.97) 


Usually the flange and tube wall are made of the same material so that 

F (T) /F (F) = 1. When the expressions for Z_ and S_ are substituted into 
ty ' ty 1 r 

Eq. (2.97) a quadratic equation for h results, 

Ah 2 + Bh + C = 0 , (2.98) 


where 

A = F^> b/6 r o , 

B = t 2 (t-t N )/2, 

C = V 2 - 1 n)/ 4 - ™F U - 

The solution for h is 

, Vb 2 - 4AC - B 

h — a 


(2.99) 

(2.100) 
(2.101) 

( 2 . 102 ) 
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If the contribution of the plastic hinge in the neck is neglected in the 
design of the flange, i.e., when and are assumed to be zero, then 

h = V 6 r o "W F £ ’ * • (2 ' l03) 

This design formula has been used previously for computational convenience 
but may result in overly conservative designs. 

Finally, a check is made of the flange height versus the bolt spacing, e, 

if s/h > 3 h = s/3 (2.104) 

2.8 FLANGE WEIGHT 

The weight added to the tube by the flange is given by computing the 
volume of the material having the cross sectional area 

= (b - t) h (2.105) 

and the centroidal radius 

r w = r.+(t+b)/2 (2.106) 

so that 

voi = 2j r r w A w (2.107) 

the actual weight is 

AW = p F vol (2.108) 
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2.9 MATERIAL DATA 

To facilitate the computation of numerical examples the properties for 
aluminum and steel commonly used for rocket propulsion systems are given 
in Tables 2-10 and 2-12. These data were taken from Ref. 15. Data for 
gaskets were compiled for some materials used in some earlier MSFC com- 
putations (Ref. 34) and are listed in Tables 2-11 and 2-13, 

Both data tables are incorporated in the computer program. They can 
be enlarged easily by including a larger variety of data. It was not the purpose 
of this study to compile all available data. 
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Table 2-10 

PROPERTIES OF METALIC MATERIALS FOR TUBES, FLANGES AND BOLTS 



Material 

E 

(psi) 

V 

(-) 

(lb/in. 3 ) 

a 

(in./in./°F) 

F f 

ty 

(psi) 

F tu 

(psi) 

1 

Af 6061-T6 @ RT 

9.9 x 10 6 

.33 

.098 

12.5 x 10" 6 

35.0 x 10 3 

42.2 x 10 3 

2 

Af 6061-T6 @ 200°F 

9.9 x 10 6 

.33 

.098 

12.5 x 10" 6 

32.2 x 10 3 

38.1 x 10 3 

3 

Af 2024-T3 @ RT 

9.9 x 10 6 

.33 

.098 

12.5 x 10" 6 

50.0 x 10 3 

62.0 x 10 3 

4 

Af 2024-T3 @ 200°F 

9.9 x 10 6 

.33 

.098 

12.5 x 10“ 6 

47.0 x 10 3 

59.0 x 10 3 

5 

347 SS @ RT 

28.0 x 10 6 

.30 

.288 

9.5 x 10" 6 

35.0 x 10 3 

90.0 x 10 3 

6 

347 SS @i 200°F 

28.0 x 10 6 

.30 

.288 

9.5 x 10‘ 6 

30.0 x 10 3 

76.0 x 10 3 

7 

347 SS 600°F 

28.0 x 10 6 

.30 

.288 

9.5 x 10' 6 

25.0 x 10 3 

68.0 x 10 3 

8 

A286 @ RT 

to 

00 

o 

X 

o 

CT' 

.30 

.288 

9.5 x 10" 6 

131.0 x 10 3 

200.0 x 10 3 

9 

A286 @ 200°F 

28.0 x 10 6 

.30 

.288 

9.5 x 10 -6 

128.0 x 10 3 

196.0 x 10 3 

10 

A286 @ 600°F 

28.0 x 10 6 

.30 

.288 

9.5 x 10“ 6 

120.0 x 10 3 

180.0 x 10 3 


Legend: E = elastic modulus 
V - Poisson 1 s ratio 
p = weight density 


a 



= linear thermal expansion coefficient 
= tensile yield strength 
= ultimate tensile strength 


LMSC-HREC TR D306492 




LOCKHEED HUNTSVILLE RESEARCH 4 ENGINEERING CENTER 


Table 2-11 

PROPERTIES OF GASKET MATERIALS 


M 

i 

tfc. 


No. 

Material 

E 

(psi) 

k g 

(psi) 

a 

G 

(psi) 

a 

(in/ in/°F ) 

M 

(-) 

h G 

(in.) 

k 

P 

1 

Asbestos l/32 in. 

44.0 x 10 3 

10.0 x 10 3 

6.5 x 10 3 

1.3 x 10' 3 

.5 

.03125 

3.50 

2 

Asbestos l/l6 in. 

44.0 x 10 3 

10.0 x 10 3 

3.7 x 10 3 

1.3 x 10" 3 

.5 

.06250 

2.75 

3 

Asbestos l/8 in. 

44.0 x 10 3 

10.0 x 10 3 

1.6 x 10 3 

1.3 x 10 -3 

.5 

.12500 

2.00 

4 

KEL-F81 

180.0 x 10 3 

8.0 x 10 3 

4.0 x 10 3 

3.8 x 10" 5 

.12 

.062 50 

3.00 

5 

ORES 321-A 

28.0 x 10 6 

40.0 x 10 3 
1 

18.9 x 10 3 

9.5 x 10~ 6 

.30 

1 

.02500 

5.50 


Legend: 


Assumed: 


E = elastic modulus = y 2 = 0.5 for No's 1,2 and 3 

K g = yield (crushing) strength ^ = y 2 = 1.0 for No's 4 and 5 

cr„ = minimum seating stress 
G 

a = linear thermal expansion coefficient 

jl = friction coefficient 

h,-. = thickness of the gasket 
G 

k = ratio of required seating stress to given pressure (see Fig. 2-9). 
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METRIC PROPERTIES OF 





E 


Material 


(N/mm 2 ) 

1 

At 6061-T6 

RT 

68 x 10 3 

2 

A l 6061-T6 

93 C 

68 x 10 3 


At 2024-T3 

RT 

68 x 10 3 


A£ 2024-T3 

93 C 

68 x 10^ 

5 

347 SS 

RT 

193 x 10 3 

6 

347 SS 

93C 

193 x iO 3 

7 

347 SS 

31 5C 

193 x 10 3 

8 

A286 

RT 

193 x 10 3 


A286 

93 C 

193 x 10 3 


A286 

315C 

193 x 10 3 


Table 2-12 


JLIC MATERIALS FOR TUBES, FLANGES AND BOLTS 


V 

(-) 

P 

(g/ram 3 ) 

a 

mm/ mm/°C 

F *Y 2 

{N/mm ) 

F. 

tu ? 
(N/mm ) 

.33 

.271 x 10' 2 

22.5 x 10“ 6 

242 

290 

.33 

.271 x IO -2 

22.5 x 10" 6 

222 

263 

.33 

.271 x 10" 2 

22.5 x 10' 6 

345 

428 

.33 

.271 x 10" 2 

22.5 x 10” 6 

324 

407 

.30 

.798 x 10" 2 

17.1 x IO -6 

242 

621 

.30 

.798 x 10“ 2 

17.1 x 10' 6 

207 

528 

.30 

.798 x 10" 2 

17.1 x IO -6 

173 

469 

.30 

.798 x IO -2 

17.1 x IO" 6 

904 

1380 

.30 

.798 x 10~ 2 

17.1 x IO" 6 

883 

1352 

,30 

.798 x 10‘ 2 

17.1 x IO" 6 

828 

1242 
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Section 3 

ANALYSIS METHOD 


The analysis method described in this section is based on thin shell 
theory and simple ring theory. These theories are not too involved algebra- 
ically to be used for hand computations. Also the approximate state of stress 
in the plastic hinge near the flange used is described. A summary of the 
formulas used in the analysis is given in Appendix B. 

3.1 SHELL THEORY 


The membrane solution for a cylindrical shell under an internal pressure 
p and a temperature differential AT is characterized by the stress resultants. 


n 

x 


p r 
r o 

2 


(3.1) 


and 


> = P r o {3 - 2) 

where n^, are the axial and circumferential stress resultants, respectively 

measured as a force per unit length. The radial expansion of the shell under 
this loading condition is 


W = 



(‘-T * 4 r o aAT 


(3.3) 


where a is the linear thermal expansion coefficient and E and v are the 
elastic modulus and Poisson's ratio, respectively. 
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In addition to this solution the edge disturbance of the cylinder, introduced 
by the flange, must be considered. It can be shown (Ref. 35) that the linear dif- 
ferential equation 


,4 

d w 
dx 4 


k 4 w = 0 


(3.4) 


where 


4 12 ( 1 - v Z ) 

1 = 2 2 
r o t 


(3.5) 


describes this behavior. This differential equation for the range of parameters 
considered, assuming the shell to be infinitely long, has the solution 

w = e k* (Cj coskx + C2 sinkx) (3.6) 


The integration constants are found from the edge conditions. The flange 
usually introduces an edge moment m Q and an edge shear q Q into the shell. 
Both are measured per unit length (Fig. 3-1). Knowing that 



(3.7) 


where 



= -B 


d 3 w 

dx 3 


x = 0 


(3.8) 


B = 


Et" 


12(l-I' fi ') 


(3.9) 
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£ 


Fig. 3-1 - Edge-Loaded Cylindrical Shell 
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Tho constants are derived using 


dw _ . -kx 
dx " ke 


(Ci -G,) cos kx + (Cj + C^) sinkx 


— = 2k^ e ^ (Cj sinkx - Ck cos kx), 

dx 


^ = -2 k 3 e ^ |(Cj - C,) sinkx - (C^ + C^) cos kx j 


It follows then that 


c i ■ 


m q km 

o . r _ o o 

2 ’ - 3 

2 k B c 2k J B 


With these constants the radial displacement is 


w 


1 

3 

2k B 


-kx 

e 


q coskx - km (coskx - sinkx) 
o o 


and the rotation (rolling) of the shell wall is 

X = “TT = — e ^ I -q„ (cos kx + sinkx) + 2 km coskx 

dx 2k c B [ o o 


For the edge where x=0 the flexibility matrix is seen to be 


w 

1 

*1 -k " 


^0 

X_ 

2k 3 B 

-k 2k 2 


m 

L °i 
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The meridional bending moment along the shell wall is 


m 


= e k* £ m 0 (coskx + sinkx) - sinkxj (3.17) 


and the meridional shear is 


= e ^ £ q o (sinkx - coskx) + 2k sinkxj 


(3.18) 


The circumferential bending moment is 


m.(p -V m 


(3.19) 


and the circumferential stress resultant is 


Et 

n (p = — w 
o 


(3.20) 


This concludes the description of the analysis of the edge disturbance. 


It remains to be shown how the stresses are computed in the elastic range 
and how the plastic state of stress is described. Three stresses exist in the 
shell, the axial stress 



12 


(3.21) 


the circumferential stress 


V 


JE t . m <P 

* £ 

12 


(3.22) 
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and the shear stress 


r 

xz 




(3.23) 


The coordinate z is measured from the shell middle surface outward in the 
normal direction. 


To arrive at an expression for the development of a plastic hinge in the 
shell it is assumed that a core (Fig* 3-2) of thickness, 



(3.24) 


is required to carry the axial force, where Yo is the uniaxial tensile yield 
strength of the material. This leaves for the plastic moment, , 


rP - 


mT 

? )/4 


and the plastic shear force, , 


P _ 


xz 




(3.25) 


(3.26) 


In order to relate the three-dimensional state of stress to the uniaxial tensile 
yield strength Y q the yield condition of von Mises is used. 




a 2 ) 2 + 




- ° 3 > 


(<r 3 - aj)' 


< Y 
— o 


(3.27) 


where cr^, are the principal stresses. 
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Fig. 3-2 - Assumed Stress Dis 
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The principal stresses for the problem at hand are 



(3.28) 


*2 




(3.29) 



(3.30) 


The expressions for the principal stresses are simplified by introducing 

ffP as a reference stress, where 
x 


cr 


P 

<P 



and 


T P 
xz 



(3.31) 


(3.32) 


then 


CT 


1 . 2 




(3.33) 



1 


a? 

x 


(3.34) 
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and the equivalent stress a of Eq. {3.27) is 


— P 
cr = a 


f 


1 + Of _ + 


a l + 


3a; 


(3.35) 


the computation of a^, and a for a given loading condition will be shown 
later. 


3.2 FLANGE THEORY 


Adding the flange to the shell requires finding the interface moment m Q 
and interface shear q Q (Fig. 3-3) in terms of given loading conditions. 


In the analysis of the flange deformations it is useful to derive an equiv- 
alent rotational spring constant per unit length of the flange (Ref. 36). Starting 
with the equation for the radial displacement w and rotation x of the inter- 
face point (A) (Fig. 3-3) caused by an applied moment m^ per unit length, 


w 


2 , I 






r r 

= - o c 

c + — c 

A 


q 

o 

0 

o 

f 


EX 




T El 


x_ 


c 1 


m 


1_ 


(3.36) 


where A is the cross sectional area and I is the moment of inertia of the ring 
cross section, an equation for m Q and q^ can be constructed by requiring 
compatibility of the displacements of point (A) on the ring and on the shell. 
Using Eq. (3.16) for the shell displacements it follows that 


1 

‘ 1 -k' 


V 

3 

2k B 

-k 2k 2 


m 



„ o 


r r 
o c 


El 




+ 


r r 
o c 


El 



m 


F 


1 


(3.37) 
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Fig. 3-3 - Ring-Shell Interface 
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which can be combined as 



+ 3 <c 2 + £) 



cp 


“1 


“ 

2k + cf3 


% 

1 + 0 


m 

o 



where 


P 


B k r r 
o c 

El 


The determinant of this equation is 


D = (1+0) 


1 , a , 2 . I . 

2 + P ( c + x) 

2k A 


-(-i+cP ) 2 


To find q and m Cramer's rule is used, 
o o 


= A. 


D * C + 2k ) m F 


and 


m o = 




+ 



The rotation of the cross section is then 
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(3.41) 


(3.42) 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306492 


X = — | (2k 2 m - kq ) 

2k (3 


hr' 2 * 2 *** I ) m F 


2k BD 


(3.43) 


The rotational spring constant is obtained from Eq. (3,43) by dividing mp by 
the rotation x> 



BD 


P 


ILL 

K A 



(3.44) 


The second loading condition to be considered in this paragraph is a 
differential radial displacement Aw between the ring and the shell where 


Aw 


W shell 


- w . 
ring 


(3.45) 


the corresponding equation for m Q and q to be solved is obtained by re 
placing the right hand side of Eq. (3.38) by 


PEI 

r r 
o c 


I 



Aw . 


(3-46) 


The solution is given by 


, =--m_ 

n o r r D 
o c 


<1 + (3 ) Aw , 


(3.47) 
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m 


= - (ik - cP ) Aw 

o c ' ' 


(3.48) 


and the rotation is in accordance with Eq. (3.43). 


X = "d < c + Aw 


(3.49) 


The same rotation can "be produced by an applied moment of 


m F = C F* 


E ( c + ?k> 

(k s + jt) 


Aw 


(3.50) 


The stresses in the ring at point s(i^ and (B) are 


a 


(P 


w 


(3.51) 


and 


<P 


(w - hx ) 


(3.52) 


where w and X are the radial displacement and rotation, respectively, of 
point (^) 

3.3 EFFECTS OF BOLTS AND GASKET 

The bolts and the gasket contribute to the elastic properties of the flanged 
connection. Both can be thought of as elastic springs (Fig. 3-4) whose spring 
constants can be combined with the equivalent spring of the flange. The gasket 
spring constant, k^, is 
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where the gasket area, A_, is 

G 


a g e g 

2 ' r o l G 


A G = 2,rr G b G 


(3.53) 


(3.54) 


E G is the elastic modulus and t G is the thickness of the gasket. Similarly, 
the bolt spring constant is 


k 


B ~ 


V 

2 jr r 


B_ 

l 


o 


B 


(3.55) 


where the total bolt area, A^, for n„ bolts is 

B B 


A 


B 


n B A oB 


(3.56) 


* s the stress area of a single bolt, and is the stress portion of the 
bolt shaft. With the radial distance 


a 


= r n 


- r. 


(3,57) 


between the bolt circle and the gasket the equivalent rotational spring is 


k G k B 
k G + k B 


The centroid of both springs is given by the radius 


r 

a 


k B 

r G + T-n^r a 


B 


k 



a 


(3.58) 


(3.59) 
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Finally, the displacements in gasket and bolts are, respectively, 


K 


6 G = u + x 


B 


K G + K G 


and 


6 b =u ' x 


k g 

k b + k g 


(3.60) 


(3.61) 


where u is the axial displacement at the centroid of both springs, and the 
corresponding changes in gasket and bolt stresses are, respectively, 


A E G 6 G 

Act g " t_ 


(3.62) 


and 


<5 

Aa = -BB 
B 


(3.63) 


3.4 SEQUENCE OF LOADING CONDITIONS 


In the preceding three paragraphs the mathematical apparatus for the 
analysis of the deformations and stresses of a flange have been presented. 
It will now be used in the step-by-step analysis of the loading conditions. 

The initial loading of the flange occurs when the bolts are torqued to 
achieve a tight seat of the gasket. This force was computed in Section 2 
based on the gasket design requirements. This initial bolt force may be 
related to the bolt torque applied when the connection is assembled, for 
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which torquing charts are available. It would go beyond the scope of this 
report to go into these torquing requirements. For the further discussion 
a bolt force, f^, per unit length of circle r , 


{ o) * 

,<o> = B^B A °B 
B 2 tt r 

o 


(3.64) 


is considered, where o is the stress in the bolts at initial torquing. For 

15 

a cantilever flange the corresponding applied flange moment is 



a f 


(o) 

B 


(3.65) 


It is evident that from Eq. (3.43) the rotation X of the cross section and from 
Eqs. (3.41) and (3.42) the interface shear q Q and interface moment m Q can be 
computed and the remainder of the analysis of shell and flange be performed 
as described in paragraphs 3.1 and 3.2. 


When the separation of the fluid system is started the internal fluid 
pressure causes an axial force in the tube 



and a force 


f 


F 



(3.66) 


(3.67) 


on the face of the flange. 


The latter force acts at a radius 

2 2 
9 r , r. r , r. 

_ 2 G+iG+i 

F " 3 r_ + r. 

G l 


(3.68) 
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The corresponding applied flange moment is 

m F > = f T (r a ' r Q > + f F (r a " r F } (3,69) 

121 

Another flange moment nip ' is caused by the differential radial displacement, 
according to Eq. (3.45) and Eq. (3.50). The radial displacement of the shell 
for the most general case was given by Eq. (3.3). The term attributed to the 
temperature differential AT is probably unrealistic for cryogenic applica- 
tions when assumed that the ring could not experience the same differential, 
i.e., both ring and shell probably experience simultaneously the same AT and 
therefore this term does not produce a Aw. The radial expansion of the flange 
ring due to internal pressure is 



(3.70) 


It is now possible to compute the rotations of the flange. Initially a 
rotation 



(3.71) 


occurs. The corresponding axial displacement is the reference position and 
taken as 


u<°> . 0 . 


(3.72) 


When m^ and are applied the rotation is 


(1) (2) 

(p) _ m F ± m F 

C E + C F 


(3.73) 
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and the axial displacement is 


(P) _ f T+ f F 
_K G +K B 


(3.74) 


The gasket and bolt deformations according to Eq. (3.60) and (3.61) are 
evaluated with and X^ . The final stresses in the flange, however, 

are computed with 


X (T) =x (o) + x (p) 


(3.75) 


for which a corresponding can be computed with Eq. (3.49). It is not 

necessary to repeat here how the stresses in the shell and the flange are 

computed from the moment m^J^ and the rotation X* T *. In summary, an 

(T) * rn 

interface moment m' and an interface shear o' ' are arrived at. Also a 

o o 

radial displacement at point (^A) (Fig. 3-3) of is computed. 


The plastic stresses at the flange neck are generated by increasing the 
(T) (T) D n 

pressure until m l and q v become m p and q p as in Eqs. (3.25) and (3.26). 

° / T V O X / r T 1 \ ^ 

At the same time \ y ’ increases to X P and w l 'increases to w P . The 

o 

stresses are then 



(3.76) 


(3.77) 
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so that 


and 


Cr 


P 

<P 


E_ 

r 

o 


w^ + V <jP 
X 




a. 


t + t 

n 

4 


c + 


2k 


2 k A 


c 

2k 


(3.78) 


(3.79) 


(3.80) 


according to Eqs. (3.31) and (3.32) . 

3.5 ESTIMATE OF THE MOMENT CAPACITY OF THE FLANGE 

The capacity of the flange to carry an applied moment of m_ is exhausted 

X 

when 


”Vu * Y o [ Z F + Z t] ( 3 - 81 > 

where Y as the tensile yield strength of the material and Z_ and Z_, are 
o F T 

the equivalent plastic section moduli of the flange ring and the tube, 

respectively. This is the same equation as Eq. (2.97). A more conservative 

assumption would be to let the stresses in the ring just reach the yield stress 

in the extreme fibers so that the elastic modulus S„ should be used instead 

IT 

of Zp. For a rectangular ring cross section with the reduced width b the 
two section moduli are 


bhi 

F "4r 

o 
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and 



(3.83) 


The equivalent plastic section modulus of the tube wall, Z-p can be expressed 
in terms of the expressions derived in Eqs. (3.76) through (3.80) and Eqs. (3.31), 
(3.32) and (3.35) when 


cr = Y. 


(3.84) 


then 


Z x 



V 1 + °i 


a 2 <* - ! 


+ a z x + 3a^ 


(3.85) 


or simply 

2 2 

Z T = & 1 4~ ^ + ^ 2 (t ' V 4 

The two dimensionless parameters are 


and 


£i = 


^ 2 “ 


\j l + a 1 + a ^ + 

3 °2 

a 2 



V 1 + a \ + 'a j + 3a^ 


(3.86) 


(3.87) 


(3.88) 
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Section 4 

COMPUTER PROGRAM 

This section describes the computer program which was developed to 
implement the design standard and verify the stresses and deformations of 
the flange. The program is written in FORTRAN IV language for use on the 
Univac 1108 Exec 8 system. The algorithms of these computer programs 
are based on the design procedure and the analysis method outlined in the 
previous two sections, A listing of the code is included. Input instructions 
for the computer program are given in Appendix C. Example problems are 
presented in Section 5. 

• PROGRAM OUTLINE 

The program consists of a main program which reads the input data, 
and four major subroutines in addition to two output routines. These major 
routines are DESIGN and ANALYS, corresponding to the design and analysis 
part of the program and PLOTF1 and PLOTF2, which are the two plot routines 
for the Stromberg-Carlson 4020 plotter. The organization of the entire pro- 
gram is shown on Chart 4-1 and the individual routines are briefly described 
in Table 4-1. The two routines DESIGN and ANALYS follow principally the 
sequence of formulas given in Appendixes A and B. The individual variables 
are easily recognizable and are therefore not explained here in detail. 

The program allows the design and analysis of cantilever flanges with 
flat gaskets and pressure energized seals. The machining of the upper flange 
surface may be with machined spot faces or with a machined grove. These 
different options can be turned on by specifying the appropriate values of the 
variable KOPT(I), as described in the User's instructions in Appendix C. 
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Table 4-1 

PROGRAM DESCRIPTION 


No. 

Symbol 

Name 

Description 

1 

DMAIN 


Design program for low profile flanges 

2 

D02 1 

METALS 

Table of metallic materials design properties 
for tubes, flanges and bolts 

3 

D022 

GASKET 

Table of gasket materials design properties 

4 

D001 

DESIGN 

Design routine for low profile flanges 

5 

D010 

BOLT 

Bolt data handling 

6 

DO 11 

BTABL1 

Bolt table for open wrenching 

7 

DO 12 

BTABL2 

Bolt table for socket wrenching 

8 

D013 

BTABL3 

Bolt table for internal wrenching 

9 

P001 

PLOTF1 

Plot routine for low profile flanges with flat 
gasket and machined spotfaces for the holes 

10 

P002 

PLOTHC 

Plot a half circle from IA to IB 

11 

P003 

PLOTLN 

Plot a line 

12 

P004 

PLOTLB 

Plot label 

13 

P005 

PLOTQC 

Plot a quarter circle from IA to IB 

14 

P006 

PLOTAR 

Plot an arrow head for different orientations 

15 

P007 

PLOTTX 

Plot text 

16 

P008 

DASHLN 

Dashed -dotted line 

17 

D100 

OUTDES 

Output of the design routine 

18 

A001 

analys 

Analysis routine 

19 

P010 

PLOTF2 

Plot of the analysis results 

20 

D200 

OUTAN 

Output of the analysis results 
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The plot routines PLOTF1 and PLOTF2 summarize the design and 
analysis. The first one plots the geometry of the flange cross section in 1:1 
scale. The second one summarizes the stresses and deformations of the tube 
wall and the flange, using a 1:2 scale for the flange geometry. The layout of 
the graphs is given on Figs. 4-1 and 4-2. The small numbers refer to x and 
y coordinate points in the code and are given here to facilitate future modifi- 
cations in the program. 

A list of the entire code is given in this section. The limited scope of 
this contract did not allow inclusion of all possible flange configurations to 
be considered in this program with the corresponding plot option. At this 
point, however, it would be possible to automate the design process further 
by combining the computer code with a different type of plotting equipment, 
allowing larger size plots. The SC 4020 plot area is limited to l\ by ~l\ inch. 

Sample computer output, printed and plotted, is presented in Section 5. 
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Fig. 4-1 - Layout of Design Summary SC 4020 Plot 
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PROGRAM LISTING 
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LIST UF ROOTlNfcS IN FLANGE DESIGN AND ANALYSIS PR u G* A M 


‘♦3‘*6Co*TPF»»LIST 1 


1 

fcHDG LIST Of 

ROOT 1 N £. S IN FLANGE DESIGN AND ANALYSIS PROGRAM 

2 

foPRT*L LI SI 1 



3 

muHDG DMA IN 

( 1 ) 

(MAIN PRQGM) 


wPKT.C DMA IN 



b 

WrtDG D u 2 1 

( 2 ) 

( metals ) 

6 

toPKT.C DL2 1 



/ 

ia li D G 0 C 2 2 

(3) 

(GASKET ) 

a 

GPRT.C Do22 



V 

fa HOG Dll 1 

( H ) 

(DESIGN! 

lu 

(aPKT.L DDLl 



1 1 

WHDG DC 10 

( S ) 

(bOLT) 

U 

la P R T t L D u 1 u 



13 

iartOG DC I 1 

( 6 ) 

( bT A6L 1 ) 

1H 

ISPKTit DLll 



lb 

WHDG DC 1 2 

( 7 > 

IBTA0L2 ) 

16 

£PKT,C DC 1 2 



1 7 

fa H D G DC 1 3 

( o ) 

{ U T ABL 3 ) 

16 

faPHT.C DC 1 3 



1 9 

laHDG P C C 1 

1 9 1 

( PLU^F 1 ) 

20 

f*PRT»C PCD J 



21 

iibHDG P C C 2 

(ID) 

(PLOThC ) 

22 

laPRTtC PwC2 



23 

ft'HDG PCC 3 

(11) 

( PLOTLN > 

2S 

(a P k T » L P o D 3 



2& 

lainOG PGCiH 

( 1 2 > 

( PIGTlB ) 

26 

(o-PRT * L PCcH 



27 

fa HOG PQCS 

< J 3) 

(PLOTwC ) 

26 

GPRT.C PloS 



29 

fehuG PuL6 

()k) 

( PLOTmR ) 

30 

laPRT.t Pll6 



31 

ghdg p d c 7 

(is ) 

(PLOTTX ) 

32 

ttPKT.L PCD 7 



33 

(»i H D G P u C P 

(16) 

(CiaSHlN ) 

3 ^ 

tPRT.L PCCB 



3b 

wHDG D 1 ul 

(17) 

(OUIDES) 

36 

laPKTiC 0 1 Co 



37 

laHDG A Jo I 

(IB) 

C A N A L Y S 1 

36 

laPR T ♦ L A w o 1 



3V 

i»< 1 1 D G P & i l 

(19) 

( PLOTF 2 1 


taPRT.L Polo 



<tl 

laHDG D2DC 

( 2D ) 

(UUT AN 1 

H2 

fa P K T # C D2yc 




«HDG UMA1N 111 (MAIN PROGM ) 


BPRT.C DMA IN 

fuRPUK 2**H|-D3/10-X**;S3 
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‘♦3‘*fc0u*TPFS,UMAIN 


1 

2 
J 
H 
b 
6 
7 
b 
9 
Jo 
J 1 
1 2 
1 3 
J *t 

15 

16 

1 7 
lb 

19 

20 
21 
22 

2 3 
2*1 
2b 

26 

27 

2b 

29 

30 

31 

32 

33 

3 H 
3 b 

36 

37 
3b 


C 

C 

c 

c 

c 


OtSUM PROGRAM for lo» profile flanges 
K.R.LE lMbACH. LGC<HtiD-HUNTSY ILdE. tXT,3b3 
21 NOVEMBER 1972 

1 FORMAT ( 12A6) 

2 FORMAT I bt lo .HI 

3 FORMAT l l 6 I S 1 

0 I MENS l ON HEAD (12) , AD (22 ) .KOpT ( 10l 
DIMENSION A(9,*n . SRfcS 1 5 . 9 > , Si R < b »H) » AP ! b ) 

0 A 1 A l A 0 ( 1 1 » I * I i22)/22*6h / 

READ (b, j ) I AD l 1 1 , 1 -X , 1 2 ) 

CALL 1 DENT ( 9 * AD | 

KC AU ( b • J ) NCASES 


1 C ASE = l 
10 READ (b , 1 > 
READ I S ,2 > 
READ ( b * 2 ) 

RE A D ( b » 3 ) 
IF(IT.EO.O) 
IFUF.EO.u) 
IF ( 1 b • L U • u 1 
IF ( iG.EC.u) 
IF ( IG.LT.O I 


(HEAD ( I) , i = 1 , l 2 1 
PiDl i T • D E L T » HT 
PF.bF ,FS.«F 
I T , IF . IB. 1G 
READ! 5 ,2 1 
REAl»1S*2) 

read ( s 1 

RE AD ( S « 2 1 
REA0<S,2) 


ET t aNUT ,HHOT , ALFAT ,F TTT ,FTUT 

EF , AnOF » R n 0 F , ALF AF , F T Y F . F V UF 

Eb.ANDB,r<rtDB»ALFAB»FTrB.FTUe 

EG.AKG.SG.aLFAG.AMUG.GaMU.GAMS.HG.SP 

Hb.bS.hR 


IF(IT.GT.O) CALL ME T AL S ( I T , E T . ANUT , KHOT , A LF A T . F 1 Y I « F I U T 1 
IF! lF.GT .0) CALL METALS C IF ,EF . AWUF . « HOF , AcF AF , F T YF ,FTUF 1 
lFUB.GT.ul CALL METaLS( I b.EB.AMub. RHOB.ALFAB.FT Yb.F TuB ) 
IF(lG.GT.u) CALL (. A SK E T ! 1 G ,E6 » AKG ,SG , ALFAG » AMUG .GAWD ,gAMS ,Hg , bP j 

RE AD ( b . 3 ) (KOPT C I 1 . 1 = 1 . 10 1 

DIMENSION TuBmTl12).FLAMTL(2) .B0LMTH2) .GaSMTLU) 

R E AD ( b * 1 ) ( T UBMTL ( I I . 1*1 .2 ) . I FLmMT L ( ! ) .1 = 1 i 2 > 

* , (bOLMTLl 1 ) . 1" 1 .2 ) . I GASMTL 1 1 > . 1*1 .21 

H E AD ( b . 3 ) NPHASE 
R E A D ( b . 2 1 DELTAT 

iai formaTiihu 


39 

102 FORNaTI 

f 

NOMINAL PRlSSUre 

He 

* • F 1 6 , J , ♦ 

PSI */ 

*+ u 

* 

t 

NOMINAL DImMeTER 

Die 

* • F 16.3,* 

INCH*/ 

i 

• 

• 

TUBE THICKNESS 

T B 

*.F|6.3»* 

inch*/ 

Hi 

• 

« 

TUBE ThICKN TOLH 

D T« 

♦ . F 1 6 , 3 , * 

inch*/ 

N3 

• 

f 

HEIGHT TO »»ELD 

ht« 

♦ » F 1 6 . 3 , * 

INCH »/ 

HH 

) U 3 FORMAT ( 

« 

PROOF FACTOR 

PF = 

* . F 10.3/ 


Hb 

* 

1 

BURST FACTOR 

bF* 

* t F 1 u . 3 / 


46 

* 

• 

safety factor 

FSs 

* iF 10.3/ 


HI 

* 

1 

gasket factor 

6F a 

1 »F 1 u.3/// ) 

Hts 

109 FORMAT < 

1 

PROPERTIES OF tube MATERIAL*/) 


49 

105 FORMAT ( 

1 

PROPERTIES OF FlAnGf 

MATERIAL* 

/ 1 

Sl 

106 FORMAT! 

« 

PROPERTIES of bolt MATERIAL*/) 


b 1 

107 FORmaT( 

* 

PROPERTIES OF GASKET 

material • 

/ 1 

be 

108 FORMAT! 

♦ 

material Table nc 

1. 

!*♦ .15/ 


53 

* 

t 

ELASTIC MOuULUS 


E»* .E 1 6 

.8 , * P 

5 4 

* 

t 

POISSON-S kATIO 


N 0 * » .Flo 

.3/ 


ss 


OENSi r Y 


NnO** (FlDtH. * Ld/CUolC-lNLM*/ 
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(■ ORmaT ( 


I I u F 0 R n a T i 
111 FORMAT l 
1 1 ? FORMAT < 


120 FOKMATl 

1 2 1 FORMAT ( 


therm exp coeff 

alf a* » 

t E 1 6 . B . * 

INCh/IN Lh/F ♦/ 

tensile yield STR 

FT y*» 

>E 1 6 • B , * 

PSI •/ 

dlTImATE Te-NS STR 

FTU«» 

• E I 6.8 . * 

PS1 *///) 

material Table no. 

I a • 

*15/ 


elastic MODULUS 

£at 

* E 1 6 » 8 . * 

PSI ♦/ 

YIELD STRENGTH 

KG** 

»E 1 6.8 , ♦ 

PSI */ 

seating stress 

SG" * 

• E 1 8.8 , * 

PSI •/ 

THERM EXP COEFF 

ALF A* » 

• E16.8, * 

INCH/INCH/F*/ 

COEFF OF FRICTION 

HU* * 

•F 10.3/ 


WIDTH COEFF I C I ENT 

G AMU * * 

, F 1 o . 3 / 


nluTH COEFFICIENT 

GAHS* ♦ 

»F 10.3/ 


U ask e t thickness 

HG* * 

•Flu. 4.* 

INCH*/ 


SEALING 5 Tk E S 5 RATE SP *«, F I G .*»/// ) 

OPT 1 ONS • / ) 

615) 

PRESSURE AlTIVaTEO SEAL*/ 

DEPTH OF THE SEaL GlAnd HS»*»FlG.3,’ I N C H » / 

WIDTH OF TnE SEaL G l A n u b5*'*FlG.3,' INCH*/ 

DEPTH OF TmE RECESS HR* » ,F 1 0«3 * ' INCH*///) 

number of phases to be considered in The analysis 

TEMPERATURt DIFFERENTIAL "*»FIG.2 i* DEg F*/l 


computeo Thickness 


T** ,FlO«<t •• INCH 


W R I T E ( 6 • Iwl ) 

WRIT £(6.102) P.dI .T.DELT.HT 
WRITE< 6 , lii3) PF.BF.FS.SF 
w R 1 T E l 6 , I o 4 > 

WR I T £ ( 6 , lob > IT.ET , a N U T t RHOT, 
«i R 1 T E I 6 . 1 u 5 ) 

ARJTEIo.luB) IF.EF » ANUF » R HOF , 
RK I TE ( 6 , 1 u6 ) 

*RIT£(6.l0b) IB.Eb.ANUb.RHOB, 
IFUG.LT.O) DU TO 2 d 
KVRI Tt ( 6 , 1 u7 ) 

ARITE( 6 , 1 DV) IG«EG*AK 6 »SG*ALf 
00 TO 25 

20 WR 1 TE ( 6 , 1 19 ) hS . bS » hr 
HG = HR 

25 COnTInuL 

ftR 1 T E ( 6 ( l l £ ) 

WR1TE(6,11I) (KOPTC 1 ) . 1*1 , 10) 
*RITE< 6 ,) 2 u) NPhaSE .OELT AT 


lull 

102) P,ol.T»DELT,HT 
1u3) PF.BF.FS.GF 
104 ) 

10b) IT.ET, ANUT.RHOT.ALFaT.FTYT.FTUT 
1u5 ) 

10B) IF.EF, AN UF.RHOF, ALF AF.FTYF.FTUF 
1 Ob ) 

I Ob ) I6.Eb,ANUb,RH0B.ALFAB,FTYB,FTUb 
•01 00 TO 2d 
107 ) 

10V) I6.E0,AKG,SG,aLfA0,AHUG,GaMU,0AHS,hG,SP 


CALL OLS ION (P ,D I , T ,u£LT ,PF ,BF ,FS ,GF 

,ET,aNUT.RHOT,ALFaT,FTTT»FTUT 
,EF .ANuF.RhOF.AlFaF.F TYF.FTUF 
*£b , ANUB • R H 0 b ,A lF a 0 ,FTYB »FTUb 
,EG,aKG»^G.ALFAG,aHUG,GAHU,GAMS,HGiSP,HS,oS 

,KOPT , A Oo , mE IGH T ,PB 

• B.H.RI ,kG.Rb,RFIl »RSP0T, DHOLE, ObPOT ,N,Bg,hT) 
gKITE(6,12)) T 


CALL PLOTF 1 1 b , H , T . R 1 ,RG,RB,RfIL,RSPOT , DHOL E ,DSPoT ,N,bG«HG,Bb,hS,Ml 
,F TYF ,FT0F .FT YB .FTUB ,SG, AKG ,P ,*£ 1GHT 
.FLAHTL.oOLMTL.GAsMTL.HEAD.KOPI) 
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D M A I N (1J ( h A 1 N P R U G M ( 


112 C 


1 1 3 
i i h 
1 15 
1 1 6 
1 l 7 
1 1 8 
I J 9 
1 2d 
121 
122 
123 

129 
1 25 
126 
1 27 
128 
l 29 

130 

131 

132 


CAUL OUTOES1HEAD,AOb.wEI&HT,k0PT,T 

* .H.RI ,KG »KB ,RF iU »f<5P0T »DH0LE .DSPClT ,N ,bo,HT) 
C 

CALL ANaLYSIP.D[,T,L*ELT,PF.BF,FS,UF 
« , t. T t A iv U T « K H 0 T »aLFaT*F1YT , F TUT 

* .EF , ANyF tRHQF , ALF aF ,F } YF »FtUF 

* , E B » A N U B , R h 0 b , A L F a 6 if T Y 6 » F T U B 

* »EG*AKfc>.3G,ALFAG iA m UG*v»A m UG,GAH5 , HG »hS *65 

* .KOPT , aOd .NPHASE .DELI A l *PB 

* , b .H , R I ,«6 ,R& »RF 1 L .KiPOT , DHOLE ,03P0T ,N * 00 , Hi 

* , a.SRES ,3TK»AP,hEaD) 

C 

CALL OUT AW ( HEAD . A .ShES »STR ) 

C 

1CA5E=1CASE+1 

I F < 1 CASE ,LE .NcASES ) GU TO 10 
C 

CALL L ND JOB 

STOP 

END 


WHDG 002 l <2 )(HETalS> 


ApftTiC 0021 

FoRPUR 2HH 1 -03/ 1 U- 1 ** : SB 
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0021 (2» (METALS) 

939600* TPF4 . 0u2 1 

1 SObROU f 1 N£ M£TALS(N,E.ANU,RHo.ALFA,FTY,FTU) 

2 C 

3 C TABLE UF METALLIC MATERIALS DESIGN PROPERTIES FuK 

9 C TUBES, FLANGES AND oOlTS 

5 C K.K.LtiMBACH, 2 h NOVEMBER 1972 


6 

c 



7 


COMMUU/PROMTX/P ( lo 

b 

c 



9 


OAT A l ( p ( I , j ) , j» j , 6 

10 


« 9.9 t +6. 

6*3 3 

1 1 


• 9.9 L + 6 . 

u • 3 3 

) 2 


* 9.9 E + 6. 

u ■ 3 3 

t 3 


* 9.9 E+6. 

0*33 

19 


• 2 B » j E+6 . 

0 • 30 

1 5 


* 2b » y t + 6 . 

u • 3o 

1 6 


+ 2b. w t+6. 

u*3o 

17 


+ 2 8 » y E+6* 

0 • 3 ti 

IB 


• 28»o E+6 . 

u • 30 

19 


♦ 28.y E+6. 

u * 30 

20 

c 



21 

c 

legend- al 

6 <j 6 i • T 6 

22 

c 

al 

2u 2 ^ T 3 

23 

c 

3^7 

SS 

29 

c 

A 2 6 6 

2b 

c 



26 


E * P C N , 1 1 


27 


A N U = P j N , 2 ) 


26 


RH0*P<N,3> 


2 V 


ALF A*P IN ,Hi 


30 


FTY=Pirt,S J 


3 1 


FTU=P(M,61 


32 

c 



33 


K E T U k N 


39 


END 


fcHDO 0022 


1 3 J (OASKET) 



WP«T#C 0(j22 

FURPUK 29h1-O3/I0-19:53 


o ) 


1=1 .lo)/ 







W*09B 

t 

12»S 

E-6 

• 3 5 » y 

E * 3 • 

i2 • U 

E + 3 . 

y . 0 9 8 

1 

12*5 

E-6 

* 32 *2 

£♦3. 

id • 1 

E + 3. 

u * y 9 8 

» 

12. 5 

E-6 

• so • 0 

E + 3» 

62*0 

E + 3. 

u « ti 9 8 

1 

12*5 

E-6 

. 9 7 . y 

£♦3. 

69*0 

£+3. 

u » 2 8 b 

> 

9.5 

E-6 

. 35.0 

E + 3 i 

*0 • U 

E + 3 . 

t * 286 

• 

9.5 

•41 

» 

Id 

. 30 « u 

t + 3 » 

7b. S 

E+3. 

y . 268 

1 

9.5 

E-6 

i 25 »y 

E + 3. 

ob • u 

E + 3, 

U*28g 

1 

9.5 

E-6 

» 1 3 1 *0 

£ + 3. 

2 u u • 0 

E + 3, 

y 268 

1 

9.5 

E-6 

. !2B.y 

E + 3. 

i 9b • Ufc +3 

w * 2 8 8 

I 

9.5 

E-6 

• 120*0 

E + 3 . 

1 do • 0 

t + i/ 

< 1 ) AT 

k r » 

i 2 ) 

AT 200E 




<31 AT 

KT i 

t 9 ) 

AT 2u0r 




1 5 > A 1 

kT 

■ < 6 ) 

AT 

2 OOF . 

(7) AT bOoF 


( 8 ) AT 

kT 

. 19) 

AT 

20uF, 

(lo)AT 6U0F 
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0022 


( 3 ) (GASKET) 


‘♦ 3 H 6 C 0 *TPFS,Du 22 

l 


2 

3 

4 

5 

6 
7 


23 

2 H 
2b 
26 

27 

28 
29 

3 L) 
3 1 

32 

33 

34 
3 b 


C 

C 

C 

C 


SUBROUTINE Gasket IN, E.AKG.SG.AUF a , A HU » 6 AMU, GAMS ,Hg , 5 P ) 

TABLE OF GASKET MATERIALS DESIGN PROPERTIES 
K • « • LE I MB ACH , 2 « NOVEMBER 1972 

C 0 MMON/PROG 5 K/P ( 5,9 J 


8 


DATA( (P( i ,j) , j S| 

• V ) • 

1 * 1 . 5 )/ 



9 


♦ • u 

L ♦ J . 

l O • 0 

E+ 3 . 

6,5 £ ♦ 3 » 

1 • 3 E • 3 • 

0*5 

10 


9 y • b 

« 

G.b 

i 

0 * 03125 * 

3 • 5 0 • 


1 1 


* M H t y 

E ♦ 3 • 

Ip.tf 

E ♦ 3 * 

3,7 E ♦ 3 * 

1.3 E- 3 . 

0*5 

12 


• u* b 

• 

U.b 

* 

0 • 0 6 2 5 * 

2 . 7 b • 


1 3 


* HHth 

t ♦ 3 t 

iiiiu 

£ + 3 » 

1,6 E ♦ 3 * 

1 . 3 E- 3 • 

0*5 

1 *1 


* 0 « 5 

t 

U • b 

t 

0*125 * 

2 • yii > 


1 5 


* 1 So • u 

E + 3 » 

8 • w 

E >3 • 

4 *o E *3 * 

3 • 8 t >£-5 # 

u *12 

1 6 


♦ 1 * u 

« 

1 . L» 

• 

0 *0625 * 

3 • y v i 


1 7 


• 2 a. o 

fc >6 . 

4 Q * 0 

E ♦ 3 • 

1 8 • 9 E ♦ 3 » 

9 . S vE "6 i 

0*30 

1 8 


* 1 • u 

• 

l • c 

• 

0 • w 2 b * 

b .5 U / 


1 9 

c 








20 

c 

LEbt-NU 

- ASBESTOb 

( 1 > 

1 / 32 , iy) 

1 / 16 * ( 3 ) 

1/8 

21 

c 


( <4 J 

KEL- 

FBI 

1/16 



22 

c 


( 5 ) 

CREi 

32 1 

-A 1/16 




E = P (N,l) 
AKG*P ( N , 2 ) 

S G * P ( N » 3 ) 
ALPAsPiN,**) 
A M U = P ( IM , 5 ) 
GAMUap(N,6) 
OAflSsP ( N , 7 ) 

hGaP < N , 8 ) 

SP*P ( iv ,9 ) 

RE TURN 
E NO 


(SHDG D 001 


( 1 ) ( OES 1 GN ) 


®PRT,C 0001 

FURPOR 24 hi - 03 / l 0 - 1 4 ; 53 



OGG 1 


( M ) f DtS I GIN ) 
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1 

2 

3 

H 

5 

6 
7 
b 
9 

I w 

I I 
J 2 
) J 


* TPF S> . Ouy 1 


2 1 
22 

2 3 
24 
2b 
26 
27 
2b 
29 

3 Li 

3 1 

32 

33 

34 
3b 
36 
3 7 

3 b 
39 
Hii 

4 1 

42 

43 

44 
4b 

46 

47 

48 
4V 
b w 

5 1 

52 

53 

54 

55 


C 

C 

C 

C 


SUBROUTINE DES I UN ( P , Di .T»0ELT.Pl> ,bF»FS|6F 

,ET,aNUT.RHoT»ALFaT,FTYT,FTUT 
.EF.ANUF.RHoFiAlFaF.FTYF.FTUF 
, £ tJ i ANUB . RHO B I A lF A b iFTYB . F T OB 

.Eu,aKG,Sg.aLFAG,AMUG,GAMU.GAM 3,HG,SP|HS.B3 
.KuPT »AoB ,WE IGh f »PB 

.B.H.KJ.RG.RB.kF IL .RSPOT .DHOLE .DSPOT ,NB,bG *hT > 

OtblGN ROUTINE FOR LOw PROFILE FLANGES 
Ktk.LEiNtiACh, 2« NOVEMBER 1972 

DIMENSION K 0 P T ( 10) 


i h 

C TUBE THICKNESS 




1 B 

JBOL T *KOPT ( 4 1 




1 6 

Ri=DI/2. 




i 7 

IF ( K 0 P T ( l > .EQ.O) 

GO 

TO 

4u 

lb 

I F < (COPT ( ) ) . t w . 1 ) 

bO 

TO 

I 0 

19 

I F ( N 0 P T ( 1 ) .£0.2) 

GO 

TO 

20 

2b 

I F ( NOPT ( 1 ) ,E0 . 3 ) 

GO 

TO 

30 


c 

c 


10 T*FS*P*r J /FT YT 
GO TO 40 
20 ALAMB*.7S 

T1=UF«P*K1/<ETUT*aLaMB)+2.*DeLT 
T2®PF«PnR I / (FT Y T*ALAMB) +2, *0EL1 
T’AMAAl C T 1 t T 2 ) 

GO TO 4C 

30 Tl=l*l»PF*p*Rl/(FTYl-.4*pF*P) 

T2»l.l*BF*P*Rl/(FTUT-.4*8F*P| 

T*aMAX 1 t n i T2 ) 

40 CONTINUE 

BOLT S12E 
OB* T 
10*0 

45 CALL bULT<DB,ETAU,EJAl*£TA2iA0B.i < ,H0LE»DSP0T , RSPuT • I Si ZE , UbOci , I y | 
I F ( NOP T (2) « EO ♦ 1) GO TO bo 
I F ( NOPT ( 2 ) ,£w . 2 I GO TO ag 
by E I *E T A 1 *06 
E2 = E T A 2 * OB 
IF( T .Gt*y.20) 

IF(T*LT*0*2l>) 

I F ( T . L I . « . I 5 1 
IF (T.LT.o. 10) 

IF (T«LT«y.yb) 

GO TO 7o 
60 Ei*ETA2»0d 
E 2 *E 1 

IF ( T .GE»0»2o) 

I F ( T » L T » y. 2o ) 

I F ( T . L F . y . I 5 ) 

IF ( T • l 1 . u • i C 1 
1 F ( T .y T .y *y5 ) 

70 CONTINUE 


RFIL*j.37s 
R F J L * O • 3 1 2 5 
RF I L*u .250 
RF IL“0. 167 5 
RF I L*y. I 25 


RF IL*o . 125 
RF lL=u. IwO 
RF I l * u . y7 5 
RF IL*u. 05& 
RF I L *u * u25 
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5b 

57 

bd 

by 

6u 

6 1 
62 
63 
6H 
6b 

6b 

6/ 

6 8 
69 
7U 

71 

72 

73 
79 
7b 
76 
7 7 
7 B 

7 9 
BU 

a i 
82 

8 3 
89 
sb 

8b 

87 

88 

89 

90 

9 l 

92 

93 
99 
9b 

96 

97 

98 

99 
100 
1 0 1 
102 
1 03 
i 0** 
1 Ob 
10b 
107 
1 08 
1 09 
1 10 
i 1 i 


C BOLT CIRCLE RADIUS 

C 1 *0* 002b 
C7 = o * ub 

lF(KGPr<2J»Eti.n K0«RI*T+C1*E» 

iFCKOPT ( 2 1 • t W • 2 ) RB«KI+T+2.*«FIL*E1 
C 

c FLANGE ftloTH 

ti*Nb*E2-R i 
C 

C G A S K E T iiv l l) T N 

IFCKOPT (3) • G T * u ) GO TO do 
RG ac #b*iKB-#b*D8 + Rl) 

IF (KtiPT(b) ,£U.0> BG»PF«P*«<j/(2»*t&AHU*AK6^GAM5*»C»*6F ) > 

1F(N0PT(5)«EQ.1) BG»PF*P«RG/(2t*lGAMUaAK6-GAMS*BP*P*PF*oFJ) 

I F { KOP T (&)«£«. 2 » BG = «B“*5*0H0Lt-Hl-2*®C2 

IF ( *» OP T 16) • E U . 3 1 GO TO 7 6 

KG=KB-»5*0H0Lt-*5*BG-C2 

ril®KG-,b*BG-C2 

hi 2 » K l 

IF(hii»(jT.H2l GO TO 7b 
RG*KH-.b*BG + C2 
74 «B=KG+,5*bG+.b»0H0LL+C2 
7b B»RB+E2-RI 
GO TO 90 

76 R G ® K I + • b * B G + C 2 
R 1 aKG+ • b • BG + C 2 
K2»RB- »b*0HUL£-L2 
1FIK1.lT.W2I GO TO ?S 
GO TO 74 
80 BGsT 

RGsK]-h,b*BG 

K 1 “K I +BG+2 . *BS 

R2oNB-.b*0H0L£ 

IF 1 R l . G T . K2 J «B»KG+.5®BG42 **bS+*5*DH0LE 
B*KB+E2-K i 
90 CONTINUE 

c 

c BOUT FORCES 

Pl»i. I 4 1 b 9 3 

9 b lF(K0Pri3>,GT.01 GO TO 100 
IF (KOP T ( S I .EG. 2 1 GO TO 98 

P B l = 2 « * P 1 * R G * B G • G A M U ♦ A K G 

lFl<OPr(S>.EQ«u) pB2®2»*PI* h G*0G*GaMS*Sg*gF-*'P1*KG* <, 2*PF*P 
I F ( iCOP T ( b I • tO • 1 ) PB2«2»*PI*Ro*BG»GA«S*SP*P*PF*GF+PI*KG # *2*PF*P 
PBsAMAX 1 < PB 1 , PB2 ) 

GO TO 1 l w 

98 PS2 = 2.*PI*RG*hG*GaMb*SP*P*PF*GF+PI*RG«'*2*PF*P 
PB 1 *P B 2 

SG2=PB2/ 12.*PI*bG*GAMU*KG> 

IF1SG2.LT. SGI Po l a2. *P WRG*B&*GAHU*SG 
PB® AM A A 1 ( Pb 1 , PB2 ) 

GO TO llu 

100 PB*P I •KG*«2*PF *P 
110 CONTINUE 

c 

c number of bolts 
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DOG 1 

1 1 2 
1 J 3 
i 1H 
1 lb 
1 1 6 

I 17 

I I a 

l 19 
1 20 
121 
1 22 
123 
12H 
125 
1 2 6 
127 
1 2a 

1 29 
) 30 
1 31 
1 32 
133 
1 3 H 
135 
1 36 

137 

138 
1 39 
1 *10 
1 H 1 
1*42 
1 H3 
1 HH 
i *45 
1 *<6 
1H7 
1 Hb 
l H 9 
1 50 
1 5 1 

152 

153 
1 5*4 
155 
1 56 
157 
1 5a 
i 59 
1 6o 
161 
162 
163 
1 6 *4 


( *4 1 (DESIGN) 

EN01*P0/(FTYB*AOB) 

ENB2=1BF/PF)*PB/(FTU6*AoB) 

ENB*AMAX 1 (ENB1 j E NB 2 1 

nb*enb 

IF1NB.6E.6) GO TO 115 
NBa6 

PB*6. *FT YB*AQB 
115 CONTINUE 
C 

c BOLT spacing 

enb=nb 

S=2 • *P I •KB/ENB 
1D=1 

SO VO=S/DB 

IF ( S 0 V 0 . G T • B ♦ G 1 GO TO 120 
1 F 1 SOVU.LT .ETA J 1 GO TO 130 
GO TO 1 HO 
120 IS! ZE»ISUE-1 

I F ( ISI2E.LT. I I GO To 1 Ho 
GO TO Hb 

1 30 ISU£=ISIZE+1 

1F1ISI2E.GT.1H) 60 TO 1 HQ 
GO TO Hb 
1 HO CONTINUE 
C 

c 

c flange height 

E»RB-RG 
R 0“ R 1 * • 5 • T 
TN»T/2. 

EMFU=FS*P0*E/<2.*PI*RG> 
BBaR=B-OHOLE*SQwT I DhOLE/S ) 

2ETA l=.,»0U 
2 E T A2»u • 1 a 

CAPAop TYF*BBAN/(6.*KU> 
CAPB*FTYF*ZLTA2*tT-lN)/2. 

CAPC«F f YF*ZETA 1 • 1 T**2-TN**2 ) /H.-EMFU 

RTSWBCAPb**2-H.*CAPA«CAPC 

KTsSUK T 1RTSG) 

H»(RT-CAPB)/(2.*CAPa) 

c 

c CHECK FLANGE HEIGHT 

SOVH*S/H 

1F1S0VH.GT.31 H b S / 3 . 

C 

C WEIGHT computation 

RrtT«.5*<2.*«I-*>T + ai 
a«t=ib-t)*h 

V0L=2.*PI*RwT*AwT 

weight*rhof *vol 

c 

return 

ENG 
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0001 ( 4 J (Dt.Sl(iN) 

!H0(i UQ10 <b) ( dOLfJ 


PKT.C D0I0 

UKPUK 2**h 1 -03/ JO- i 4 J S 3 
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U 0 iu <b> ibuii) 


3H6uU«TPFS 

1 

1 

3 

H 

b 

6 

7 

d 

9 

iU 
1 1 
1 2 
13 
19 
1 b 
I 6 
1 7 

1 3 

19 

20 
21 
22 
23 
29 
2b 
26 
27 

2 B 


PC 1 0 

S U a rt 0 U I INt BOLT(0 ,ETAC*£Ta) » £ T A 2 . APB , OHOL£ » USPO T , K Sr* 0 T , 1 S i L c. 
• . JdOu T « 101 


C 

C 


c 


c 


COilMON/BOL rPT/0X(l9,51 

IF ( J60LT.EQ. 1 ) CALL BTAblI 
I F ( J BOL T « EG) • 2 ) CALL BTAdL2 
IF ( JBQLT.ep. 3) CALL BTABL3 
iFlJdOLT.LE.P.Orf.jBuLT , GE , 9 1 STOP 

1 F ( ID.GT.O) GO TO Id 
1»1 

15 D I *PX (1,1) 

QD*P-DI 

I F ( PD • l£ » 0 ) GO TO 2o 
I a 1 ♦ 1 
GO TO lb 
id i = ism 

20 D*OA (1,1) 

ETAo*OX (1,2) 

ETA l »DX 11,3) 

E T A2»0 A (1,9) 

A03 = DX 11,5) 

OHOl£=U+.O 0S 
05P0T = 2 » • E T A) *0 
rtSP0T=.C62 
IF(ID.EO.O) ISl2E=I 


29 

30 

31 


HE rUH(M 
ENO 


)HOG P011 


(o) (dTABU) 


)pKT , C Poll 
JHPUK 29H1-C3/ 1C-19J53 
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4346Qu* TPFS. 00 1 l 


1 SUBROUTINE BTA0L1 

2 C 

3 C BOlT TABLE FOR OPEN DRENCHING 

4 C 14 SUES FROM ,25 In TO |.S jN NOMINAL DIaMETEK 

5 C 

6 COfiMON/BOtTO l /OH ( I 4 .5 1 

i cOmmon/boltot/o ( 14 .5 i 

8 C 

9 DAT A ( ( OX { I , J ) » J* I «S 1 , l ■ 1 , 1 4 1 / 


J 0 

• » 25 gu» 

3 • 00 

1 1 

* * 3125 . 

2.60 

12 

* • 3 7 bo i 

2*67 

13 

* • 4375 > 

2.57 

1 4 

* • 5030 » 

2.50 

15 

• *5 625 . 

2.4 b 

1 6 

♦ * 6250 * 

2 . 4 q 

1 ? 

* • 75 jQ, 

2.33 

18 

* • 87 bo . 

2-35 

19 

* 1 * 0000 * 

2.25 

20 

• 1 • 1 2 iv) i 

2.22 

21 

* 1 • 25 gq • 

2.25 

22 

* i • 3 7 So * 

2.23 

23 

* 1 • 5000 « 

2.17 

2 4 

c 


25 

00 10 1 = 1,14 

26 

00 1 u J* i , 5 


2 7 

10 0 ( 1 . J 1 = 0 X < I , 

J 1 

28 

RETURN 


29 

tUD 




2*00 

* 

1.50 . 

.03182. 


1*80 

• 

1 * 40 » 

• 0524 J » 


1.67 

♦ 

1*33 . 

0*1)774? • 


1*57 

• 

1.29 . 

• 1 0 * 3 i » 


1 • 62 

• 

1.24 , 

. 1 4 1 9 0 * 


1*56 

• 

1*22 » 

. I8l9i, 


I *50 

• 

1*20 . 

.22 600. 


i .49 

* 

1*0“ » 

. 3344q * 


1.43 

« 

1 » 0 7 . 

.46173* 


1.37 

• 

i .o 6 » 

• og 5 7 4 » 


1*33 

• 

1*00 » 

.763 27 . 


1.40 

• 

l*0o » 

• 929o» . 


1.36 

• 

1*00 * 

1 * 154flo« 


1.33 

8 

i » 

1 * 4yS2a/ 


iHOO 0012 </) (BTABL21 


SPRT.C 0012 
' URP UN 24h 1 -03/ 1 0- 1 4 • S3 
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0012 (7) ( d T A dL 2 ) 


939600MPFS..DU 1 2 

1 SUBROUTINE B T A Bl 2 

2 C 

j c bolt table for socket wrenching 

h c 1 9 sizes from .25 in to i.b in nominal oumeteh 

s c 

6 c0MM0N/t»0l_T02/0X l 1 9 » 5 > 

7 COMMOn/BOL TDT/O ( 19 , 5 ) 

a c 


9 

0 ATAIIDxIl 

. 0 1 . J * 

10 

<» »2S00* 

2 , 7 6 

1 1 

* • 3 1 2 5 . 

*1 

in 

• 

12 

• • 3 7 50 i 

2-37 

13 

• *9375. 

2.26 

19 

• * S OuO * 

2.18 

is 

• *5628. 

2. 2o 

1* 

♦ • 6 2 8 o • 

2.22 

1 7 

t • 7 SOO * 

2* 1 2 

IB 

* .8750. 

2*28 

1? 

• 1*0000. 

2.19 

20 

• 1 * 1 250 • 

2.1** 

2 l 

• 1 • 2 S do » 

2-09 

22 

• 1 * 3 7 Sq ■ 

2*00 

23 

• 1 * 5o00 * 

2 « u2 

29 

C 


25 

DO 1 0 1=1. 

19 

Zb 

DO 10 3=1. 

5 

27 

10 0 ( I »0 1 *0X » I . J 1 

2 a 

RETURN 


29 

END 



UNO G 0013 (b) (BTABL3I 


•51.1*1.191/ 



1.60 

• 

1*90 » 

.03102 . 


1 .50 

• 

1.20 » 

.05293 . 


1 *33 

• 

1 * 2 d » 

•07799 . 


l • 2 5 

« 

1*19 » 

. 1 U 6 J 1 . 


1 *20 

• 

1*10 • 

. 19190 * 


1*22 


1 . 1 1 • 

. 16199 . 


1.25 

• 

1 • 1 2 . 

. 22600 * 


1 * 1 7 

• 

1 * 0 7 » 

. 33996 . 


1.31 

• 

1 • 1 9 . 

• 96173 . 


1.25 

f 

1*10 • 

. 605 / 9 . 


1*22 

• 

l * 0 7 • 

•76327 . 


1 • 1 ti 

1 

l * 09 • 

• V 29 u 5 . 


1.16 

• 

1 » u 2 • 

1 . 1 59 o 8 » 


1 * 1 3 

• 

1*00 * 

1 . 90525 / 


SPRT.C 0013 

FURPUR 2HHl-03/lO-ll:Si 
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*<3H60a*TPFi .Do I 3 
J 

2 C 

3 C 

■* C 

s c 

6 
7 

0 C 

9 


10 
1 1 
12 
l 3 
l 4 
15 
1 6 
1 7 
10 

19 

20 
21 
22 

23 

24 C 

25 

26 

27 10 

20 

29 


SUBROUTINE BTABL3 

BOLT TABLE FOR INTEkNAL W«ENcHIN& 

n sizes from *25 in to i.s in nominal diameter 

COMMON/B0LTD3/DX 1 l 4 » 5 ) 

COMMON /B0LT0T/DU4.&I 



UA TA I < OX ( I 

. J ) » Js | 

.5) 

.1*1.14)/ 





• 2 5 >j 3 * 

1 .92 

t 

1.16 

• 

0.96 

9 

* 0 3 1 0 2 i 


*3125. 

1*06 

t 

1 .0’ 

• 

u • 9 3 

• 

» 0 5 2 4 3 • 


• 3750. 

1.79 

1 

1 * u ^ 


u*9i 

• 

•07749. 


• 4375 . 

1.00 

1 

Uu3 

• 

u • 9 1 

• 

• 10631 » 


• SdOO • 

1.70 

t 

1 • go 

• 

u • 9 u 

9 

• 141 V0» 


•5625. 

1.76 

• 

a« vd 

| 

u* 89 

t 

• 16194. 


• 6250 • 

1.75 

• 

o*96 

• 

u * 88 

t 

• 22600 • 


•7500. 

1 .60 

♦ 

0*91 

• 

u*8H 

• 

• 33446 i 


* 0 7 5 0 . 

1.69 

1 

ii * 9 0 

• 

u * 8S 

• 

•46173. 


1 • Oil JO . 

1.67 

9 

0*89 

• 

u • 8 *♦ 

• 

•60574. 


1 ♦ l 25 j • 

1*86 

1 

u* 96 

• 

g • 9 2 

t 

.763*7 . 


l » 2500 . 

1.67 

t 

0 • 8 7 

• 

u * 83 

1 

• 9 29oS • 


1 *3750. 

1 • 00 

• 

*2*93 

* 

g • 8 9 

» 

1 .15400 . 


1 • 5 0 u 0 * 

1 . 65 

• 

0 • 8 5 

• 

0*82 

• 

1.40525/ 


DO l J 1*1.14 
DO 1 j J»1 ,S 
0 ( I ,U)=OX I I ,J) 
RETURN 
END 


iHOG P001 1 9 ) (HLOTF 1 ) 


iPRT.C P001 
URpUR 24H1-03/10-14 J53 
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>001 


(9) IPL0TF1) 


19600* TPFS .PJG1 

1 

2 i 

3 - 
9 

5 C 

6 c* * * * 

7 C**«* 

a c**»* 

9 C 

1 0 
1 1 
12 
i a 

19 
I b 

1 6 C 

1 7 
Id 

19 

20 

2 1 
22 

2d C 

29 
2b 
26 
27 
2B 
29 
do 

31 

32 

33 
3 9 
3b 

36 2 10 

37 
3tt 
39 
90 

M| 2 1b 

92 

93 
99 
9b 

96 

97 
9B 
99 
bu 
bl 
b2 
bd 
59 
bb 


SUBROUTINE PUOTfl(B»H.T.Rl ,« 6 »HBt^PlL.RSP 0 T, 0 H 0 LE» 0 bP 0 T.N.Bo,HG 

► .BS.HS.HT 

> ,FTYF,FTUFiFTYb,FTUB,SG»FYG,P.«ElGHT 

, .FLAMTL.BQLmTL.GASMTL.HEAD.KOPI ) 

PLOT KOU TINE FOR L0»* P r 0F1lE FLANGE WITH FLAT GASKET AND 

machined spotfaces for The bolts 

K«R«LE IMBACH, 7 NOVlMBER 1972 

DIMENSION FLAMTl ( 21 .B 0 LMTL( 2 ) »GASMTLt 2 ) ,HEAO{ 12* 

. OA T A l I i 1 » OAT A2 i l 1 i DAT A3 < 3 l »0ATA4 ( 3 > . 0ATA5 l 3 1 

t 0ATA6(4) ,0ATA7< 1 1 1 .DAI AB( U ) .0ATA9U 11 tOATAlOl U > 

.DAT A 1 1 i 5 I .DAT A 1 2 lb 1 

» X < 100* » Y 1 1 DO ^ • I X < loOl « l T t 1 00 ) 

, KOP T ( 1 0 I 

call frameviui 

CALL SCRECT l 3 1 . 3 l .99 I ,99 | ) 

CALL PK1NTV (72, HEAD. 91 . 1003) 

A * 3 • 7 b 

CALL XSCALV ( -A , A ,0 .u I 
CALL YSCALV ( -A * A ,0 . ^ ) 

I 1 6»0 
D=DMOLE 
A 1 l >“8/2. 

Y « 1 I **H/ 2 . 

X<2)»Xll) 

Y < 2>»-Y ( 1 ) 

[ F < KQP T ( 2 1 »EQ . 2 ) GO TO 210 
X ( 3 1 *- X ( 1 ) +T + RF 1 L 

Y < d 1 *Y < 2 ) 

Xl9)=X(d>-RFlL 
Y 1 9 ) »Y t 3 ) ♦RF IL 
GO TO 2 l b 
E1«R1*B-KB 

X < 3 ) *X ( 2 ) -2 . *E 1 

Y ( 3 1 * Y ( 2 ) 

X ( 9 > =X (d > -2. *RF 1 L 
Y ( 9 1 «Y ( 2 ) 

continue 

X ( S ) s A l 9 ) 

Y < b ) * Y ( 2 ) +H T 
X (6 )»X lb >-T 

Y 1 6 > * Y 1 b 1 
X ( 7 1 s X < 6 1 
Y ( 7 >»Y l 1 > 

X(d>*Xl7>*RB-Rl+D/2. 

Y (d >»Y C 1 1 
X i 9 ) »X (B ) 

Y < 91 «Y < 2) 

X 1 lo >=A 1 B >-D 
Y ( lo ) «Y < 1) 

X( 1 1 )*X( 10) 
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POOI ( 9 1 (PlOTFJJ 


S6 


XI 1 2 ) *X < 8 ) - 0/2 * 


57 


Y (121 »Y ( 1 1-.2S 


5a 


X ( 1 3 ) B X ( 12) 


59 


Y< l 3 1 a Y ( 2 ) ♦ • 25 


66 


Y ( 1 4 ) * Y ( 2 1 


6 1 


X ( 1 4 ) » X ( 131-0SP0T/2. 

♦RSPoT 

62 


X ( 15) ax ( 14 ) -K SP 0 T 


63 


Y < lb ) ay < 2 ) +RSP0T 


64 


X ( 16 ) aX ( 15 > 


65 


IF(X<l6)*i_T*X(3l) Gu 

TO 95 

66 


X < 1 4 ) =X < 2 1 


67 


X < 1 41 ax ( 3 1 


o a 


A ( 1 5 ) a X ( 3 1 


69 


X < 1 6 ) *X ( 3 ) 


7 u 


Y ( 1 4 ) » Y ( 3 ) 


71 


Y ( 1 5 ) = Y ( 3 ) 


71 


Y ( 16)ay (3) 


7 3 


1 1 6= 1 


74 


GO TO 102 


75 

95 

COnTI hue 


76 

C 



77 


0X»M JJ-XI 16) 


76 


1 F l OX .GE.NFIL ) GO TO 

U1 

79 


QYaSQrlT 4RF1L*«2~DX*«2) 

ao 


OYBAKsKF I l-QY 


a i 


Y( 1 6 ) a y ( 2 ) + D Y 8 Aw 


82 


GO TO 102 


83 

10k 

X( 16 ) aX ( 4 1 


84 


Y ( 1 6 ) a Y (41 


85 

102 

CONTINUE 


86 

c 



07 


X(l7)aX(7)+RG-Rl+8G/2. 

88 


Y ( 1 7 1 * V < l i 


89 


X ( 1 8 1 «x ( I 7 ) -BG 


90 


Y ( l 8 ) a y ( j ) 


9 J 


X ( 1 9 ) *x 1 1 7 ) 


92 


Y < 19)ay< i j-hg 


93 


X < 20 1 »X { 1 8 1 


94 


Y(2U)=Y< 19) 


95 

c 



96 


X(21 )®x( 18) 


97 


Y ( 2 l ) = Y < 20 ) -. 1 25 


98 


X ( 22 ) »X ( 18) 


99 


Y ( 22 > = Y l 2 1 ) -.375 


100 


X ( 23 ) ax ( 1 8 ) 


lul 


Y ( 23 ) ay ( 22 ) + . 1 25 


102 


X(24|ax(6)-2.0 


103 


y < 24 ) a y ( 23 ) 


104 


X < 2 5 ) a x ( 17) 


10b 


y(25)ay (21 ) 


1 06 


X ( 2 6 ) a A 1 17) 


107 


Y(26)sy ( 2 b 1 - . 7 50 


108 


X< 27 ) ax ( 17 ) 


109 


Y ( 27 1 a y ( 26 )«•. 1 25 


1 16 


X ( 28 ) ax ( 24 ) 


1 1 1 


Y C 2 8 ) a Y ( 2 7 ) 
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A ( 2 9 ) a 

Y ( 29 I = 
X < 3ui ■ 

Y < 30) » 
X ( 3 l I * 
Y ( 31 ) * 
X<32>= 

Y < 32)* 
X < 33 ) = 
Y(33)* 
X C 3** > » 

Y ( 3 9 ) a 
x t ; b ) = 

Y ( 35 ) s 
X ( 36 ) « 
Y(36)a 
X< 3? ) « 

Y ( 37 ) s 
X<3«) ! 

Y < 38 ) « 
X < 39 ) » 

Y ( 39 ): 
X < 9U ) > 

Y < 90 ) 1 
X(9l ): 

Y (9 l ) : 
XC92) 

Y C 92 > ' 
X C H 3 ) = 
Y< 93) 
X<99); 
Y(99) 
X(**S) 
Y(95> 
X ( 9 6) 
Y< 96) 
X ( 97 ) 

Y (97 ) 
X ( 98) 

Y < 9b ) 
X ( 99 ) 

Y < 99 ) 
X ( SO ) 

Y ( bO ) 
X • S i ) 
Y(S1 ) 
X<52) 

Y ( 52 ) 
X ( 5 3 ) 

Y ( S3 ) 
X ( 59 ) 

Y < 59 I 


3 A(l )♦• 1 25 
3 Y < l ) 

“X ( 29 ) ♦•750 

* Y ( 1 ) 

*X ( 29 ) ♦. 375 
«Y ( 1 ) 
s X ( 3 1 ) 

= Y ( 1 > ♦H/2-. I 25 
3 X ( 3 ) ) 

3 Y ( 32 )+. 2S 
*X131) . 

= Y ( 2 ) 

= X ( 3 1 ) ♦ . 1 2b 
= Y < 2 ) 

= X ( 2 I ♦• 1 25 
= Y (2) 

*X( 30)-. 1 2S 
3 Y ( 1 ) 

3 X < 37 ) 

* Y ( 1 ) + (H + HT)/2. 
*X ( 37 > 
=Y(38)+.25 

“ X i 3 7 ) 

*Y(SI 
= X < 30 > 

3 Y ( b > 

=X ( 5 ) + . 1 25 
*Y 15) 

*X(3l > + • 12b 
*Y ( 19) 

3 A ( 19 ) +* l 25 
*Y < 19 ) 

M ( 3 1 ) 

3 Y ( 19) 

3 X ( 3 l ) 

3 Y 1 19 ) -• 375 
»X { 29 > 

3 Y l 1 ) ♦. 125 
= X ( l ) 

= Y ( 97 1 
= X ( 29 ) 

= Y ( 2 ) -• 1 25 
=X ( 1 2 ) 


3 Y C 9 9 1 
*X ( 29 ) 

»Y ( 5 ) -• 1 25 
3 X ( 6 ) 

3 Y ( 5 l I 
*X < 5 ) 

3 Y ( 5 l ) 

*xl 5 ) + *375 
3 Y < 5 l I 


DK*>KF1L*< 1 • -SORT ( 2 • 
X(Sb> 3 X(5)+0tf 
Y ( 55 ) 3 Y ( 2 ) +0R 


••125 


1/2.) 
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POD 1 ( 9 1 ( P L Q T F 1 ) 


» 68 



1 FI K0PT( 2 1 .£<3 . 2 1 Y ( 55 ) »Y { 2 » = 0* 

) 69 



X 1 56 1 *x ( 55 1 ♦. 25 

1 7u 



Y < 56 1 «Y 1 55 ) + . 25 

1 7 1 



X < 57 1 =x 1 90 1 +. 1 

1 7 2 



Y < 5 7 J * Y 1 5 6 1 

I 73 



X C58 > »X C 12) 

179 



Y 158 1 »Y ( 2 1 

175 



X<59)bx(121+.1U 

l 7 6 



Y <SV)*Y 1 2) + . 10 

177 



X (60)=X (90>+. 1 

1 76 



Y 1 60 ) *Y C 59 ) 

1 79 



IF (KOP T ( 3 ) , E Q . 0 ) 60 TO 310 

1 8u 



X 1 6 1 )»x< 17) 

181 



Y ( 6 l 1 »Y l l 7 ) -HS 

1 82 



X C 6 2 ) M l 1 7 1 +BS 

183 



Y 1 62 > =Y ( 6 1 1 

189 



X(63)«X(62) 

1 85 



Y C 6 3 1 = Y (171 

1 86 


310 

CQNT I NUE 

187 

C 



I 86 



NP = 6 3 

1 89 

C 



l9U 



00 10 1=1 .NP 

191 



CALL X5CL V 1 (X( I 1 * I X 1 I 1 • 1 £ RN J 

192 



CALL YSCLV1 < V t I 1 » I Y ( I ) # 1 ERR 1 

193 


1C 

CONT I NUE 

199 

c 



195 



CALL PLOTLN < 1 ,2 . 1 X . 1 Y ) 

1 96 



CALL PL0TLNl2,3,lX,iY) 

197 



1 F ( KOPT ( 2 ) . EQ • 1 1 CALL PL OT QC ( 3 , 9 , 1 X i 1 Y > 

196 



IF ( X 0 P T ( 2 1 * E Q « 2 ) CALL PlOTHC < 3 ,9 , 1 X , I Y 1 

l 99 



CALL PLOTLN ( 9 ,5 . I X , i Y 1 

200 



CALL PLOTLN ( 5 ,6 , IX . I Y > 

201 



CALL PLOTLN < 6 ,7 . I A . I Y 1 

202 



IFUQPT (3) .EQ. I 1 GO TO 320 

203 



CALL PlO T L N ( 7 • 1 • 1 X • 1 Y 1 

209 



GO TO 325 

205 


320 

CALL PLOTLN l 1 8 , 20 . I A , I Y ) 

206 



CALL PL0TLN(2tl. 19, I A , I Y ) 

207 



CALL PLOTLN ( 19.I7.IA.IY) 

2 08 



CALL PLOTLN < 1 7 . 6 1 . I A , I Y ) 

209 



CALL PLOTLN l 61 , 62 , I A , I Y ) 

210 



CALL PL0TLN(62,63. IA, lYl 

21 1 



CALL PLOTLN (63.1. I A * 1 Y ) 

2 1 2 


325 

CON T 1 NUE 

213 



CALL PLOTLN ( 8 .9 . I X . 1 Y ) 

219 



CALL PLOTLN ( 10. 1 l . I A . 1 Y ) 

215 



IF (KOPfl 2 1 .EQ.2 1 GO TO 220 

2 1 6 



CALL PLOTLN! 3, 19, IX , I Y ) 

217 



IF ( 1 lo.EU. 1 1 GO TO 96 

2 t 6 



CAUL PLOTUC 119 , 1 5 . IA , l Y 1 

219 


96 

CONT I NUE 

220 



CALL PLOTLN! 15» 16.IA, IY) 

221 


220 

CONT INUE 

222 

c 



223 



CALL OASHLN ( 1 2 « 1 3 , I A . 1 Y 1 
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224 C 


225 

I F ( KOP f C 3 ) .EQ. 

1 ) GO 

TO 330 

22 b 

CALL 

PLQTLN ( 1 7 

. 1 9. IX 

• IY) 

22 7 

CALL 

PlOTLN < 1 6 

• 20* I A 

• IY ) 

226 

CALL 

PLQTLN! 19 

• 20 . I A 

• IY) 

22 9 

330 CONTINUE 



230 

CALL 

PLOT L N ( 2 ) 

• 22, I A 

• IY) 

231 

CALL 

PLOTLN 1 23 

• 24, IX 

. IY) 

232 

C ALL 

PLOTLN ( 25 

• 26.IA 

• IY) 

2 33 

CALL 

PLOTLN! 27 

• 28. I A 

• IY) 

234 

CALL 

PLOTLN ( 29 

♦ 30 . I A 

• m 

235 

CALL 

PlOTLN ( 3 1 

. 32 . I A 

, IY) 

236 

CALL 

PLOTLN ( 33 

• 34. I A 

. IY) 

237 

CALL 

PLOTLN I 35 

• 36 . I A 

. IY > 

2 3d 

CALL 

PL0TLN(37 

.38 . I A 

• IY) 

239 

CALL 

PLOTLN 1 39 

.40. I A 

. IY) 

240 

CALL 

PLOTLN t 4 1 

• 42. I A 

• IY) 

24 1 

CALL 

PLO TLN ( 43 

• 44. IA 

, 1 Y ) 

242 

CALL 

PLOTLN ( 45 

• M 6 . I A 

, I Y ) 

243 

CALL 

PLOTLN ( 47 

• 4B.IA 

• IY) 

244 

CALL 

PLO TLN ( 49 

• 50. I A 

• IY) 

245 

call 

PLOTLN ( 5 1 

.52. I A 

• IY) 

246 

CALL 

PL0TLNI53 

.54 . I A 

. I Y ) 

247 

CALL 

PLOTLN ( 55 

,56, I A 

• IY) 

246 

call 

PLOTLN ( 56 

.57 . 1 A 

• IY) 

249 

CALL 

PLO TLN ( 58 

• 59 , I A 

• IY) 

250 

CALL 

PL0TLN(59 

.60* J * 

. IY ) 

25 1 

C 




252 

CALL 

PLOT AR ( 1 . 

23 * IX . 

I Y ) 

253 

CALL 

PLOTAR ( 1 . 

2 7 • I X , 

IY) 

254 

CALL 

PLOTAR ( 1 , 

48, 1 X . 

IY) 

255 

CALL 

PLOTAR! 1 , 

50. IX , 

1Y> 

256 

CALL 

PLOTAR ! 1 . 

S 2 « I X i 

I Y ) 

257 

call 

PLOTAR ( 2 . 

53, IX . 

I Y ) 

258 

call 

PL0TART5, 

ss.ix, 

I Y) 

259 

call 

PLOTAR < 4 . 

3 l « I X , 

I Y ) 

260 

call 

PLUTaR (3. 

3H, IX, 

IY) 

26 1 

call 

PLOTAR l 4 , 

37, IX, 

IY) 

262 

call 

PLOTAR < 3 . 

40 . I X , 

I Y ) 

263 

call 

PLO T AR ( 3 . 

R5.IX , 

IY ) 


264 

265 

266 
2 67 
266 
2 69 
270 

27 1 
27 2 
27 J 
274 
2 75 

276 

277 


C 


0A r A ( OAT A I ( I) , I ■ 
DATA ( 0ATA2 < I) . 1 = 
DATA( UATA3 ( 1) ♦ 1 * 
QATA(0 aTA 4< I ) , la 
0ATA(0ATA5( I ) • I a 
0 ATA < DAT A6 ( I ) . I * 
DATA ( DaTA7 < 1 ) . I = 

* K S 1 , FfUa 
0ATA(0aTa 6( I ) . la 

* KS J » FTUa 
0ATA»UATA9( 1 ) . ! a 

* JNG STRESS" 

OATA (DATA lot I ) . I 


. I ) / 6H 01 A/ 

. I ) /6H « / 

• 3 ) / 1 6 h DlA. HOLES/ 

,3)/l8ri 0 1 A SPOTFACE / 

♦ 3 ) / 1 8h R FILLtT / 

,4)/24h PRESSURE PS1G 

,m/66H flange material 
ksi / 

.11) /66H BOlT MaTEH I AL 
KSI / 

.!i)/66H GASKET MATERIAL 
KSI / 

1 i 1 I J / o6H 


276 ALO STRENGTHa KSI / 

279 DAT AIOATA 1 1 I I ) . I » I .5 ) /3uH WE I 6H T OF FLANGE 
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POu 1 

280 

28 1 
282 

283 
289 
285 

284 

287 

288 

289 

290 
29 1 

292 

293 
299 

295 

296 

297 

2 9 8 

299 

300 

301 

302 

303 
309 
308 

306 

307 

308 

309 

3 1 0 
31 1 
312 
313 
319 
3 lb 
3 1 o 
3 1 7 
318 
3)9 
3 2 u 

321 

322 

323 
329 
32b 

326 

327 

328 

329 

330 

331 

332 

333 
339 
33b 


(91 (PlOTFI) 

DAT a l O aT a 1 2 l I } , I »1 ,S 1 /30H PRESSURE ENERGIZED SEAL 
C 

DG0=2.*RG+B6 
06 J «D60-2.*BG 


01=2 

. *K I 




08 = 2 

. *NB 




OF 0 = 

0 I +2 . *8 




HttELD»H.HT 




call 

PLO TLB C 28 

.29.5 

» I X 

» I Y » DGQ > 

CALL 

PlOTTX (28 

,56,5 

.IX 

*IY .DATA1 ,6) 

CALL 

PlOTLB < 29 

,24.5 

.IX 

» 1 Y , DG I ) 

CALL 

PLOTTX < 29 

,56.5 

• IX 

» JY.0ATA1 .6) 

CALL 

PLOTLB (97 

,29.5 

. I x 

, 1 Y , OFO ) 

CALL 

PLOTTX <97 

.56,5 

.IX 

» I Y .DAT Al ,6 ) 

CALL 

PLOTLB ( 99 

,29.5 

,IX 

, I Y ,08) 

CALL 

PLOTTX! 99 

,56,5 

. I x 

, 1 Y , 0 A T A 1 ,6) 

CALL 

PLOTLB ( b 1 

,29.5 

. I x 

, 1 Y ,0 I ) 

CALL 

PLOTTX (SI 

,56.5 

• IX 

,1 Y ,0ATA1 ,61 

CALL 

PLOTLB (53 

.29,5 

. IX 

. IY ,T) 

CALL 

PL0TLB(57 

.20,5 

.IX 

. lY.KFlL) 

CAUL 

PLOTTX ( 57 

.52.5 

. I X > I Y .DATA2 .6) 

call 

PLOTLB ( 60 

.29,5 

.IX 

, I Y .DHOLE t 

call 

PLOTTX! AO 

>56.5 

• IX 

. I Y , 0 A T A 3 , 18 ) 

1 P L = 

1X160)^120 




JPL* 

i Y l 60 ) +b 




ON = N 





IF < KOPT ( 2 > ,£Q . 

2) GO 

TO 

390 

CALL 

PL0TL8<60 

.29,- 

15, 

I X , 1 Y ,0SP0T > 

CALL 

PLOTTX (60 

,56.- 

15 , 

IX,1Y,DATA9.18) 

call 

PlOTLB ( 60 

,29,- 

•*5 » 

IX, iy.rspoti 

CALL 

PLOTTX ( 60 

,56.— 

35, 

I X , 1 Y » DA T A5 » 1 8 1 


390 CONTINUE 

CALL PL0TL8I32. -29,10.1*. IY*N) 

CALL PLOTLBl30,-29,lO.lX,IY»H«ffcLD) 

CALL PLOTLB I 93 , 29 ,-20 . I X » I Y ,*G 1 

CALL PLOTTX ( 26 ,-300,-b0. 1 X , 1 Y ,DA TA7 ,66 » 

CALL PLQTTX126, -306.-70. 1X,Iy,DATAB, 661 

IF < KQPT < 3 ) ,GT .0 I GO TO bo 

CALL PLOTTX! 26 , -300. -9 0 ,1X,Iy.DaTA9»64) 

CALL PlQTTX (26 , -300. -I 08, IX »| Y ,0 AT A 10,66) 

GO TO b i 

50 CALL PLOTTX < 26 ,-300*- 108 , I X , 1 Y ,OaT A 12 ,30 ) 

51 CONTINUE 

CALL PLOTTX C 26.-156,-50. I X , lY .FLANTL , 12 1 
CALL PlOTTX(26»-156,-7o.1X,1y .BOLMTL, 12) 

1 F ( KOP 1 ( 3) .GT .0 1 GO TO 60 

CALL PLOTTX (26,-156,-90, IX »lY .GASMTL, 1 2] 
40 CONTINUE 
C 

ftyf=ftyf/ iOUO. 

FTUFsFTUF/ 1000. 

FTYB=FTY8/1000« 

FTU 8 =FTUB/ 1000. 

SG=SG/ i 000 • 

F YG = F YG/ 1000. 
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Pool 

( 9 ) 

(PlQTFI ) 

336 


CALL PL0TLB1 2 6 .-12 ,-50» 1 X . I Y .FT YF J 

337 


CALL PL0TL&(26ill6»-5O»lX,lY,FfUF) 

338 


CALL PL0TLB<26»-l2,-7p»lX.lY,FTYB) 

33V 


CALL PL0TLB126* 116.-70* IX, JY.FTUB) 

3R0 


IF(K0PT<3>.GT.Q> 60 TO 70 

3 V 1 


CALL PLOTLB l 26 , l t 6 »-90 , l X » I Y ,S6 1 

3R2 


CALL PlOTLB ( 26 , 1 1 6 » - 1 0? » 1 X , l Y »F Y 6 > 

3V3 


70 CONTINUE 

2HH 

C 


3V& 


CALL PL0TTXl6,*2u0.5c,lA,IY.0ATA6, 

3R6 


JXP«IX<6 1-120 — ■ 

3**7 


lYP*IY(6)+50 

3Vb 


CALL LABL V ( P , I XP. 1YP.V, 1 , V I 

3RV 


CALL LABLV ( DN, I PL . JPL , 3 . 1 , 3) 

3&0 

C 


35 J 


CALL PLOTTX 16,-200.30. U.I Y. DAT All 

352 


I XW= U t 6 ) -56 

353 


I Yin* I Y ( 6 ) +30 

35V 


CALL LABLV( HEIGHT, IXW.IY6. 6, 1.21 

355 

c 


356 


RETURN 

357 


ENO 


ftHDG P002 ( 1 4> I iPLOTHCl 


URKT.C P0D2 

FURpUR 2HHl-03/10-lR:53 
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PO U2 llljl (PlOTHC) 


^3H6uO*TPF S 
1 

2 

3 

H 

b 

6 

7 

& 

9 

iU 
\ 1 
] 2 
1 J 

1 H 
i b 
i 6 
1 7 
1 b 

1 9 
2Q 

21 

22 

23 

2 H 
2b 
2 A 

27 

2b 

29 

3b 


Pb02 

SUBROUTINE PLOThC ( I a , 1 B , 1 X , I V ) 

C 

C • • * • PLOT A HALF CIRCLE FROM ja To lb 
C • • • * K * K »LE I MBACH , Q NOVEMBER 197? 

C 

0 I MENS ION 1 X ( 1 ) . 1 Y < 1 ) 

N-2U 
A N = N 

AP=3.lHl&926/AN 
1R» 1 IX ( 1 AJ-IX I lb) )/2 
II *1 X ( 1 A 1 
J i a I Y C 1 A ) 

N P 1 * N ♦ 1 

00 lb J»l ,NP1 

AU»U 

AP J* I AU- l • ) * A P 
CP J®COS < AP J » 

SP J*S 1 N ( AP J ) 

M = I R 

0 X = K • ( 1 ,-CPJ) 

DY®R*SPJ 

1 D X “0 X 

1 0 Y ®0 Y 

I2« I X ( I A ) - I OX 
J2=IY( IA1-I0Y 
CALL L I N£ V ( i 1 o J 1 , I 2 . J2 > 

I 1 ® 1 2 
J1»J2 

lb CONTINUE 
C 

return 

END 


«HD& PC03 


U 11 1 P L 0 T L N > 


®PRt,c P003 

FURPUR 2 < *Hl-b3/lC-lH;b3 
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H3HaQu*TPFS.P063 

1 

2 c 

3 C***« 

H c***» 

b c 

6 

1 

8 C 

9 

10 


S U b K 0 U I I N E PL0TtNllA.l8.IX.lv) 


PLOT A LINE 

K.K.LElMbACH, fl NOVEMBER 1972 
DIMENSION IX( 1) . ivm 

CALL LlNEVnXCiAl.miAl.IXUBI.lYClB)) 

return 

END — 


8HD6 POOH 


(12) (PLOfLB) 


»PRT.L POOH 

PURPUR 2HHl-C*3/10-lH;b3 
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POO 4 < i 2 ) l P LOT Lb ) 




1 

z 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
I 2 


C 

c • ♦ * ♦ 
c* * • * 

c 


c 


SUBROUTINE PL0TLBII.NX»NY,|*,IV*Z> 
PLOT LABEL 

K.K.LE1MBACH, 8 NOVEMBER 197? 

DIMENSION IM1I.1YII) 

I XP* l X ( II + NX 
I Y P * I Y ( I ) + NY 

CALL LABLV < l , 1 XP , I YP ,7 » I .3 ) 

RETURN 

END 


BHOG PQG5 


IUI (PLOTuC) 


BpRt ,C PQuS 

FURPUK 2 MHI-CJ/lO-n;S 3 


4-31 



LMSC-HREC TR D3 06492 


PQ05 

113) IPtoTWC) 

< *3‘«6Cu*TPFS.PuOS 


i 


SUBROUTINE PLOToC < I A . I B , 1 X , 1 Y ) 

l 

C 


3 

C*«** 

PLOT A QUARTER CIRCLE FROM U 

H 

c**»* 

k.r«leimbach» a November 197 ? 

5 

c 


6 


DIMENSION 1 X ( 1 ) . I Y ( 1 ) 

7 


N° l 0 

B 


an*n 

9 


A P * 3 « 15926/(2. *AN) 

1 u 


1 R = I X ( 1A)*IX( IB) 

1 i 


1 l = 1 X ( I A) 

12 


J 1 a 1 Y ( 1 A ) 

1 3 


NP 1 *N + 1 

1M 


00 10 0 = 1 i NP 1 

15 


A Ja J 

1 6 


AP J» ( A J- 1 . ) ♦ AP 

17 


CPJ*CQS l AP J ) 

IB 


SP J a S 1 N ( AP J ) 

1 9 


R a 1 K 

2U 


dx=k*spj 

2 1 


0 Y»R* ( i * -CP J 1 

22 


l 0 X “D X 

23 


1 D f=DY 

2*i 


I 2“ 1 X 1 1 A ) - 1 OX 

2b 


J 2 = I Y ( 1 Al + IOY 

2 6 


CALL L INEV ( 1 1 , J 1 » 1 2 • 02 > 

27 


I 1 a 1 2 

2 6 


J l *02 

29 

10 

CONTINUE 

3D 

c 


31 


RETURN 

32 


END 


IShOo P006 (In) (PLUfAR) 


®prt.c pou6 

FUKPUR 24H1-0J/ i£>-n :b3 
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(19) (PloTaK) 


939600*TPFS.P006 


1 

2 

3 

9 

5 

6 
7 
6 
9 

10 
1 1 
1 2 
1 3 
l 9 
lb 
16 
1 7 
18 

19 

20 
21 
22 
22 
29 
25 
2o 

27 

28 

29 

30 

31 

32 

33 
39 
3b 

36 

37 

38 

39 
90 
9 1 

92 

93 
99 

95 

96 

97 

98 

99 

50 

51 

52 


C 

c • ♦ • • 

C * * 

c 

c 

c 

c 

c 

c 

c 

c 


10 


20 


30 


HQ 


SO 
1 00 


c 


SUBROUTINE PLOTaH 1 lURNT , i , I X , | Y ) 

PLOT an ARROW HEAD FOR DIFFERENT ORIENTATIONS 
K.R.LElMBACH, 7 NOVcMBEr 197* 

1 ORNT POINTING 

1 RIGHT 

2 LEFT 

3 UP 

9 DOWN 

5 DOWN-LEFT 95 DEGREES 

DIMENSION IXU).lY(l) 

I 1 s 1 X ( I ) 

Jl = l Y ( I 1 

IFlIORNT.NE.il GO To lo 
12 = 11-10 
02=Jl+5 
13=12 
J3 * J l -5 
GO TO 100 

IF < 10RUT • N E » 2 I GO To 20 
12 = 11+10 
J2«Jl+5 
13=12 
J3= J 1 -S 
GO TO 100 

IF ( IQRNT.NE.3) GO To 36 

12=11+5 

J2 = J 1 - 1 c 

13=11-5 

J 3 = J 2 

GU TO 100 

I F (I ORNT . NE . 9 ) GO To 9u 

I 2= I 1 +5 

J2»J 1 + 1 0 

13=11-5 

J 3 = J 2 

GO TO 100 

1 F ( 1 ORNT .NE .5 1 GO TO 5u 

12 = 11+11 

02= J 1 +3 

I 3 = 1 1+3 

J3» J l + 1 l 

GO TO 100 

RETURN 

CALL L 1 NE V 1 1 1 . J I ,12,J2) 

CALL L1NEV ( II . Jl , 13. J3l 
CALL L INEV (12.J2.I3.J3) 


RETURN 

END 



LMSC-HREC TR D306492 

poo* iim iplotahj 

«HD6 P0O7 (lb) ( PLOTT X ) 

SPHT.C PQu7 

FURPUR 2HHl-U-3/lO-lH:i33 
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l 

i c 

3 C * • * • 

H C • * • • 

b c 

6 

7 C 

a 

9 

lu 

I 1 C 

I I 
13 


SUBROUTINE P|lOTTXCI,NX»NY,|X,IY,aK,NTEXT) 


PLOT TEXT 

K.R*LElMbACH, 8 NOVEMBER 197? 

DIMENSION I X 1 1 ) , I Y ( 1 ) • AR ( 1 ) 

IPLT»1a( I ) +NX 
jPLT* I Y ( I ) + N Y 

CALL PK INTV < NTEXT . AK« IPLI , JPlT ) 

RETURN 

END 


WHD(j P008 


(lb) (DASHLN) 


»PRT.C P008 
FURPUR 2*)Hl-03/ lo-l*) 
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POOO (16) (OASHLN) 


♦ 3460lJ*TPFS. 

Po08 


1 



S U 8 K 0 U T l N £ OaSHLN ( 1 A • IB . I X 

2 

c 



3 

c 


K • K • L • - 11/30/72 

4 

c 



6 



01 MANSION I X ( 1 ) . 1 Y ( 1 ) 

6 

c 



7 



1 1 » 1 X 1 1 A > 

a 



JI«I Y ( 1 A ) 

9 



j2 a J 1+26 

1 u 


1 0 

CALL L 1 N£ V ( I 1 ,J1 , I 1 .02> 

1 1 



Jl*J2*25 

12 



J2*J 1 *25 

J 3 



IF ( J2.LT • I Y ( I B ) ) GO TO 10 

14 



IF < J1 .LT • 1 Y( IB) ) GO TO 15 

lb 



GO TO 20 

1 6 


15 

J 2 = I r 1 IB) 

1 7 



60 TO 10 

18 


20 

continue 

19 

c 



2U 



J1=I Y< 1 A ) + 3 5 

21 



J 2 = J 1 *5 

22 


30 

CALL L 1 N E V ( Il.Jl.lt . J2) 

23 



Jl»J2+45 

24 



j2 » J 1 ♦5 

2b 



1FIJ2.LT. i Y C IB) 1 GO TO 30 

26 



1 F ( J 1 ,UT • l Y ( I B > > GO TO 35 

27 



GO TO 40 

28 


35 

J2-1 Y 1 IB ) 

29 



GO TO 30 

30 


40 

CONT l NOE 

3 l 

c 



32 



RE1 UkN 

33 



EHO 


®HD6 DluO 117) (OUTOES) 


SpKT.C t>lUU 

fUKPUK 24H1-03/10-14J53 
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DlUO 

3 4 6 (JO ♦ 
1 
2 

3 

4 

5 

6 
7 
U 
9 

1U 
1 1 
1 2 
1 3 
1 4 

15 

16 
1 7 
IB 

1 9 
20 

2 l 
22 
23 
29 

25 

26 
27 

2 B 

29 

30 

3 l 

32 

33 
39 
3 S 

36 

37 
3b 
39 
90 


U 7 > (QUTUES) 

TPFS.DlQO 

SUBROUTINE OUTDES < HEAD »A0B ,*El &HT iKOPT.T 
• .B.H.Ri , RG.RB.RFil.RSPGT. DHOLE, D>PUl .NB.bG.nT ) 

C 

C*m« OUTPUT OF DESIGN PARAMETERS 
C**«* K R L 1 / 8/73 
C 

DIMENSION HEAD ( l 2 I .KOPT ( 10 I 
1 0 I FORMAT ( I 2 A6 > 


D I = * i F 1 0 • 9 • • IN*/ 
OGI** .Flu, 4, • IN'/ 
0 G 0 a • , F 10.4, ' IN*/ 
qB» ' . F 1 u » 4 » * IN'/ 


FORMAT ( • 

AUB* ' »F 1 U . 9 • » 

SQ» I N * / 


• it EIGHT*' »F10«4,' LB 

'/ 


« 

B«* .Flu. 4* » 

IN'/ 


1 

H=« ,F !u.4 . ♦ 

IN'/ 


f 

R I * * ,F 10.4. » 

IN' , » 


t 

RG*» .Flu. 4 . ' 

IN' . ' 


• 

V 

RB« • »F lu.4 » * 

IN* .' 


• 

RF IL»» »F lu.4 . ' 

IN'/ 


1 

R S P 0 T “ ' .F10. 4 i ' 

IN*/ 


1 

DHOLE = * »F lu.4 . • 

IN'/ 


• 

DSPOTa* . F 1 u . 4 . ' 

IN'/ 


f 

NB« ' * I 1 0/ 



♦ 

b&= • »F lu. 4 i * 

IN'/ 


« 

HT» • .F 10.4 • ' 

IN'/) 


103 FORMAT 18F 1 0 • 4 ) 

01=2. *KI 
0GI=2.*RG-BG 
DGQ*2 . *RG+BG 
DB=2.*KB 

rtRlTEl6,lul> (HEAOU ) . 1*1 » l 2 l 
IF (KOP r ( 3 ) .GT .01 GO TO 50 

WRITE (6, 1 o2 ) aOb*w£»GHT , B • H , R I iUl » RG , DG I .dGO.HB.Do 

• ,RF IE ,R6P0T iDHGLE |U5pOT »NB»BG,HT 

50 WHITE 16, 1 U 3 I AGO, WEIGHT, T»B»H*K 1,01 .RB.DB.RFIL.kSPOI , DHGLE.uSPO 

• « N b * H T 
RETURN 

END 


fclHDG AUDI 


(lb) (ANAEY5) 


BPK T ♦ C AOO l 

FURPUK 29Hl-03/10^|9:S3 
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tie) (analyS) 
H3160G*TPF$. AUGl 


l 



SUBROUTINE ANALYStP.OI .T .DELT.PF ,BF,FS,GF 

2 



• ,E f , ANUT ,RHoT • ALFA! .ftyt .ftut 

i 



* ,EF , ANUF ,«HoF . ALFAF ,FTYF .FTUF 

M 



* .Ed.ANuB.RHoB.ALFAb.FTYB.FTuB 

5 



• f E8,2Kb.SG,ALFAG.AMUG,SAMU,GAfib,H(j,HS*bS 

6 



« .KuPT,A0B,NpHASE.D£LTAT,P8 

7 



* .B.rtiRJ ,RG»RB,RFlL.RSPQT,DHOLE,OSPDT,NB,B<i 

a 



♦ ,A.SRES»STR.AP.MtAt)> 

9 

c 



lu 

c * * * ♦ 

STRESS AND UEFORrtATiQN ANALYSIS 

1 1 

c 


K.R.LEIM&ACH, 5 JANUARY J973 

l 2 

c 



1 3 



DIMENSION A ( 9 , M ) * SR l S ( 5 » M ) , S T R l 5 * M ) * AP t 8 ) t HE A 0 ( 12) .KOPT ( ID) 

i M 



PI 3 3 • 1 M l 59 

15 



RO 3 Kl ♦ T/2. 

1 6 



FX»P*Ko/2« 

17 



FR=P»h*R I /R j 

l b 



RS«RG-BG/2. 

1 9 



FP»P* ( RS* + 2-R 1 ** 2 ) / t 2 • *R0 ) 

Zu 



RP 3 (RS**2+Rl*RS+Kl*+2)/< 1.5*<RS + RI ) ) 

Zl 



AG=2 . »P I •KG*aG 

22 



IF(K0PT(3) *GT.O) GO TO Stf 

23 



£KG*AG+£G/ ( 2»+PI+R0*HG) 

2H 



GO TO bl 

2b 


5C 

£KG = AG*EF/(2.*Pl*r*G*<HG + HS)) 

26 


51 

COnTI nde 

27 



ELB = H 

28 



enb = nb 

29 



AbaENB* AOB 

3u 



EKB 3 AB*Eb/(2.#Pl *RO*E lB ) 

3 1 



£»KB-RG 

32 



RAaKG+EKB»E/tE<a+EKu) 

33 



CE»EKB*EKG*£*+2/ (EKB+EKG ) 

3 H 



AFad *H 

3b 



AIF=AF*H**2/12. 

36 



C = H/2. 

37 



RC»R I +8/2 > 

38 



6EN0 = ef*r**3/<12.+ (i •-ANUT**2 > ) 

39 



AKM*i2.*(|.-ANUT**2l/(Ru»*2*T**2) 

Mu 



AK2 3 SURT < AKM ) 

M 1 



AKoSOR T ( AK2 1 

M2 



BETAe8END*AK+R0+RC/tEF*AIF) 

M3 



0X»U .+8ETA) •( .5/AK2+BET A*(C**2+AIF/AF )l-{C*eETA-.5/AK ) **2 

MM 



CFa8END*DX/(BETA*(BtTA*AlF/(AK*AF)+»25/(AK*AK2))) 

M5 



C WF * * 5 / ( A K 2 + A K * t* E N D 1 

M 6 



BET AOXaBETA/UX 

M 7 



CMF= BET AD X • ( . 5 / A* 2 + b£ T A * A I F / AF ♦ • 5 • C / AK ) 

MB 



COF= BtTADX * ( C + « 5 / AK ) 

M 9 

c 



5u 


101 

J= 1 

51 



BFC»Pd 

52 



EAM*£*aFC/(2.«Pl+RG) 

53 



AOFLau. 

5 M 



Crt I Q sEAM/CF 

55 



FRUTsCHIG 
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AUG 1 

US] 

1 ( AnALYS ) 

56 



BSTR0=BFC/AB 

5 7 



GSTR|J»BFC/AG 

SB 



BSTKSbBSTKO 

59 



GSTRS»GSTRO 

60 



ENX»u* 

61 



WP*0. 

62 



EMX=CMF*EAM 

63 



EQX=CuF*EAM 

6S 



RDFLa + CWF* (E(JX-aK*EMX 1 

6 5 



ENY= ET*T •ROFL/RO 

66 



EMY= ANUT *EMX 

67 



GO TO 201 

66 

c 



6 V 


102 

J»2 

70 



EAMlsFP* (RA-RP) +FX» IRA-Kdl 

71 



DW=P*R u »*2*(l.-ANuT/2.)/(ET*Tl +R0* ALF A T *DEL T A T -«0*RC *FR / l t F* AF ) 

7 2 



EAM2*CF *BETA* ( C*« B/AK J ♦Oift/OX 

73 



AOFLn l p P + F X ) / (E<G+Ea6) 

79 



CHIP* 1 E A M 1 ♦ E A M 2 1/ (Cfc + CF 1 

75 



FRUT»CHI0+CHIP 

76 



og=aofl+chip*ekb«e/ie<b+ekg) 

77 



OB»ADFL-CHlP*EKG*E/lEKB+E«G) 

76 



aSTRS=«BSTKO + EB*DB/EuB 

7V 



6STRS=GSTR0-EG*OG/HG 

00 



BFC«B5TRS*AB 

01 



E M X “ F X 

82 



*P«P*ftO**2»< 1 .-ANUT/2. >/ (ET*T ) +RO*ALFAT*oELTAT 

83 


1 12 

EAh=CF*FROT 

QH 



emx»cmf*eam 

8b 



EOX«CQF*EAM 

86 



RDFL» + t.WF*IE^X-AK*EHX)+wP 

87 



ENY* EI*T«RUFL/RO 

80 



EMY= ANU f*EMX 

8 V 



GO TO 201 

90 

c 



91 


103 

J = 3 

92 



D«i»t>w-K0*ALFAT*0ElTA t 

93 



E AM2«CF*BET A* ( C+ .&/AK 1 *Ow/OX 

9 4 



CHIP* ( EAM l +EAM2 1 / < Cfc+CF 1 

9b 



FROTaCH I 0 + C M I P 

96 


1 1 3 

0G«AdFL + CHIP*E<B*E/1EKB + £)<G) 

97 



Db=AOFL-CHIP*EKG*E/lEKB+EKG) 

96 



BSTRS»BSTR0+EB*0B/EuB 

99 



6STRS=GSTR0-EG»OG/HG 

1 00 



bfc*bstrs*ab 

1 0 1 



lFU.tw.3) WP«1YP-R0*ALFAT*0ELTAT 

102 



I F ( J. EG « H ) *pBpF*P*Ro**2*< l»-ANuT/2*)/(ET/T) 

1 03 



60 TO U2 

1 U 1 * 

c 



10b 


1QH 


l 06 



CH1P*PK*CHIP 

1U7 



FROT «CHI O + CH I P 

108 



aofl«pf *aofl x 
£NX«PF*FX 

109 



1 1 0 



GO TO 113 

1 1 1 

c 
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AGOI ( 16 ) (ANALYS) 


1 2 
i 3 
IH 

15 

16 
17 
1 B 
1 9 
2G 
21 
22 
23 
2H 

25 

26 
27 
2b 
29 
3 G 
3 1 

32 

33 
3H 

35 

36 

37 
3 B 
39 
HO 
HI 
H 2 
H 3 
H H 
Hb 
H 6 
H 7 
H 6 
H 9 
bu 

51 

52 

53 

5 H 

55 

56 

57 
56 

59 

60 

6 1 
62 
63 

6 H 
65 
6 1 > 

07 


201 KTOPbKDFL 
ftBQT*ROFL-H*FROT 
SFTOP«£F *«TOP/RO 
SFBQT*EF»«BO T/Ro 

IF (K0P1 < 2 > ,£Q. 1) 60 TO 202 

S X I 6 • • E MX / T**2*ENX/T 

5X0 a *6»*EMX/T**2*ENA/T 

SYI a -6.*EMY/T*»2*ENY/T 

SY0 a «-6.*£MY/T*«2*ENY/T 

TXZ = 1 .5*£0X/T 

GO TO 203 

20 2 fX*T+»5*RFlL 
SXl=-6.*EMX/TX*»2+EwX/TX 
SXO* 6.*£MX/TX**2 +EnX/TX 
SYI=-6.*£MY/TX*#2+EnY/Ta 
S YO® 6. *EMY/TX**2+E«Y/Ta 
TX2«1 •5*EQX/TA 

203 CONTINUE 

c 

A 1 1 , J ) *bFC 
A ( 2 • J )<*EAM 
A ( 3 . J ) =ADFL 
A ( H , J ) sRqFL 
A ( 5 . J >=FROT 
A<6»J)=asTRS 
A ( 7 tJ)=G5TRS 
A ( 6 t J ) =SF TOP 
A<9 »J)*SF6QT 
C 

SR£S ( 1 , J ) =ENX 
SR£5 ( 2 , J ) «ENY 
SRES ( j , j ) aEriX 
SRE.S ( H , J ) «£MY 
SRE5(S, J)aE<SX 
C 

STR < 1 , J > aSX 1 
5 Tk ( 2 , J » *SY I 
STR(3,J)«SX0 
ST« ( H , J ) YO 
5TR < 5 » J ) «T XZ 
C 

AP ( 11 =bFC 
AP ( 2 ) s£ AM 
AP ( 3 ) cAOFl 
AP ( H ) «WP 
AP ( 5 ) *>FROT 
AP ( 6 ) *lNX 
AP C 7 ) *£ M X 
AP(6)a£ux 
l PHASE* J 

CAUL P LOTF 2 ( D I , b , H , R B . RF 1 L . T , H T , ALF A T , DEL T AT 

* , I PHASE , AP ,£T , ANUT *AK ,£F 

* ,HE AO ,KOPT ,P f PP) 

J* J + 1 

IF (J.GT oNPHASE ) GO TO 3to 
I F ( 3 • £ 0 « 2 ) GO To 102 
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1 6b 


IF ( J, 

£0*3) 60 TO )03 

169 


I F ( J » 

EU.V) GO TO 10R 

i 7U 

C 



i 7 1 

c 

END OF STRESS AND DtFORMATlON ANALYSIS 

1 11 

c 



173 

c 

ultimate moment capacity 

i / H 

300 

TN = 

0 • S * T 

1 7b 


ALFA | 

a 1 T * • 2 -TN**2)/4, 

176 


ALFA l 

» ALF A 1 • ( BETA*A I F/ ( Ak*AF ) ♦ 0 • 28/ ( A*2*AK ) > «£F *RQFL/rtg 

1 7 7 


ALF A 1 

* ALFA1/(bEND.(C.2S/AK2 ♦ &ETA*AIF/AF ♦0«S*C/AK)> 

178 


ALF A I 

« ALF A 1 + ANU1 

1 7 9 


ALF A2 

* (o«25*(T + TN)*(C + o«S/aK))/Iq»S/AK2 + BETA ♦ A1F/AF +v.t» 

] 8G 


1 

*C / AK ) 

181 


ALBAR 

= l.o ♦ ALFA1 + ALFA1**2 ♦ 3.U*ALFA2**2 

1 62 


ZETAl 

= 1.0/ SORT ( ALBAR ) 

163 


2ETA2 

* ALFA2*ZETAl 

1 6 H 


SXX a 

ti*H**2/4»0 

18b 


EMFU 

a FTYT*(SXx/RO ♦ 0*2&*ZeTAI*(T**2 - TN*»2) ♦ It. r k2m ( T- 1 N ) *C 

1 66 

1 ODD 

FORMAT!* MFU* *.Ei6.8,* JN-LB/JN » / * ZETA1« *»tl6,8/ 

187 


1 

* 2 E T A 2 = * * E 1 6 » 8 / ) 

18B 


*.RITE<6,1000) EMFU. ZETAl ,ZETa2 

1 69 


RETURN 

1 90 


END 



iHDG PQ10 (IV) (PL0TF2) 


spkt.c Pc i o 
URPUK 2RHl-C3/|{)-l4:S3 
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PO 1 0 

(19 * 

(PL0TF21 

H3H60u*TPF*. P010 

1 


SUBROUTINE PL0TF2 (Dl ,B»h,RBiRFlL»T , MT .ALFAT.DELIAl 

2 


• , lPHASt,AP»ET,ANUT,AK »EF 

3 


• .HEAD ,KOPT ,P , PF) 

S 

5 

c 

c***« PLOT routine: FOR SUMMARY OF STRESS AND DEFORMATION ANALYSIS 

6 

C • * * * FOR LOR PROFILE FLAwOES 

7 

< . R . L£ 1 MB ACH , 19 DECEMBER 1972 

8 

c 


9 


DIMENSION X ( 1 DO 1 . Y ( 100> * IX ( loO ) . 1 Y ( 100 ) 

10 


# .HEAD ( 121 .KOPTtlQ) 

1 1 


DIMENSION A P ( 8 ) 

1 2 

c 


1 3 


CALL FRAMEV(O) 

IS 


CALL SCRECT ( 31 .31 .99 1 .991 > 

lb 


CALL PR 1 NT V < 72 .HEAD .S 1 * 1003 1 

16 


A 1 » 3.5 

1 7 


A2 = S.O 

1 8 


A 3 = 3 » 7 b 

19 


CALL XSCALV ( ” A 3 .A3 . u » 0 1 

20 


CALL YSCALV C - A i .A2.U.0) 

2 1 

c 


22 


SCALE-2 • 0 

23 


BS=a/SC ALE 

2 S 


H5=H /SCALE 

2b 


KFILShRFIL/SCALE 

26 


ts=t/scale 

27 


HTS » 2.0 

2B 


R l a 0 I / 2 « 

29 

c 


30 


X ( 1 1 ®bS/2 • 

31 


Y ( 1 >*-HS/2. 

32 


X ( 2 > »X ( l ) 

33 


Y 12 >=-Y ( 1 1 

3H 


I F ( KOPT ( 2 1 ,EQ .2 1 00 TO 2 10 

3& 


X < 3 1 *=-A ( 1 ) +TS + RF 1LS 

36 


Y ( 3 > “Y ( 2 > 

37 


X < S 1 «X ( 3 ) -RF ILS 

38 


Y < S 1 »Y ( 2 ) +RF ILS 

39 


(.0 TO 2 1 5 

so 


210 CONTINUE 

SI 


£ 1 aR 1 + a - R 8 

S 2 


ElS»El/SCALE 

S3 


X(3)»X(2)-2.*ElS 

SS 


Y ( 3 1 a Y ( 2 ) 

Sb 


X< S ) =X ( 3 ) - 2 * • K F ILS 

S 6 


Y ( S ) =Y ( 2 1 

S 7 


215 CONTINUE 

S3 


X ( 5 1 »X ( S 1 

H9 


Y(5)»Y(2)+HTS 

50 


X(6)*X(5)-TS 

51 


Y ( 6 1 =Y IS ) 

52 


X < 7 1 »X ( 6 1 

53 


Y ( 7 1 a Y ( 11 

5s 


X (8>*X l 7 1 

55 


Y(8)“Y(7)"»1 4-42 
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119) (PuOTF2) 


5 6 

A ( 9 1 * X ( 7 ) 

57 

Y ( 9 ) = Y < 7 ) -.6 

5 b 

X< IU)=X( 1 ) 

59 

Y < 10 ) » Y < 8 1 , 

A Li 

iinisim 

6 ] 

Y < 1 1 ) a Y ( 9 1 

62 

X ( 12)»X (7 ) 

6 3 

Y < 12 )=Y ( 7 >-.5 

6*) 

X < 1 3 ) =X ( 1 ) 

6 5 

YU3)*Y ( 12) 

6 6 

X 111) =X < 1 ) ♦ . 1 

67 

Y l 1 9 ) *Y ( 1 ) 

Ab 

X ( l 5 ) »X ( | ) ♦ « 6 

69 

Y< 15)»Y ( 1 ) 

?a 

X (1 6 ) = A ( 1 9 ) 

7 1 

Y < 16 )=Y (2> 

72 

X ( 1 7 ) =a < 1 5 ) 

7 i 

Y ( 1 7 ) = Y ( 2 ) 

7 9 

X< 1 8 ) * X ( 1 )♦.!» 

75 

Y ( 1 8 ) s Y ( I ) 

7 o 

X ( l 9 ) *X ( 1 8 ) 

77 

Y ( 19 ) =Y ( 2 ) 

78 

X ( 20 1 = A < S > 

79 

Y < 20 ) *Y ( 5 )+. 1 

8 j 

X ( 2 1 ) “ A ( 5 ) 

8 1 

Y ( 21 )»Y (5) +.6 

82 

X ( 22 ) »A ( 6 1 

63 

Y ( 2 2 ) * Y C 2 0 ) 

89 

X ( 2 3 ) a A ( 6 ) 

85 

Y ( 23 > = Y C 21 ) 

66 

X ( 29 ) »A 1 5 ) 

87 

Y ( 2 9 1 = Y ( 5 > + 0.6 

88 

X ( 2b ) *X ( S ) ♦ .2 

89 

Y < 2 b ) = Y (29) 

9a 

A ( 26 ) s x ( 6 ) 

9 I 

Y( 26 )»Y (29 ) 

92 

X (27 )*X( 6 ) -.2 

93 

Y ( 27 ) *Y ( 29 ) 

99 

X ( 28 ) s a ( 71-1.5 

95 

Y (28 ) «Y ( 1 2 > 

96 

X ( 29 ) *A ( 26 ) -. 5 

97 

Y(29)»Y( 12) 

96 

X ( 30 ) » A ( 29 ) - # 5 

99 

Y(3jl»Y( 12) 

) 00 

X(31 )»X(29)-.25 

iOt 

Y ( 3 1 ) «* Y ( ) 2 1 + l.o 

102 

X ( 32 ) * A ( 3 1 ) 

106 

Y ( 32 ) * Y( 12) - 0,5 

109 

X ( 33) =A( 2) ♦. 1 

105 

Y I 33 ) * Y ( S ) ♦* 1 

106 

X ( 39 ) » a ( 33 ) * 2 • 

107 

Y ( 39 ) a Y ( 33 ) 

1 06 

X < 3b ) a x ( 33) ♦1,0 

109 ■ 

Y ( 35 ) a Y ( 33 > 

1 10 

X ( 36 ) ax < 35 ) 

1 1 1 

Y ( 36 ) *Y ( 2 ) 
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(19) ( PLOtF 2 ) 


12 


A (3/1 

* X ( 6 ) -U.l 

1 3 


Y ( 3 7 ) 

* Y < 5 ) +0 • 1 

1H 


A < 3« ) 

3 X l 3 7 ) *» 2 » 

lb 


Y < 38 ) 

- Y ( 37 1 

16 


X < 39 ) 

« X ( 3 7 ) - 1 • 

1 7 


Y < 39 ) 

3 Y < 3 7 1 

1 d 


x ( ho) 

3 A { 39 ) 

19 


Y(Hu) 

- Y 1 2 1 

2l) 


xm ) 

*X < 3b > 

2 i 


Y < HI ) 

« Y U ) - . 1 

22 


X ( H 2 ) 

3 X< H 1 ) *2,0 

23 


Y ( H2 ) 

3 Y ( H l ) 

2H 


X ( H 3 ) 

3 X ( H l > + 1 . J 

2b 


Y(H3) 

* Y ( 1 ) 

26 


X ( HH ) 

=X(H3) 

27 


Y (HH) 

3 Y ( 2 ) 

28 

C 



29 


NP»HH 


30 

c 



31 


DO lo 

1 *1 .NP 

32 


CALL 

XSCLV 1 C X 1 I ) ♦ IX 1 I ) • l ERR » 

33 


CALL 

YbCLV 1 < Y ( I) , l Y i l ) , l ERR ) 

3 9 


1C continue 

38 

C 



3 e> 


CALL 

PLOTLN ( 1 , 2 , I X , 1 Y 1 

37 


CALL 

PL0TLN(2.3.IX, 1Y) 

38 


I F < KOP r ( 2 } ,EQ . 1 ) call PlOT«C(3.H 

39 


IF( KOPT ( 2) *E«*2 ) CALL PlOThC(3,H 

HU 


call 

PlOTLN < H ,5 , I X , 1 Y 1 

HI 


CALL 

PlOTLN < 5 ,6 . I X , l Y > 

H2 


CALL 

PLUTLNl6.7.IX.iY) 

H 3 


C ALL 

Pl 0TLN(7.1 < I X « i Y ) 

HH 


CALL 

PLOTLN < 8 ,9 . IX , 1 Y ) 

H S 


CALL 

PLOTLM 10, 1 1 . IX , IY 1 

H 6 


CALL 

PLOTLN ( 12*13. J X , 1 Y l 

H 7 


CALL 

PlOTAR(2.12.IX,1Y) 

H 8 


CALL 

PlOTAR (1*13* IX. I Y 1 

H 9 


CALL 

PLOTLNl l*»» IS. IX* IY j 

5U 


CALL 

PLOTLN (16. 17, IX, IY) 

SI 


CALL 

PLOTLN (I8.19.IA.1Y1 

S 2 


CALL 

PLOTARIH , la. IX . IY) 

S3 


CALL 

PL0TAR(3»19,IX.IY) 

5 H 


CALL 

PL0TLN(20.21 .IX.lYl 

SS 


CALL 

PlOTLN(22.23,Ia,IY) 

S 6 


CALL 

PL0TLN(2H,25,IX.IY) 

S7 


call 

PLOTLN ( 26,27 . I X , l Y 1 

Sd 


CALL 

PL0TAR(2.2H,IX.I Y) 

S 9 


CALL 

PlOTaRU ,2ft, IX. IY) 

6U 


CALL 

PlOTLN ( 1 2 .28 , I X , I Y ) 

6 l 


CALL 

PLOTAH <1,12, IX. IY) 

62 


call 

DbHLNV i I X ( 28 ) »IY(28) , J X ( 29 ) 

63 


CALL 

PLOTLN ( 29 ,30 . I X . I Y j 

6 H 


call 

DA8HLN131 ,32. » A , I Y ) 

6 S 


CALL 

PLOTLN ( 33 ,3H , I X . 1 Y ) 

6 6 


CALL 

PLOTLN ( 38 • 36 . IA , 1 Y > 

67 


call 

PlOTLN< 37 ,38 . IA, IY > 
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PG 


(19) (PLOTFZ) 


1 6 8 



CALL P00TLN(39.40,IA,lY) 


(69 



CALL PL0TLN1H1 »92,IA,IY) 


1 70 



CALL PlOTLN ( 93 ,99 . I X . I Y > 


1 7 1 

c 




l 72 



ISCALEal 


1 7 2 



X l * X ( 39) 


179 



Y l»Y (39) 


175 



Y 2 * Y 1 - . 0 5 


176 



1 0Y*5 


177 


20 

CALL YSCL V 1 ( Y 1 . 1 Y 1 , I E R R ) 


1 7a 



CALL Y8CLVI lY2.IY2.iERR) 


l 79 



STR£SS=9d« 


1 80 



1 YS=I Yl + IDY 


i a 1 



00 2 1 1 * ( .9 


182 



X2»X 1 


i a j 



CALL XSCLVI (XI , 1X1 , iERH) 


189 



CALL XSCLV 1 ( X2 , 1X2 . 1ERK 1 


1 85 



CALL L 1 NfcV ( I X 1 , I Y 1 . 1 X2 » I Y2 ) 


186 



|XS=IXl-8 


187 



CALL LABLV (STRESS, IXS. I YS.3. 1 .2 > 


1 88 



X 1 = a 1 - . 25 


189 



STkESS*STRESS-10. 


1 90 


21 

COnT INUE 


191 



I SCALE* I SCALE* l 


1 92 



IF< iSCALE.EO.2) GO TO 22 


1 93 



IF( ISCALE.EQ.31 GO TO 23 


1 9*1 



IF ( ISCALE.EQ.9 1 60 TO 29 


1 95 


22 

X 1 »X (37 > 


1 96 



Y l ■ Y ( 3 7 1 


197 



Y2*Y l - » 05 


l 98 



I D Y *5 


1 99 



60 TO 20 


200 


23 

X 1 *X ( 92 1 


201 



Y 1 *Y ( 92 > 


202 



Y2«Y1*.0& 


20 3 



I0Y*-2i 


209 



60 TO 20 


205 


29 

CONTI NUE 


206 

c 


0 I MENS ION 01(3). 02(31, 03(31, d9(3)»d5 (3) 


207 



.06(3) 

208 



• .0 7 (3), D8(J). 09(3). 010(31.011(3) 

209 



* .012(9) ,013(9) ,019(9) .015(9) 


2 1 0 



DIMENSION 0101(31 . 0102(3) . 010^(3) , 

0 1 G*t t 3 ) 

2 1 1 

c 




212 



0 A T A ( 01 ( I ) ,I«1 .31/18H PHASE 

/ 

2 1 J 



0 A T A ( 02 < 1 ) . !■ 1 .3 1 / 1 8 H BOLT FORCE 

/ 

2 1 9 



D A T A ( 03 ( 1 > , I *1 .3 ) /( 8H (LB) 

/ 

2 1 5 



U A T A ( 09 ( 1 ) , I = 1 .3 ) / 18H APPLIED MOMENT 

/ 

216 



OATAt Ob( I 1 , 1 • 1 . 3 )/ 18 H < ] N<*Lti 1 

/ 

2 1 7 



0 A T A ( 06( I ) .1*1 , 3 ) / 1 8H AXIAL OlSPLACEMT 

/ 

2 1 8 



0 A T A ( 07i 1 1 , I ■ l , 3 ) / ( BH (IN) 

/ 

2 1 9 



OATAt 03 ( I > , I *1 , 3 > / 1 8H ROTATION 

/ 

22u 



DATA) 09(i),I*)»3)/18H ( RAD ) 

/ 

221 



OA T A ( 0 10 < I > . 1* 1 . 3 ) / 1 8 H 


2 22 



DATA ( D1 1 ( 1) , I * 1 ,3 ) / i8H 


223 



DATA(012( I) . 1=1 .9) /29H STRESSES ON INSIDE / 
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PQlQ 

22 * 

22b 

226 

22 / 

22b 
229 
2 JO 
23 1 
2 32 
2 33 
238 

235 

236 

237 

238 

239 
28g 
28 1 
2 8 2 
283 
288 

285 

286 

287 

288 
289 
2b J 

251 

252 

25 3 
2 S 8 

255 

256 

257 

258 

259 

260 
261 
262 
263 
268 

265 

266 

26 7 
2 68 

269 

270 

271 

272 

273 
278 

275 

276 

277 

278 

279 


(19) (PL0TF2) 


C 


C 


C 


c 


c 


0ATA(D13(I ) » I « I . 8 ) /28H STRESSES ON OUTSIDE / 
0ATA(D18< 1 ) .I-1.8I/28H BENDING STRESSES IN F LG/ 
DATA < D 15 < I ) , I *1 .8 > /28H SUMMARY OF ANALYSIS / 


CALL PrtlNTV < 18 
CALL PR I NT V < 18 
CALL PRINTV( 18 
CALL PR JNTV ( 18 
CALL PR I NTV ( 18 
CALL PR I NTV ( I 8 
CALL PR I NT V US 
CALL PR I NTV ( 18 
CALL PN1NTVM8 
CALL L 1NEV ( l 90 
CALL L 1 N£ V ( 350 
CALL L I N E V ( 5 I 0 
CALL LlNEV (670 
CALL L I N£ V ( 830 


01.190.151) 

0 2 * 350. 151) 
03.350*131) 

08.510.151 ) 

05.510. 131 ) 

0 6 , 6 7 o , 15) ) 
07. 670.131 1 
08.830* 151) 

09 .830. 131) 

1 t 1 »V9) » 1 1 1 ) 
161 .350.31 ) 
161 .510*31 ) 

1 6 I .670 .3 i ) 

) 6 l .830.31 ) 


I XP»I X ( 38 ) *20 
I YP*1 Y (38 ) +30 

CALL PR1NTV128.D12.IXP.1YP) 
I XP» I X ( 33 I *20 
l YP 3 I Y ( 3 3 ) +30 

CALL PR1NTV128.013.1XP.IYP) 
1 X P « J x ( 8 1 ) + 20 
1 YP» I Y ( 8 1 ) -80 

CALL PR I NTV l 28 .0 I 8 , 1XP . I YP ) 
CALL PRINTVC28. 015. 800.961 ) 
CALL LINEV(800. 951 .568.951) 
CALL LINEV(800.987.568,987) 


DIMENSION DATA! I l ) 

DATA (DATA l (I) . I«1 . 1 J /6H DJa/ 


CALL PlOTLB ( 28 .0 .5 , 1 X . I Y ,D 1 ) 

CALL PL0TTX< 28,3 2,5, 1X.IY.DAT A l ,6) 
CALL PLOTLB ( 1 8 . 1 0 .25 , IX , I Y ,H) 

CALL PL0TLB<9,l0.-2i»,lX,IY,B) 

CALL PL0TLB(21 . 10. lu. IX.IY .T) 


BF C 3 A P ( 1 ) 

E AM» AP ( 2 ) 

ADFL» AP ( 3 ) 
m P 0 3 A p l 8 ) 

CH I 0 s AP ( S 1 
ENQ* AP ( 6 > 
EMO«AP ( 7 > 

EQO = AP l 8 ) 
SSCALEs. 000025 
rtSC ALE« JO. 

RO"R I * 7 / 2 . 

AL 3 NTS *2.0 
OL *o . u5 
NL* AL/DL 
I 3 I PHASE 
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POJO (IV) (FL0TF2 ) 


2eo 
28 1 
282 

28 J 
2 8 9 

285 

286 

287 

288 
289 
2 9U 

29 1 

292 

293 
299 
29b 
296 

2 97 

298 

299 

300 

301 

302 

303 

3 0 9 

305 

306 

307 

308 

309 

310 

31 1 
3 1 2 
313 

319 
3 1 5 
3 1 6 
3 1 7 
3 1 8 
3 l 9 

320 

32 1 

322 

323 
329 

325 

326 

327 

328 

329 

330 

331 

332 

333 
339 
335 


AK2* AK*AK 

BENQ<*cT*T«* 3/< 1 2 • • < i.-ANUT**2) ) 

CWF».5/(AK2«AK«'8EnD) 

00 100 K * 1 , N L 
£K«K 

AKa ( EK-1 . ) *UL 

S X a S J N ( A K * X K ) 

CXaCOS ( AK*XK ) 

EXaEXP ( -AK*XK ) 

WX«CrtF*EX* (EOO*CX*An*EMo* ( C X »5 X ) )+HlP0 
EHX= Ex*(EMO*ICa*>SX>-EQo*5X/aK) 

E OX * Ex* (E8 u*1Sa»CX )+2**AK*£mO*SX) 

IF (K.EO. 1 ) ftTQPaftX 
E NX *E No 

£NY=ET*T*rtX/RO 

!F( 1PhaSe»E0.2J ENYbenY-eT*T*ALFAT*UELTaT*P*R0*ANUT/2 

IFUPHASE.E0.3l ENY»ENY+P*RO*ANur/2* 

IF(1PHASE.E0.9) EnY»EnY+PF#P*R0*ANUT/2, 

ErtYa ANUT*£MX 

IF (KOPT (2 ) «EU. 1 ) (iO TO 202 

SXla-b»*EMX/T**2*ENx/T 

SX0“*6.*£MX/T**2*£NX/T 

SYI«-6.*EMY/T**2+ENY/T 

5roa+6.*EMY/T**2+ENY/T 

T M A X « 1 «S*EQX/T 
GO TO 205 

202 IF(XK.LE.RFIL) GO TO 203 
T X a T 

GO TO 209 

203 XPaRFIL-XK 

Y P*5QR T <RFIL**2-Xp**2) 

T X a T *RF I L-YP 
209 CONTINUE 

SX1“-6»*£MX/TX**2+ENX/TX 
5X0“ 6, •EMX/TX**2+EnX/TX 
SY1*-6.*£MY/TX**2+EnY/TX 
SYO= 6.*emY/TX**2+EnY/Ta 
TMAX a) »5*EQX/TX 
205 CONTINUE 

YS=Y (2)+XK/SCALE 
X5X1 ax ( 3V ) ♦SSCALE-SX 1 
XSY |«X ( 39 1 +SSCAlE*SY I 
XSXOaA ( 35 ) +SSCALE*SXO 
XSY0«X(35)+SSCA0E*SY0 
XTMAX=XC6)-SSCALE*TMAX 
A ft X a X ( 6 ) •*■ ft SCALE* ft X 

CALL YSCLVl (YS, 1YS.1ERR) 

CALL XSCLV1 (XSXI ,1X5X1 ilERR) 

CALL XSCLV1 (XSYi , IX6YI ,l£RR) 

CALL XbCL V 1 l XSXo . I X5X0 . 1 ERR 1 
CALL XSCLV1 IXSYO, iXsYO* IERR) 

CALL XSCL V 1 < XTMaX , I XTMAX , I ERr 1 
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p 0 1 u 

336 

337 
3 3 8 
33V 
340 

34 1 

342 

343 

344 

345 

346 

347 
34B 
349 
3SQ 

35 1 

352 

353 

354 

355 

356 

357 
35« 
359 
36u 

36 1 

362 

363 

364 

365 

366 
367 
366 

369 

370 

37 1 

372 

373 

374 

375 

376 

377 

378 

379 
3SU 

38 1 
3 82 

383 

384 

385 

386 

387 

388 
389 
39Q 
391 


1191 { PL 0 TF2 ) 


C 


c 


c 


CALL XSCL V 1 tXrtX.lXWA,jEKR) 

IF CK.EU. 1 ) GO TO 98 
I F (K ,ty.6» GO TO 98 
lFU.Ey, 1 U GO TO 96 
GO TO 99 

98 CALL PLOTVl IXSXJ »IT5,29,0> 
CALL PLOTV C IXST I , 1 YS.30.Q) 
CALL PLQTVl IXSX0,lY5,2?,0) 
CALL PLOTVUXSTo. ITs.30,0) 
CALL PLOTV I IXTMaX.XTS.2S ,0) 
CALL PLOTV! IXWX, I YS,28,0> 

99 CALL PLOTV C IXSXI , l YS.3&.0) 
CALL Pi_OTV< IXSYl ,1Y5,35,0) 
CALL PLOTV ( IXSXO, I YS. 35,0) 
CALL PLOTV U XSYO . 1 YS , 35 ,0 ) 
CALL PLOT V < 1 XTMaX , J TS ,35 ,0) 
CALL Plot V < I XWiX , 1 YS .35 ,0 ) 

100 CONTINUE 


* 8 0 T « T 0 P - H * c H 1 0 
XAiBOT »X ( 7 ) ♦WBOT 
AVKTOP«X < 7) +WTOP 
SlGT0P*EF*WT0P/«0 
SI GBOT»EF»WBOT/NO 
YSTOP* Y < 2 ) 

YSBOTsYW) ' 

XST0P*XCH3)+SSCALE*SI6T0P 

XSB0T=X( 43) +SSCALE*S1GB0T 

CALL XSCLV 1 t XSTOP , I ASTOP , JERr ) 

CALL YSCLV 1 ( YSTOP. I YSTOP , JERR ) 

CALL XSCLVl (XSBOT, 1XSB0T , JERR ) 

CALL YSCLVl (YSBOT,IYSB0T,IERR) 

CALL XSCLV1 IX»B0T. IXWbOT.IERr) 

CALL X5CLV l ( XWTOP . I XWTOP , I ERR ) 

CALL PLOTV( 1XST0P, 1YST0P,30,0> 

CALL Plot V { I XSBuT , l TSbOT , 30 ,0 ) 

CALL LINEV1IXC44) ,1 Y{44) , IXST0P.1YST0P) 

CALL L 1 NE V ( l X ( 43 ) , |V( 1 ) » I XSBoT , I YSBOT) 

CALL L 1 NE V { 1 X S T 0 P , I T S T OP , I XS 8 0 T , 1 Y SBOT 1 
CALL PLOTV ( IXWBOT , I T (1 > .28 ,0> 

CALL DSHLNV ( I XWBOT , 1 Y ( I ) , I XWTOP , X Y < 7 ) ,S, 15) 


0 A T A ( 0 J o 1 l I ) * 1*1 ,3 ) / 1 8 H BOLT-UP / 

0ATA(0102 l I) » 1 » 1 ,3 ) / 16H START-UP / 

DATA (0103 ( I) , I® 1 , 3 ) / 1 8H OPERATION / 

0 AT A < 0 104 U ) , l « l , 3 > / 1 8H SHUT-DOftN / 


lFUPHASE.EO.il CALL PR 1 N T V < ) 8 , 0 1 0 1 » I 9 0 . 7 6 ) 
I F ( I PHASE • EQ • 2 ) CALL PR 1 NT V < 1 8 . 0 l 02 » l 9 q , 76 ) 
IF ( 1PhASE.E0.31 CALl P« 1 NT V ( 1 8 ,0 103 . I 9Q , 76 ) 
1F« lPHASE.Eu.4) CALl PR i NT V ( l 8 . U 1 04 . 1 90 , 7 6 ) 

CALL L ABLV < BFC , 38p * 76 ,- 6 , 1 , 1 > 

CALL LABLV < EAM , 540.76 . - 6 , 1 , 1 ) 

CALL LAdLVlAOFL. 700.76. -6. 1 . 1 ) 
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1191 ( P L 0 T F 2 ) 


392 

CALL LABtVtCHltl. 860.76 

393 C 


3 9 S 

RETURN 

395 

END 

®HD(i D 2Q0 

(20) (OUTAN) 


BPRTeC L»20u 

PURPUR 29HJ-Q3/ 1G-JH;S3 
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2U0 ( 20 ) (OUTAN) 

H3H6uo*TPf S » D 200 

1 SUBROUTINE OUTAN (HEAD .A »SR .5 ) 

2 C 

3 <:*••• PRINT-OUT ANALYSIS KESUtTS 

‘t C KtR.LElMBACH, H JANUARY 1973 

5 C 

6 DIMENSION A < 9 , H> , SR ( 5 , H J , S < S , H > , HE AO ( I 2 > 

7 C 

B 101 FORMAT(lHl) 

* 102 FORMATI 1 2 A 6 ) 

10 103 FORMAT(/t OUTPUT OF THE STRESS AND DEFORMATION ANALYSIS RfcSuLTS«/‘ 

H 1QH FORMAT!/* VARIABLE *,30*.» BoLT-UP * 


1 2 

♦ 



,» StAKI-UP * 


1 3 

* 



,» OPERATION * 


IH 

♦ 



,» Shut-down *//) 

IS 

201 

FORMAT l 1 6 H 

BOLT FORCE <LB) ,2mX,HE16.8/ 


1 6 

* 

29H 

EQU 1 V 

APPL MOMENT 

( 1 N-LB/ IN ) , 1 1 X ,HE 16.8/ 

1 7 

♦ 

22H 

AXIAL 

deflect I ON 

( IN) , lax ,HE16.8/ 


i 8 

• 

23H 

RADIAL 

DEFLECT I ON 

(IN) » I7X .HE16.8/ 


1 9 

• 

2 5rt 

flange 

ROTATION (RADIAN) ,1SX,HE16. 

8/ 

20 

• 

1 8h 

bolt STRESS IPSI) 

* 22 A ( HE 1 6 * 8 / 


2 1 

» 

2 Qrt 

GASKET 

STRESS ( PS 1 ) .2UA,HEI6,8/ 


22 

* 

27H 

STRESS 

In FLANGE 

TOP (PS1 ) . I3X .HEl6*6/ 

23 

* 

3CH 

STRESS 

In f Lange 

BOTTOM ( P S I ) , 1 OX , 

HE 1 6 • 8 / ) 

2H 

2 02 

FORMAT ( HCh 

STRESS 

RcSULTaNTS 

ILB/IN) 

NX«,Ht|6.8/ 

2 S 

* 

HOH 



(LB/IN J 

«V* t He 1 6.8/ 

26 

* 

HOH 



1 l N-LB/ 1 N > 

MX" » HE 1 6 * 8/ 

27 

« 

HOH 



( IN-LB/IN) 

MY" .HE16.0/ 

28 

♦ 

HOH 



(LB/IN) 

UX*»HEl6*8/) 

29 

203 

FORMAT ( HOH 

STRESSES AT NeCK 

(PS1) inner 

SUX»,He 1 8.8/ 

30 

* 

HOH 



t PS I > INNER 

SlGY* i HE 16.8/ 

3 1 

* 

HOH 



(PSI) OUTER 

St GX" • HE 16 .8/ 

32 

0 

HOH 



(PSI) OUTER 

S I g Y s i *• L 1 6 • 8/ 

3 3 

0 

HOH 



(PSI) MAX 

TAU B »H£16.8/ ) 


3H C 

3 b WRITE(6. 101 ) 

36 lf«R I TE ( 6 , lu2 » ( HEAD ! 1 > , I = I , 1 2 ) 

37 IAR1TE(6, 103) 

38 wRlT£<6,luH) 

39 WRITE (6, 2o lH<A<l.J),J*l»H>» 1*1,9) 

HO ARnEl6,20 2M(SftlI,J»,J»|.H),l*l«5) 

HI wRlTE(6,203)<<S<l.J),J a l.H)»I*l.S) 

H 2 C 

Hi RETURN 

HR END 


or j N 
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Section 5 

NUMERICAL EXAMPLES 


In this section one example is presented that has been computed by- 
hand. Corresponding computer results of this and of additional examples 
are also given. 


5.1 EXAMPLE: STEEL FLANGE WITH STEEL GASKET 


Given: 

Safety Factors: 


Tube and Flange 
Material: 

Bolt Material; 

Gasket Material; 


Nominal pressure 

Nominal diameter 

Tube thickness 

Proof 

Burst 

General 

Gasket 

347 SS steel 
A286 Steel 
CRES 321 -A 


p = 1500 psi 
d. = 8.00 inch 
t = .438 inch 


1.5 

2.0 

1.5 

2.0 


The material data are given in Tables 2-10 and 2-11. Following the outline 
in Appendix A, the following results are obtained: 


(a) Tube thickness given as t = 0.4375 inch 

The thickness based on Eq. (2.23) would be 


t 


1.5 x 1.5 x 10 x 4.000 


35 x 10" 


9.0 x 10 3 
35 x 10 3 


0.257 inch 


For a more accurate thickness computation Eq. (2.18) 
should be used, giving 
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6,0 x 10 3 

23.3 x 10 3 - 0.75 x 10 3 


For higher internal pressures this difference is more 
distinct. 

(b) Initial guess of bolt size 

dg = 0.4375 inch, size number 4 (Table 2-3) 

Machined spot faces 

e l = 1.03 x 0.4375 = 0.4506 inch 
e 2 = 0.91 x 0.4375 = 0.3981 inch 

Hole diameter 

d ho ie = 0.437 + 0.005 = 0.442 inch 

Spot face diameter 

d . = 2 x 0.4506 = 0.9012 inch 

spot 

Fillet radius for spot face 

r , = 0.062 inch 

spot 

(c) Bolt circle radius 

r B = 4.00 + 0.4375 + 0.062 + 0.4506 = 4.9501 
On the plot appears the diameter of the bolt circle 

(diam)g = 2 x 4.9501 = 9.900 inch 

(d) Flange width 

b = 4.950 + 0.3981 - 4.000 = 1.348 inch 
On the plot appears the outer diameter of the flange as 

(diam)p^ = 8.000 + 2.696 ~ 10.696 inch 
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(e) Gasket width and gasket radius 
Estimate for gasket radius 

r Q = 1/2 (4.950 - 0.221 + 4.000) 
= 1.2 (8.729) = 4.3645 inch 


Gasket width, calculated on the assumption that the gasket 
is initially stressed to the yield strength K^, but under proof 

pressure it is allowed to the pressure dependent seating stress 

a G = k p P * see A * 5 ( b )) 


1.5 x 1500 x 4.3645 


2 [l.O x 40.0 x 10 3 - 1.0 x 1.5 x 5,50 x 1.5 x 10 3 x 2.5] 

x 10 * 9 J2 ,, x 10 3 _ 0 .322 inch 

80.0 x 10 - 49.5 x 10 3 30.5 x 10 3 


The gasket will be located close to the bolts, allowing a 
tolerance of c 2 = 0.05 inch 


r G = 4.950 - 0.221 - 0.161 - 0.050 = 4.518 inch 

The inner radius of the gasket is 

r Q = 4.518 - 0.161 = 4.357 inch 
i 

and the corresponding diameter appearing on the plot 
(diam) G = 2 x 4.357 = 8.714 inch 

i 

The outer radius of the gasket is 

r Q^ = 4.518 + 0.161 = 4.679 inch 
and the corresponding diameter 

(diam) G ^ = 2 x 4.679 = 9.358 inch 
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(f) Required bolt force 

P^ L) = Zn x 4.518 x 0.322 x 1.0 x 40.0 x 10 3 
= 6.28 x 1.455 x 40.0 x 10 3 
= 9.1374 x 40.0 x 10 3 = 365.5 x 10 3 lb 

Pg ^ = 7 r x (4.518)^ x 1.5 x 10 3 x 1,5 

+ Zir x 4.518 x 0.322 x 1 .0 x 1 .5 x 1,5 x 10 3 x 5.50 x 2.0 
= 3.14 x 20.41 x 2.25 x 10 3 + 6.28 x 1.455 x 4.5 x 10 3 x 5.50 
= 63.24 x 2.25 x 10 3 + 9.075 x 24.75 x 10 3 
= (142.3 + 224.6) x 10 3 = 366.9 x 10 3 lb 

(g) Number of bolts 

_ 366,9 x IQ 3 

131 x 10 3 x 0.10631 

= 366 ' 9 X = 26.3 ~ 26 bolts 

13.93 x 10 

(2.0/1. 5) 366.9 x 10 3 ~ 488 x IQ 3 

200.0 x 0.10631 ~ , in 3 

2 1 .3 x 10 

= 22.9 ~ 23 bolts 
= 26 bolts are required. 

(h) Bolt spacing 

2a- 4.518 _ 6.28 x 4.518 
S 26 26 


28.37 

26 


= 1 .09 inch 




n 


B1 


n 


B2 
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Minimum allowable spacing is 

s . = 1.81 x 0.4375 = 0.792 inch < 1 .09 inch 

min 

Maximum allowable bolt spacing 

s « 8 x 0.4375 = 3.500 inch >1.09 inch 
max 


{i) Flange height 

Internal lever arm 

e = 4.950 - 4.518 = 0.432 inch 

Radius of the shell middle surface 

r = 4.000 + 0.219 = 4.219 inch 
o 

Thickness required to carry axial force 


t N = 0.438/2 = 0.219 inch 

Untimate moment to be carried 

1.5 x 366.9 x 10 3 x 0.432 
m Fu = 2?r 4.219 

550.4 x 10 3 x 0.432 
6.28x4.219 


237,7 x 10' 
26.5 


= 8.97 x 10 3 in-lb/in. 


Effective flange width 

b = 1.348 - 0.442-^0.442/1.09 
= 1.348 - 0.442 -^0 .406 
= 1.348 - 0.442 x 0.637 
= 1.348 - 0.282 = 1.066 inch 
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Assume 

i l = 0.8 

C 2 « 0.18 

The coefficients of the quadratic equation for h are 

a , r ,«3 1.066 37.31 , „3 , „„„ , „3 

A = 35x10 7 , _ . a - . x 10 = 1.474 x 10 

6 x 4.219 25.31 

■D n I- , n 3 n ,o 0,438 - 0.219 op o rti n-7 , n 3 

B = 35 x 10 x 0.18 x 2 ~ 35 x 0.0197 x 10 

= 0.6895 x 10 3 

C = 35 x 10 3 x 0.8 x ^- 438 > ~ t 0 . * 21 ?)- - 8.97 x 10 3 
= 35 x 10 3 x 0.2 x (0.192 - 0.048) - 8.97 x 10 3 
* (1.008 - 8.97) x 10 3 = -7.962 x 10 3 

R 2 = [(0.6895 ) 2 + 4 x 1.474 x 7.962jx 10 6 
= [o.4754 + 5.896 x 7.962 jx 10 6 
= [o.4754 + 46.944o]x 10 6 = 47.42 x 10 6 

R = 6.886 x 10 3 


h = 


6.886 - 0.690 
2 x 1.474 


6.196 

2.948 


= 2.102 inch 


If the contribution of the plastic hinge is neglected, i.e, 
if = C 2 = 0 is assumed, then 


A = 1.474 x 10 3 
B = 0 

C = -8.97 x 10 3 
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R 2 = 4 x 1.474 x 8.97 x 10 3 
= 5.896 x 8.97 x 10 3 
= 52.89 x 10 3 

R = 7.273 

h - m - 2 - 467 inch 

This is 0.3 inch more than the previous result. 

The same result would have been obtained by talcing the 
old formula from Ref. 1, 

h _ / 6 x 4.219 x 8.97 x IQ 3 

if 35.0 x 10 3 x 1.066 

V 25.31 x 8.9T _ / 227.0 

37.31 ^37.31 

= -^6.084^ = 2.467 inch. 


The weight savings accomplished by considering the plastic 
hinge is therefore approximately 10%. 


(j) Flange weight 

Weight area 

A = (1.348 - 0.438) x 2.102 
w 

= 0.910x2.102 - 1.913 in 2 


Centroidal radius 


r 

w 


4.000 


. 0.438 + 1.348 
+ 2 


= 4.000 + 0.893 = 4.893 inch 
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V olume 

Vol = 2 it x 4.893 x 1.913 = 58.78 in. 3 

Weight 

Aw = 0.288 x 58.78 = 16.9 lb 

f 

On Fig. 5-1 the design geometry is summarized. Figures 5-2 and 5-3 

show the stresses and radial displacement at initial torqueing and at proof 

pressure, respectively. The axial stresses (7^ are indicated by the curve 

labeled "X, " the circumferential stresses Oy by "Y" and the transverse shear 

stresses, T , by "T. " The radial displacement is shown as "W. " At the 
xz 

bottom of the plot the total bolt-force and the applied moment in (in-lb/in), as 
well as the axial displacement and the rotation of the flange are given. A sample 
printout is given for verification. 

5.2 EXAMPLES FOR WEIGHT COMPARISON WITH CONVENTIONAL FLANGES 

Before some weight comparisons with conventional flanges are made it is 
instructive to discuss a series of designs computed by the program. This series 
points out the need for judgement in the selection of the design parameters and 
materials. 

Figures 5-4 through 5-6 show a flange which was designed to meet the 
algorithms for minimum tube wall thickness and minimum gasket width. Pos- 
sibly the tube wall thickness is less than the minimum gage requirements for 
handling and accidental impact loads. The gasket width should be selected to 
fill out the available space between the inside of the tube and the inside of the 
bolts, including some dimensional tolerance. Possibly a thinner gasket should 
be designed. An algorithm is available in the program to automatically compute 
the gasket width to make use of the available space. Figure 5-4 shows a design 
with 6 bolts, which is the minimum. Figures 5-7 through 5-9 show a design 
similar to the previous one with double the pressure. In both this and the 
previous designs the stresses are well below the allowable ones. This is due 
to increased flange height based on bolt spacing allowing h not to be less than 
s/3. This requirement is based on experience since it is difficult to assess 
it analytically. The intent is to avoid waviness of the flange between the bolts. 
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test flange ISOO PSI. 8 IN 01 A 


PAE5SURE ISOO PSI 6 
IEIGHT OF FLANGE 16.979 LB 



0.HA70TA, 26 MIXES 

0.9010IA SAOTPACE 
0.0629 FILLET 


FLANGE FlATEAIAL 3975S AT AT FTYs 35.000KS1, FTU= 90.000KS! 
BOLT flATEAIAL A286 AT AT FTY=I3l .OOOKSt » FTU=200.000KSI 
GASKET FlATEAIAL CAES321-A SEATING STAESSs 18.900KSI 

VIELO STRENGTH* 10. OOOKSt 
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Fig. 5-2 - Stresses at Initial Torquing 


X =a x axial stress {ksi) T = transverse shear stress (ksi) 

Y = C(p circumferential stress (ksi) W = radial displacement (10-fold 

magnified) 
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lt$T FLANGE 1500 PSt, 8 IN 01 A 



Fig. 5-3 - Stresses at 1500 psi (proof pressure) 
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f LANGE PBfiflPtTHIC CASE 


DIAs AIN PRESS = 100 PS1 


PRESSURE 100 P51G 
tftGHT OF FLANGE O.MA LB 




Fig. 5-4 - Flange 1, Design 
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FLANGE PARAMETRIC CASE l DIA= AIN PRE5S = 100 PSI 






5-13 









LMSC-HREC TR D306492 



5-14 











LMSC-HREC TR D306492 


FLANGE PARAMETRIC CASE 2 DIA= AIN PRESS = 200 PSI 


0.I25R 

0.25501 A. 0 HOLES 
0.5S00I A SPOTFACE 
0.062R FILLET 


06? 


FLANGE MATERIAL AL606IT6 RT FTV= 35.000KS), El 0; H2.000KSI 
BOLT MATERIAL AL202AT3 RT FTVs 50.000KSI. HU* 62.000KSI 
GASKET MATERIAL ASBEST1/16 SEATING STRESS? 3. TOOKSI 

VIELO STRENGTH? 10.000KSI 


Fig. 5-7 - Flange 2, Design 


PRESSURE 200 PSIG 
hEIGHT OF FLANGE 0.700 LB 
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fLftwtt PARAI tTfilC CASE 2 0[A= SIN PRESS = 200 PSI 

SUnrnARv op analvsis 



PHASE BOLT FORCE APPLIED FUnENT AXIAL OISPLACEP1T ROTATION 

<LB) <!N-L6> (IN) <RA0> 

BOLT-UP S.SA600XIO* 03 I.62969XI0* 02 O.OOOOOXIO" 01 B.ABSmiO -03 

Fig. 5-8 - Flange 2, Stresses at Initial Torquing 
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FLANGE PAfiAfnETRIC CASE 2 01 As AIN PRESS s 200 PSI 
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Figures 5-10 through 5-12 show a design in which the required gasket 
width is controlling the width of the flange. A decrease in the minimum 
seating stress at proof pressure would reduce the required gasket width. 

A different gasket should be used in this case. Since the flange height was 
selected based on the strength requirements the peak stresses as shown on 
Figs. 5-11 and 5-12 are close to the allowable ones. 

Figures 5-13 through 5-15 again show a flange design controlled by the 
bolt- spacing -to-height ratio of 1/3. Consequently the stresses are low. Figures 
5-16 through 5-18 show a flange with more balanced proportions. It has the same 
inner diameter as the previous one but the pressure is doubled. The third flange 
with this diameter is again controlled by the gasket width requirements (see Figs. 
5-19 through 5-21). For this flange a different gasket should be selected. 

The flange shown on Figs. 5-22 through 5-24 is well proportioned and 
the stresses are well under the allowable stresses, although here as before 
bolt spacing is the controlling factor. Figs. 5-25 through 5-26 show a flange 
designed for the same inner diameter but twice the pressure. This is a 
strength-controlled design. 

Finally, Figs. 5-28 through 5-30 show three typical low profile designs. 

The last two are again partially controlled by the width of the gasket although 
only slightly. 

Table 5-1 presents a comparison of flanges designed with conventional 
and low profile contours. The weight savings are impressive even using the 
unfavorable configuration with the gasket located toward the inside of the tube. 
Figures 5-31 through 5-45 present plots of the low profile flanges with the 
conventional contour indicated by a dashed line and shading. The saving in 
space requirements is obvious when the outer diameters of these flanges are 
compar ed . 
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FLANGE PARAMETRIC CASE 


Dl As AIN PRESS = 500 PSI 


PRE55URE 500 PStG 
LEtGHT OF FLANGE 0.503 



0.25501 A. 12 MOLES 
0.580DI A SPOTFACF 
0. 062R FILLET 


FLANGE MATERIAL AL6051T6 RT FTV= 35.000K5I, FTU= A2.000KS1 
BOLT MATERIAL AL202AT3 RT FTYs 50.000KS! , FTUs 62.000KS! 
GASKET MATERIAL ASBEST01/14 SEATING STRESSs 3. 700K5I 

YIELD STRENGTHS 10.000KSI 
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SUWftRY OF ANALYSIS 


STRESSES ON INSIDE 


STRESSES ON OUTSIDE 


-W. -30. -20, -10. 0, 10. 20. 30. MO 

BENDING STRESSES IN FLG 


BOLT FORCE 


APPLIED (WENT AXIAL OlSPi ACEN1 ROTATION 

(IN-LB) (IN) (RAO) 


1,90075X10 


H. 60911X10 


0.00000X10 


1.16959X10 
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FLANGE WAT^RI AL ALG061T6 AT FTVs 35.000KSI, FTU: AJ.OOOKST 
BOLT MATERIAL ALZ03AT3 AT FTY = 50.000KSI, F|U= W.OOOKS! 
GASKET WATER! AL ASBESTt/lS SEATING STRESS s 3.700KS1 


YIELD STRENGTH: IO.OOOKSI 



Fig. 5-13 - Flange 7, Design 
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Fig. 5-14 - Flange 7, Stresses at Initial Torquing 
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FLANGE PARAMETRIC CASE 7 OIA= SIN PRESS = 100 PS1 



Fig- 5-15 - Flange 7, Stresses at Proof Pressure 
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PRESSURE 200 PS1G 
iiEIGHT OF FLANGE 1.120 U 



0.1 25R 

0.25501 A, 13 MOLES 
0.58001 A SPOTFACE 
0.062R FILLET 



FLANGE MATERIAL AL6081TG RT FTVr 35.000KSI * FTU= W.OOOKSI 
60LT MATERIAL AL2025T3 RT FTYs 50-000KSI. FTU= M.000KS1 
GASKET MATERIAL ASBEST0U16 SEATING STRESS^ 3.700KSI 

VIELO STRENGTHS IO.OOOKS1 
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FLANGE PARAMETRIC CASE 6 D1A= 81 N PRESS = 200 PS1 





Fig. 5-18 - Flange 8, Stresses at Proof Pressure 
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Fig. 5-20 - Flange 9, Stresses at Initial Torquing 
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SUWWRV OF ANALYSIS 


STRESSES ON JNS10E 


ir 


STRESSES ON OUTSIDE 



*10. *30. -20. -10. 0. 10. 20. 30. NO. 

BENDING 5T BESSES IN FLG 



BOLT fOHCE 


APPLIED ItrnENT AXIAL DISPLACED ROTATION 

UN-LB) (IN) (RAD) 


T. 52752X10 


1.03511X10 


1 .13665X10 


I.5N919XI0 
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PLWG£ PARAMETRIC CASE 13 0IA=12IN PRESS = 100 PSI 


PRESSURE 100 PSIG 
kEIGHT OF FLANGE 2.37R LB 



FLANGE MATERIAL Al 606 IT 6 RT FTVs 3S.000KSI, FTU- A2.000KSI 
BOLT MATERIAL AL202RT3 RT FTV= 50.000KSI, FTU= 62.000KSI 
GASKET' MATERIAL ASBE5T1/16 SEATING STRESS^ 3. 700KSI 

VIELO STRENGTHS 10.000KSI 


Fig. 5-22 - Flange 13, Design 
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FLANGE PARACETRIC CASE 13 0IA=12IN PRESS = 100 PS1 


SUWIARV OF ANALYSIS 



PHASE 

BOLT FORCE 

APPL1E0 (WENT 

AXIAL OISPLACem 

ROTATION 


(LB) 

UN-LB) 

(IN) 

(RAO) 

BOLT-UP 

i.msmxio 40 * 

I.17I77X1Q* 02 

0.00000X10* 01 

6.MYIOOXIO -03 


Fig. 5-23 - Flange 13, Stresses at Initial Torquing 
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FLANGE PARAItTRIC CASE 13 01A=I2IN PRESS = 100 PSI 


SUWARV OF ANALYSIS 


STRESSES ON INSIDE 




STRESSES ON OUTSIDE 




BOLT FORCE 



-NO. -30. -20.-10. 0. 

10. 20. 30. 90 

J 

BENDING STRESSES 

IN FL6 

n 

0.623 



APPLIED CDCIENT 

J AXIAL DISPLACED 1 

ROTATION 

t IN-LB) 

1 UN) 

(RAO) 


OPERATION 


2.90635X10 


2.33837X10 


I.7VW9X10 


1.29533X10 
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FLANGE PARAMETRIC CASE U DIAxl21N PRESS - 200 PS! 


PRESSURE 200 PSIG 
bCIGHT OP FLANGE 2.030 



0.1 07R 

0.25501 A. 28 HOLES 
0.58001 A SPOTFACE 
0.062R FILLET 


FLANGE MATERIAL AL8061T6 RT FTY- 35.000KS1, MU= R2.000KSI 
80LT MATERIAL AL202RT3 RT FTY= 50.000K51 , FTU- 62.000KS* 
6ASKET MATERIAL ASBEST01/18 SEATING STRESS- 3.700KSI 

YIELD STRENGTH = 10.000KSI 



5-34 




LMSC-HREC TR D306492 


FLANGE PARAPETAtC CASE JH 0[A=12IN PRESS = 200 PSI 



Fig. 5-26 - Flange 14, Stresses at Initial Torquing 
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FLANGE PARAMETRIC CASE 19 0IA=22IN PRESS = 100 PS! 


PRESSURE 100 PSIG 
WIGHT OF FLANGE 3.94! LB 



FLANGE ARTERIAL AL60MT6 RT FTY= 3S.OOOKS1, Fill- 92.000KS! 
BOLT MATERIAL AL2029T3 RT FTY= 50.000KS1, FTU- 62.000KSI 
GASKET MATERIAL A58EST1/14 SEAUNG STRESS- 3.700KSI 

YIELD STRENGTH' tO.OOOKSt 


Fig. 5-28 - Flange 19, Design 
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FLANGE PARAMETRIC CASE 20 OJA-221 N PRESS = 200 PSI 
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FLANGE PARAMETRIC CASE 25 DIArRSIN PRESS = 100 PSI 


PRESSURE 100 PSIG 
hEIGHT OF FLANGE 18.686 LB 



FLANGE MATER1 AL 
BOLT MATER! Al 
GASKET MATERIAL 


AL606116 RT 
AL202RT3 RT 
ASOESTI/lfc 


FTV: 35 . 000KS1 ■ FTU: A2.000KSI 
FTY- 50.OO0K5! , FTU= F2.000KS1 
SEATING STRESS: 3.700KSI 

Y1EI.0 STRENGTH: 10.000KSI 


Fig. 5-30 - Flange 25, Design 


250R 

.25501 A, 179 HOLES 
.58001 A SPOTFACE 
. 062R FILLET 
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Table 5-1 - Comparison of Flanges 
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Fig. 5-31 - Flange Comparison 1 

This number refers to the list on Table 5-1. 

Legend: The contour of the conventional flange is 

indicated by dashed lines and shading. 
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LIQUID OIVGEN HO PSI, 12 IN DM 


PRESSURE HO PSI6 



O.J12R 

O.H»tA, 56 MOLES 
0.68I0IA SPOTFACE 
0.062R FILLET 


FLANGE MATERIAL At 
DOLT MATERIAL AL 
GASKET MATERIAL ALLPAX 


FTY= 35.000KSI, FTU# H2.000KSI 
FTV= 50.000KSI, FTU= A2.000KSI 
SEATING STRESS: 3.TOOKSI 

YIELD STRENGTH# IO.OOOKSI 
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Fig. 5-33 - Flange Comparison 3 
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Fig. 5-34 - Flange Comparison 4 
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Fig. 5-35 - Flange Comparison 5 
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MXOUS OXYGEN 900 PSI, t.S IN DIA 


PRESSURE 300 PSIG 



0.3I2A 

0.3T3D1A, 38 MOLES 

9.6810IA 901FACE 
0.067R FILLET 



FLANGE MATERIAL AL 
BOLT MATERIAL AL 
BASKET MATERIAL ALLPAX 


FTV= S.OOOKSI, FTU* A2.000KSI 
FTV= 50.000KSI, F1U= M.OOOKSt 
SEATING STRESS* 3.T00KSI 
VIELO STRENGTHS 10.000KSI 
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Fig. 5-39 - Flange Comparison 9 
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Fig. 5-40 - Flange Comparison 10 
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GASEOUS OXYGEN «0 PSI, 7 IN DIN 



FLANGE MATERIAL AL FTYr 3S.000KSI, FTU» A2.000KSI 

BOLT MATERIAL AL FTYr BO.OOOKS1, FTW* 6E.OOOKSI 

BASKET MATERIAL ALLPAX SEATING STRESSa 3.TOOKSI 


YIELD STRENGTHS IO.OOOKS! 



Fig. 5-41 - Flange Comparison 11 
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PRESSURE 60 P5I0 



0.912* 

0.27S0IA, 17 HOLES 

0.580DIA SDTFACE 
0.062* FILLET 



FLANGE MATERIAL AL 
OOLT MATERIAL AL 

gasket material all pax 


FTVc 35 .OOOKSI, FTU = 52.000KSI 
FTVr 50.000KSI t FTU= 62.000KSI 
SATING STRESS* 3.700KSI 
VIELD STRENGTH^ 10.000XSI 
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Fig. 5-45 - Flange Comparison 15 
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8ELT.UIL D A T A i i , 0 A T A 3 

El T 006-RU B67- I o o3/ 1 J- I H : 52 : 5 1 
cycle iogi 


OpQOOl 

000002 

QnnnaJ 

GOO 

Cuu 

000 

6X(jT 0AB5 

hardcopy 

2 


THESI MTA Wrtf u*EI> 
pot EX AM fie 1 IN sect. r. 

00000** 

Duo 

test flange 

1500 

Pil t 8 IN Dl A I 

opooos 

GOO 

1 .5 +3 

3 • u 

. R 3 7 s 

1.75 

QQ0O06 

Guu 

1 .S 

2*0 

1.5 

2 «u 

ooooo? 

000 

& 5 

8 

b 


000008 

Quo 

u 1 

u 

3 1 


000009 

GOO 

3H7SS AT RT 

3 1 7 SS 

AT RT A28* 

AT HT CRES321-A 

0000 1 0 

000 





OOOOl 1 

00 J 

•BOO . 




0000 1 2 

Guu 

SPACE SHUTTLE BAIN ENGINE FLANGE BOOG PSI. 16 1 ^ 01A 

0000 1 3 

UJU 

8 • + 3 

16*0 

I • 

2. 

0000 1 1 

uou 

i .& 

2.0 

1.5 

2,0 

0000 ] s 

uuO 

8 8 

8 

- 1 


0000 1 6 

UJU 

• 3 

* 6 

• Q 5 


0000 1 7 

Guu 

u 1 

i 

3 


GOGO 1 8 

OoO 

A2B6 AT RT 

a286 

AT RT A28* 

AT RT K-SEAL 

0Q0019 

uuu 

M 




000020 

Guu 

♦ 6uU • 



- 
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NOMINAL pressure p* 
nominal DIAMETER 0 l * 
tuse thickness t* 
tube thIckn tol« dT* 
HEIGHT TO flELD HT* 


PROOF FACTOR PF» 

BURST FACTOR fc>F» 

safety factor fs» 

GASKET FACTOR GF ■ 


PROPERTIES OF TUdE 

material 



MATERIAL TABLE NO, 

I* 

B 


ELASTIC MODULUS 

£■ 

* 28O0O0O0 + 0B 

PSI 

PO I SSON'S ratio 

NU ■ 

■ 300 


OEnS 1 T Y 

RHQ ■ 

• 2 6 80 L0/CUB 1 C- 1 NCH 

THERM EXP COEFF 

ALFA* 

.9SOOOCOO-OS 

INCH/ 1NCH/F 

TENSILE YIELU STtt 

FTY* 

•3SOOOOOO+05 

PSI 

ULTIMATE TENS STk 

FTU* 

. 9000u0u0+05 

PS 1 


l&OQ.OOO PSI 
B « uOO INCH 
.938 INCH 
.000 INCH 
I.7B0 INCH 


I .BOO 
2.000 
I .BOO 
2*000 


PROPERTIES OF FLaNGE MATERIAL 


MATERIAL TABLE Nu. 
ELASTIC MODULUS 
POISSON'S RATIO 
density 

THERM exp colff 
TENSILE YIELD STh 

ultimate tens st« 


I* 

B 


E* 

.28000000*08 

PSI 

NU* 

• 300 


RmO* 

.2880 LB/CU01C-INCH 

ALFA* 

. 9SOOQOU0-0S 

INCh/INCH/F 

FTY* 

• 3 5000000+05 

PSI 

FTU* 

.90000000*05 

PSI 


PROPERTIES OF BOlT material 


material table no. 
elastic mooulus 

POISSQN-S RATIO 
OENS I T Y 

THERM exp coeff 
tensile yield STk 

ULTIMATE TENS SfR 


1* 

8 


E * 

• 2 0UOGOi>o + u8 

PS I 

NO * 

• 300 


RHO* 

.2aap LB/CUBlC-lNCh 

ALFA* 

•9SOOOOUO-OS 

InCh/Inch/F 

FTY* 

■ 1 3 1 DDOuu *06 

PSI 

FTU* 

.20000000*06 

PSI 


PROPERTIES of gasket MATERIAL 

MATERIAL TABLE no. I* B 

ELASTIC MODULUS E* .28000000+08 PSI 
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YIELD STRENGTH 

seating stress 

THERM £XP cOEFF 

coeff of friction 
width coefficient 

AlDTn COEFFICIENT 
GASKET THICKNESS 
SEALING STRESS RaTE 


KG* , HOOOOOOO+05 PSI 

SG* .1890yOuO*OS PSi 

ALFA- « 95000000*05. INCH/INCH/F 

MU— >300 

GAMU* 1*000 

GAMS- 1*000 

MG* • 0250 INCH 

SPs S.SOOO 


OPT IONS 

0 io 3 ‘ ooo°° 

NUMBER OF PHASES to be considered IN the ANALYSIS 

TEMPERATURE DIFFERENTIAL * -SOO.0O DEG F 
COMPUTED THICKNESS T» .437S INCH 

test flange isoo PSI, a in dia 


aqb* 

. 1063 

SQ-iN 

weight* 

16,979/ 

LB 

B» 

1.3487 

IN 

H* 

2. 1044 

IN 

R I* 

4,U000 

IN DI« 

RG* 

4*5183 

IN OGI* 

D G 0 * 

RB* 

4*9506 

IN OB* 

RF I L* 

.3750 

IN 

rspot* 

.0620 

IN 

dhole* 

.4425 

IN 

ospot* 

.9012 

IN 

NB* 

26 


bg* 

, 322 1 

IN 

HT* 

1 , 7S0G 

IN 

MFU* 

, 1 1 1 7oH8S+0S IN 

Z£ T A 1 * 

, 7 697617 6 + 00 

ZET A2» 

,2421 1 833*00 


8*0000 IN 
8.7146 IN 
9*3567 IN 
9.9012 IN 


IN-LB/ IN 
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TEST FLANGE 1500 PSI , 8 IN 01A 

OUTPUT of The sTkess and deformation analysis results 


VaR I ABLE 


BOLT-UP 

staRt-up 

OPtRAT IO n 

bolt force (lb* 


.37061 2S6+G& 

.98850169*06 

. J6ol 1698+06 

E QU I V APPL MOMENT (IN-LB/IN* 


•6 0 99o856*o9 

-•2767i795+ 0 5 

.91552638*09 

axial deflection cin* 


• OoOOOoOu 

.9529S9B9-0& 

.95295989-05 

RADIAL OEFlECTIOm (in* 


•75853278 ”G 3 

-.21659999-01 

♦30019199-02 

FLANGE ROTATION (RADIAN* 


• 1330 1 1 1 5*02 

-.60896899-02 

.20197832-02 

bolt stress < ps i > 


. 13**u6267 + q6 

.17673339+06 

• 13028552 + 06 

GASKET STRESS (PSI * 


«9o5323S9*05 

.9271 J927+Q5 

.28673023+05 

STRESS IN FLANGE TOP (PSU 


.50399101*09 

-• 1 93759S3 + 06 

•19920502*05 

STRESS IN FLANGE BOTTOM (PSI* 


- . 13593958*05 

-.58698989*05 

-.82202766*09 

STRESS RESULTANTS (LB/IN* 


•OoOoOOOG 

.31690625*09 

.31690625*09 

(LB/IN* 

NY- 

.22025599*09 

— .628926ii6 + 08 

•87152195+09 

( I N-LB/ I N > 

MX* 

• 1 3975q52*09 

-.61693181*09 

.2091 1 3u 1 *09 

( IN-LB/IN) 

MY- 

. 9q925 | S 5*03 

-. 18507959*09 

•612339u9*03 

ilb/in* 

OX- 

.25839098*09 

-.1 1827670*05 

*39132061+09 

STRESSES AT NECK (PSI* INNER 

SI G X " 

-.2o697o79*05 

.9982 J22S+u5 

-.26289259*05 

(PSI* INNER 

S I G Y» 

-.26852167+09 

-.72199953*05 

•953B8235+09 

(PSI* OUTER 

SIGX- 

.2o697b79*o& 

-.89698225+05 

•56919258+05 

(PSI* OUTER 

SIgY- 

.97333907*09 

- • 1 2906639*06 

•23399879*05 

n (PSI* max 

TaU* 

♦62001? 15+09 

-.28386906+Q5 

.93916995*09 


ShOT-DOWn 


•3S98692Q+Q6 

• 1 o 7 1 0®5 3 + yS 

• 1 9299396-09 
• 1876uM9-Q2 
.23571190-02 
. 12838699*06 
.22793359*05 
• 129S1992+Q5 

-.20H7y799*05 

.97960938*09 

.59975279*09 

.23879926*09 

.71636279*03 

.95781067*09 

-.29085098*05 
-• 22875950*09 
• 99272S98*q5 
.19719689*05 

• 10987956*05 
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Section 6 
CONCLUSION 

The foundations have been laid for a simple but comprehensive design 
procedure for low profile flanges with a subsequent stress and deformation 
analysis. The algorithms have been programmed and the format for the 
basic output, i.e.» a summary of the flange geometry and a summary of the 
analysis results have been established. 

The computer program is set up for relatively few options of flange con- 
figurations within the class of low profile flanges. The amount of programming 
was limited by the number of man hours available for this contract. 

From the accompanying stress analysis it becomes quite obvious 
whether a design is sound or whether some basic design parameters need to 
be changes, such as the type of the gasket. The program is not automatic in 
the sense that it makes selective design decisions, which normally originate 
in the designer's mind based on his experience. Such a design procedure falls 
under the category of design optimization from the operations research stand- 
point (Ref. 37). The method described in Ref. 37, however, could be 
automated and combined with the current design/analysis program. The few 
material data and bolt geometry data currently incorporated in the design/ 
analysis program would then have to be expanded to large varieties. This can 
be done with the current program without any modifications to the existing 
logic. The current lists would just be longer having more entries. 

Further work is needed in verification testing. A test procedure to 
verify the moment carrying capacity of the flange, covering the entire range 
from elastic stresses to the formation of the plastic hinge in the flange neck, 
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is needed. These tests should he carried far beyond the initial yielding. 
Permanent strains in the highly stressed regions as well as permanent 
rotations should be measured versus applied moment. These tests can be 
carried out with an unpressurized test fixture since the entire loading can 
be expressed as an equivalent externally applied moment on the flange. 


The test should be carried with highly instrumented specimens of the 
following diameter sequence: 4, 8, 12, 22 and 45 inches, as shown on some 
of the examples in Section 5. Three pressure levels should be considered 
which have yet to be defined. High pressure levels for the small diameters 
and lower pressure levels for the large diameters are recommended. 

The specimens should first be tested without bolt holes, then with bolt 
holes, but without spotfaces, finally with spotfaces. Thus the weakening 
effects and the stress -raising effects of both the bolt holes and the spotfaces 
could be measured. Finally the fillet should be machined off and a groove 
as described in Section 2 be established. 


The instrumentation should include strain gages on the inside and the 
outside of the shell wall and the flange to verify essentially the stress dis- 
tributions shown on the plots labeled "Summary of Analysis." The strain 
gages should be mounted between bolts and in line with the bolts. Further 
instrumentation is needed to measure bolt force and gasket contact stress. 
Finally the deformation measurements, rotation and axial displacement, 
require some optical devices, possibly mirror systems. 

Further analytical work should proceed along the lines of a three- 
dimensioned elasticity solutions for a typical slice of the flange (see Fig. 6-1), 
requiring a three-dimensional finite element network. Lockheed -Huntsville's 
structural network analysis programs have not been made operational to 
include this type of analysis, although it would take only moderate further 
development effort to make the appropriate program modifications. 
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Hopefully the computer program delivered under this contract will 
be useful to the designers who are meant to use it. As more user's experience 
is accumulated it will definitely be necessary to make changes and improve- 
ments. The accompanying documentation in this report is provided for this 
purpose. 
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Section 8 
NOTATION 

coefficients for quadratic equation for h 
cross-sectional area of the flange 
amplitudes of the stresses 

stress area of the bolt 

total bolt area 

total gasket area 

flange cross-sectional area used for weight computation 

inner and outer tube radius, respectively 

radial lever arm between gasket and bolts 

bending rigidity of the shell wall 

creep constant 

width of the flange 

effective width of the flange 

gasket width 

effective gasket width 

width of the seal gland 

integration constants of the shell equation 

axial lever arms between centroid of flange and flange 
neck 

distance of spotface from shell outer surface 
constants in creep law 

equivalent rotational spring constant of the gasket and 
the bolts 

equivalent rotational spring constant of shell and flange 
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D 

d 


B 


hole 


e b> e g 


e 

e 


e i’ e 2 

F F 
tu’ ty 

f f f 
B’ T F 


3 

G 


I 

K 

K 


K , k 
P P 


B 


L. 

L 

i 


B 


determinant of the coefficient matrix of the shell-flange 
flexibility equation 

nominal diameter of the bolt 
diameter of the bolt hole 
elastic modulus 

elastic modulus of the bolts and of the gasket, respectively 

base of the natural logarithm 
radial lever arm 
radial spacing 

ultimate tensile strength and tensile yield strength, 
respectively 

bolt force, tube force, and force on flange face, respectively 
per unit length of radius 

height of the flange 

conductance parameter 
gasket thickness 

depth of the recess 

depth of the seal gland 

moment of inertia of the flange cross section 
strength 

equivalent constant in law governing interface leakage 
slope of sealing-f orce-vs-pressure curve 
shell parameter 

equivalent spring constant of the bolts 

equivalent spring constant of the gasket 

wave length of axial variation of stresses 
length of the leak channel 
strained length of the bolt 
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pipe bending moment 

exponent in interface leakage law 
applied flange moment 

edge moment 

meridional bending moment 

circumferential bending moment 

creep exponent 
number of bolts 

axial stress resultant 

circumferential stress resultant 

creep parameter in Lawson- Miller creep law 
bolt force 

total force on flange required 
gasket force 

force caused by internal pressure 

pressure 
volume leak rate 
edge shear 

meridional shear stress resultant 

stress ratio for fatigue design 
radius of a point in the shell wall 
radius of the shell wall middle surface 

equivalent radius of gasket and bolts spring constant 

bolt circle radius 

radius of the flange centroid 

fillet radius on the upper surface of the flange 
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r 


F 



r . 

1 


r 

s 


r 

spot 

r 

w 



s 

T 

t 


n 


t 

rupt 

u 

vol 

W 


w 

w 

X 



a 

a, (3 


radius of the application point of the force acting on the 
flange face 

gasket radius 

inner radius of the tube 

radius of the seal contact surface 

fillet radius for the spotface 

radius used for weight computation 

pipe shear force 

elastic section modulus of the flange 

line load on the gasket 

circumferential spacing 
tempe rature 
wall thickness 
thickness of the gasket 

part of tube wall thickness required to carry axial force 
due to pressure 

time to rupture 

axial displacement of the flange 

volume added to the tube by the flange 

weight leak rate 

width of the leak channel 

radial displacement 

axial coordinate 

tensile yield strength 

plastic section modulus of the tube and the flange, 
respectively 

linear thermal expansion coefficient 
creep constants 

stresses relating V and to when yielding occurs 
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y 

Vv^z 

Ap 

AT 

At 

AW 

Aw 

A (T-rjj A Of 


G 

e 

^ 1 * ^ 2 
r] 0 ,t) r n 2 

p 

p 

p 

<T 

a B 


X- xz 

a, a e 

<P 

X 

4» 

(F.s.) 


dimensionless parameters containing shell and flange 
stiffnesses 

angle between cylinder axis and weld 

coefficients relating the physical gasket width to the effective gas- 
ket width at yielding and under operating conditions, respectively 

internal pressure difference 

temperature difference between flange and tube 

tube wall thickness tolerance 

weight added to the tube by the flange 

difference in radial displacement between flange and shell 
change in stress of bolts and gasket, respectively 

deflection of the bolts 

deflection of the gasket 

strain 

coefficients related to the state of stress in the flange neck 
design parameters for bolt spacing 

coefficient of friction 
Poisson's ratio 
density of the material 
shell parameter 
stress, contact stress 
stress in the bolts 

contact stress of the gasket 

axial circumferential and radial shear stress, respectively 

equivalent stress 

circumferential angle 
roll angle 

weld weakening factor 
factor of safety 
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Appendix A 


A.l TUBE THICKNESS t (in.) OR (mm) 


(a) Given 


(b) 


t = 


where 


(F.S.) p r. 


ty 


(F.S.) 

P 

r. 

l 



a factor of safety 

the design pressure (usually proof pressure) 

radius of the inside of the tube, equal one-half 
of the inner diameter, d., of the tube 

tensile yield strength of the tube 


(B.F.)p r 

(c) t = T ■ ■ + 2 At 

f. a 4> 

tu 


t 


2 

where 


(P.F.) p r 
c 1 + 2 At 

f a 4* 

ty T 


(B.F.) 

(P.F.) 

* 

At 


p l 

p 2 


safety factor for burst condition 

safety factor for proof condition 

0.70 ... 1.00 (weld reduction) 

manufacturing tolerance for the tube wall 
(approximately 0.01 in) 

burst pressure 
proof pressure 

ultimate tensile strength of the tube 
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(d) 


t 


2 


F 


T 

ty 


t 


tensile yield strength of the tube 
maximum of t^ and t^ 


1.1 (B.F.) r. 
F* - 0.4 (B.F.) p L 


1.1 (P.F.) p 2 r. 
F^ u - 0.4 (P.F.)p 2 


where the symbols have the same meaning as in (c) 


.2 BOLT SIZE d B (in) OR (mm) 

Initial estimate dg -= t 

from bolt table for the proper wrench clearance and the given 
d„ find from Tables 2-1 through 2-3 or Tables 2-4 through 2-6 

the following quantities 

V ^1’ ^2’ A oB’ d hole’ r spot’ Size 

(a) Machined Spot Faces 

6 1 = ^l d B 
e 2 = T7 2 d B 

r fi j (fillet radius) according to Table 2-7. 

(b) Machined Groove 

e l = ^2 d B 
e 2 = e l 

r fil (groove radius) according to Table 2-8. 
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A. 3 BOLT CIRCLE RADIUS r g (in) OR (mm) 

(a) Machined Spot Faces 

- r. + t + c, + e, 

B 1 11 

where 

Cj = 0.0625 in or 1,5 mm 


(b) Machined Groove 


r 


B = 


r. + t + 2 r,.. + e. 
i fil 1 


A. 4 FLANGE WIDTH b (in) OR (mm) 

b = r B + e 2 - r i 

A. 5 GASKET WIDTH b Q AND GASKET RADIUS r Q (in) OR (mm) 


Estimate for gasket radius, r^: 

, , d hole , _ 4 

r G = 1/2 (r B T~ + T i ] 


Gasket Width, b^: 


(PF)pr G 

l a) b a * 2[yj k g . r 2 o a (FG)j 


or 


(PF) P r G 

(b) b G = 2 [y L K G - y 2 Kp(PF) P (GF)] 
whe r e 

y^ = a width factor for the gasket under initial deformation 
K g = the yield (crushing) strength of the gasket 
y^ = a width factor for the gasket under operating condition 
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= seating stress of the gasket 
(G.F.) = gasket factor 

kp = ratio of seating stress over pressure 


(c) b_ = r 0 - r. - 2 c_ (the available space is used) 

(jr Jd 2 1 2 


where c^ is a tolerance, c^ = 0.05 in or 1.0 mm. 


The force required to pre-form the gasket is 


P G ’ = 2 * r G b G n K G 


and the force required to keep a zero-leak connection is 


P G ’ “ P p + Z * r G b G r 2 °G <G - F - ) ° r 
P G 2) = P p + 2Tr G b a y 2 k p (P.F.)p(G.F.) 


where 


P p = * r G p (P ' F -) 


With the width b^ computed above the condition 


p(l) _ p (2) 
G “ G 


should be met. 


Gasket Radius, r • 
Lj 


(l) Gasket Close to Inside of Tube 


r G = r i + b G /2 + c 2 
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Check for space: r = + b^.,/2 + c^ 

r „ = r „ - d, . / 2 - c_ 

2 B hole' 2 

if Tj > r 2 set 

r G * r i + b G /Z+c 2 

r B = r G + b G /2 + d hole/ 2 + c 2 

(2) Gasket Close to the Bolts 

r = r - dhole - ^9 _ c 

G B 2 2 2 

Check for space as under (a). 

.6 PRESSURE ENERGIZED SEAL, EQUIVALENT b Q AND r Q (in) OR (mm) 


Estimate b„ = t 
Lr 

find 

r G = r l + b G /2 

check if space for seal gland is sufficient: 


r l * r i + b G + 2 b s 

r 2 = r B d hole/ 2 
if rj > r^ set 


‘B = r G + b G /Z + 2 b s + d hole/ 2 


where b g is the width of the seal gland. 

In both cases (A. 5 and A.6) the width of the flange has to be 
recalculated using the new bolt circle radius, 


b ' r B + e 2 ' r i 


as under A. 4. 
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A. 7 REQUIRED BOLT FORCE P g (lb) OR (N) 


(a) Flat Gasket (see A. 5) 


Pk 1 ’ - r G b G n k g 


k 2> * I p (P.F.) + 2, r G h Q r 2 o G (G.F.) 


or 


P< 2) = n r 2 p (P.F.) + 2jr r Q b Q y 2 (P.F.) k p (G. F. ) 


Pg = maximum of p B ^ and pj^ 


(b) Pressure Energized Seal 


P B = * r G p (P,F,) 


A. 8 NUMBER OF BOLTS, n 


B 


n 


B 


B1 “ F B A 

F ty oB 


n 


(B.F./P.F.) P B 

B2 = 

tu oB 


where 


n B = maximum of n BJ and 


"R 

F = tensile yield strength of the bolt 
ty 

F B = ultimate tensile strength of the bolt 
tu 


A. 9 BOLT SPACING s (in) OR (mm) 


s — 2 jr r B / n B 
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if s/dg > 8 decrease bolt size (if possible) and go back to A. 2. 

if s/d„ < n increase bolt size (if possible) and go back to A. 2. 

' B 'o 

A. 10 FLANGE HEIGHT h (in) OR (mm) 

6 = r B ' r G 

r = r. + t/2 
o 1 ' 

‘n = f / 2 

Ultimate moment to be carried 

(F.S.) P B e 
m Fu = 2 ITT 

Subtract effect of bolt holes from the flange width 

b = b ' d hole y d hole/" a " 

Assume 

tj = 0.8 

C 2 = 0.18 

and compute 

A = F f y Vfcr o 

B = F fy l 2 (*' - W/ 2 

C = F f y L <* 2 - l N >/ 4 - m Fu 
R 2 = B 2 - 4 AC, R = VS? 
h = (R - B)/2A 
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where F = tensile yield strength of the flange. 


A. 11 CHECK FLANGE HEIGHT 


if s/h >3 h = s/3 


to prevent waviness of the flange when too thin. 


A. 12 WEIGHT AW (lb) OR (kg) 

2 r. + t + b 
r w = —hi — 

A = (b -t) h 
w 

Vol = 2jr r A 
w w 

AW = p F • Vol 

where p^ = weight density of the flange. 

NOTE; Materialdata for some flange, bolt and gasket materials are 
given in Tables 2-10 through 2-13. 
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Appendix B 

SUMMARY OF THE ANALYSIS METHOD 




LOCKHEED- HUNTSVILLE RESEARCH & ENGINEERING CENTER 

i ■ ! , • . 

■; 1 . , ; . * ■! I . 1 ' 



LMSC-HREC TR D306492 


Appendix B 


B.l APPLIED FORCES 


= p r Q /2 (axial force in the tube wall) 
f r = phr./r^ (radial force on the flange) 


f 

P 




)/2 


r 

o 


(pressure force on the flange face) 


where 


P 

r 

0 

r. 

1 

h 


r 


G 


applied pressure 

radius of the tube wall middle surface 
inner radius of the tube 

flange height 
gasket radius 


B.2 SPRING CONSTANTS 


(a) Gasket 


A G = 2,r r G b G 
K G - A G V 2 * r o h G 


(Linear spring constant for 
flat gasket) 


K_ = A r E /Z-n r (h + h ) (Linear spring constant 
' J 1 r o R s for cantilever flange with 

pressure energized seal) 

where 

Aq = gasket area 
h^ = gasket thickness 
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h = depth of the recess 

h = depth of the seal gland 
s 

E^ - elastic modulus of the gasket 
G 

E.^ - elastic modulus of the flange 
F 

(b) Bolts 

A B ” n B A oB 

K-r, = A_ E_/2 ir r j (L inear spring constant for 
B B B ° B the bolts) 

where 

- total bolt stress area 

A -n = stress area of one bolt 
ors 

n fi = number of bolts 

2 = stressed length of the bolt. 

B 

(c) Equivalent Rotational Spring of Bolt6 and Gasket 
r B - r G 

r G ■+ K b °/< k b + K g ) 

K„ K e 2 / (K r + K r ) (Rotational spring constant 
B G for bolts and gasket) 

where 

e = lever arm between bolt circle and gasket circle 

r = radius of centroid of combined springs 
a 
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(d) Equivalent Rotational Spring of the Flange 



c 

r 

c 

B 

k 

P 

D 


bh 


A f h 2 /l2 
h/2 

r. + b/2 
E t 3 /l2(l -V 2 ) 
12(1- V 2 )/r 2 t 2 

Bk r o r c/ E F *F 



+ P (c Z 




2 



(Rotational spring constant for 
bolts and gasket) 


where 

A„ = cross-sectional area of the flange 
F 

b = flange width 

T = moment of inertia of the flange cross section 
F 

r = radius of the centroid of the flange 
o 

B = bending rigidity of the tube wall 
k = shell parameter 
(3 = flange parameter 

D = determinant of coefficient matrix of equation 

for interface bending moment and interface shear 
force at the flange neck 


(e) Constants for the Determining w q , m Q and q Q at the Flange Neck 

c = l/2k 3 B 
w 
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c 

m 


d ( 2 1 2 + p 4 + 2 ° k ) 


q ’ + 

B.3 INITIAL TORQUING (o) 


m F°’ = e V 2,r 

r Q (Applied flange moment) 

x <01 - 4°’/^ 

(Flange rotation) 

4°’ * v a b 

(Bolt stress) 

°G 0) * V A G 

(Gasket stress) 


Variables at the flange neck: 

n x°^ = ® (Axial force) 

= c m m^ (Meridional bending moment) 

(Shear force) 

= c (q^°^ - k m^°)) (Radial deflection) 
w x x 

n^°) = t w^°V r 0 (Circumferential force) 

m^°^ = V m^°^ (Circumferential bending moment) 
v x 


B -4 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306492 


PRESSURIZATION (p) 

AD 


m: 


f (r - r)+f (r - r ) (Applied flange moment 

P a p x a o . » 

v from pressure) 


Aw = 


m 


( 2 ) _ 


p r Q (1 - v/Z) 
E^Tt 


3 (c + l/2k) Aw 
D 


r r f 

+ r cirr* AT - ° C A - (Difference in 

° .Ej.p radial deflection 

between tube and 
flange) 

■ (Equivalent applied 
flange moment due 
to Aw) 


^ = (f^ + f v )/(K r + K b ) (Axial displacement) 


(p) 

F F 


C E + C F 

*< T) 

= x (o) + x (p) 

6 g 

= n<P> + X <P) % 

6 b 

= u< P) - X (p > K g 

a< T > 

G 

■ 4°’ + e g V 


(Additional flange rotation due 
to pressurization) 


_(T) _ (®) j. r ft /t 

ct B “ a B + E B 6 B' £ B 


k B' 


(Total stress in the gasket) 
(Total stress in the bolts) 


Variables at the flange neck: 
V (P) 

x x 


n_’_* ' = f (Axial force) 


(T) 

m F 

- c F X< T » 

(T) 

m 

- . m (T> 

— c m 

X 

m F 


r ™( T > 

- c q m F 


(Shear force) 
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(T) 

w 

= 0 (q<°> - 

W X 

( T ) 

E t . w< T > 

n 

y 

r 

o 

m (T) 

m 

y 

(T) 

- V m 
X 


,(<>)) 

X 


p r* {1 - V /Z) 

r -hr a T AT (Radial 

■“T 1 o l deflection) 


- a,j AT +p j (Circumferential force) 
v (Circumferential bending moment) 


B.5 STRESSES IN THE FLANGE 


w top 

w bottom 

a top 

^bottom 


w 

w - h X 


E F w to|/ r o 
E F w bottom^ r o 


B.6 STRESSES IN THE FLANGE NECK 


a 

x 



6 m n 

= + 1 +JL 

1 t 2 fc 


XZ 


= 1.5 q/t (max) 


where 

m = meridional bending moment 

m^ = circumferential bending moment 

n = axial force 
x 

n = circumferential force 
Y 

q = shear force, 
x 
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B .7 VARIATION OF THE SHELL VARIABLES ALONG THE TUBE 


n x (x) = n x (o) 


w (x) = c e 

X w 


n 


_(o) coskx - k m x (o) (coskx - sinkx)J 


P *1 0 -|) 

+ Lt »,°t m 


r <U°) . "1 

m x (o) {coskx + sinkx) ^ — sinkx 

q (x) = e x |^q x (o) (sinkx - coskx) + 2 k m x (o) sinkxj 
n (x) = Et w x (x)/r o - r o o< T AT + p r Q v/l 


-k 


m (x) = V m (x) 
y' x 


The stresses at a point x are then computed according to 
paragraph B.6. 


B.8 PLASTIC HINGE 


a, 


t 2 - t 2 F 
n E 


£ k + 

k A F 4k 3 


J— f -L i w , . 

: 7 i \ w + ^ 


t + t 


ou - 


n 


c + 


2k 




a = yi + a x + aj + 3 a 
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1/a 


a 2 /a 


b h 2 /Z 


FS x *1 2 

r-'-t <* 

o 


*n» + E 2 - * n ) 7 


(Ultimate applied 
flange moment) 
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INPUT INSTRUCTIONS FOR DESIGN AND ANALYSIS PROGRAM 
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Appendix C 


Format 

Column 

Description 

Units 

12A6 

1-72 

Instruction to plotter operator 

— 

15 

1-5 

Number of cases 

— 

12A6 

1-72 

Title of the plots 


E 10.4 

1-10 

p = pressure 

pai (N/mm^) 


11-20 

d. = inner diameter of the tube 

in. (mm) 


21-30 

t = thickness of the tube 

in. (mm) 


31-40 

At = thickness tolerance 

in. (mm) 


41-50 

h^, = height of tube frustum 




being part of the flange 

in. (mm) 

El 0.4 

1-10 

P.F. = proof factor 

— 

E 1 0.4 

11-20 

B.F. = burst factor 

— 

E10.4 

21-30 

F.S. = factor of safety 

— 


31-40 

G.F. = gasket factor 

— 

15 

1-5 

i^, = tube material number 




{see Table 2-6) 

— 


6-10 

ip = flange material number 




(see Table 2-6) 

— 


11-15 

ig = bolt material number 




(see Table 2-6) 

— 



iQ = gasket material number 




(see Table 2-7) 




Only if i T = 0 




Material properties of tube 


E 10.4 

1-10 

E^. = elastic modulus 

psi (N/mm^) 


11-20 

Vrj, - Poisson's ratio 

— 


21-30 

P T = density 

lb/in^ (kg/mm 
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Card 

6 


7 


8 


9 


10 


Format 

Column 

Description 

Units 

E10.4 

1-^ 

1 

o 

a^, = thermal expansion 
coefficient 

in/in/°F q 
( mm/mm/ C) 


41-50 

T 

F^ = tensile yield strength 

psi (N/mm 2 ) 


51-60 

T 

F^ = ultimate tensile strength 
Only if i F = 0 

psi (N/mm 2 ) 



Material properties of flange 
(same format and description 
as card 6) 




Only if i B = 0 




Material properties of bolt 
(same format and description 
as card 6) 




Only if i G = 0 

Material properties of gasket 


E10.4 

1-10 

E_ = elastic modulus 
Li 

psi (N/mm 2 ) 

E10.4 

11-20 

Kq = compressive strength 

psi (N/mm 2 ) 

E10.4 

21-30 

= seating stress 

psi (N/mm 2 ) 


31-40 

a G = thermal expansion 




coefficient 

in/in/°F 


41-50 

= friction coefficient 

(mm/ mm/°C ) 



Only if i G < 0 
Pressure energized seal 


E10.4 

1-10 

h = depth of the seal gland 

in. (mm) 

E10.4 

11-20 

b g = width of the seal gland 

in. (mm) 


21-30 

h^ = depth of the recess 

in. (mm) 

1015 

1-50 

Options : 

Option 1=0: 

Read tube thickness from 
card 3: 

— 



Option 1 = 1: 

Tube thickness computed 
according to Appendix A, 
paragraph A. 1(b): 
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Card Format Column Description 

10 1015 1-50 Option 1 = 2: 

Same, but paragraph A. 1(c) 

Option 1 = 3: 

Same, but paragraph A, 1(d) 

Option 2 = 1: 

Machined spot faces 

Option 3 < 0: 

Flat gasket 

Option 3 > 1: 

Pressure activated seal 

Option 4 = 1: 

Open wrenching (see Table 2-1 
or 2-4) 

Option 4 = 2: 

Socket wrenching (see Table 2-2 
or 2-5) 

Option 4 = 3: 

Internal wrenching (see Table 2-3 
or 2-6) 

Option 5 = 0: 

Gasket width according to 
paragraph A. 5(a) 

Option 5 = 1: 

Gasket width according to 
paragraph A. 5(b) 

Option 5=2: 

Gasket width according to 
paragraph A. 5(c) 

Option 6^3: 

Gasket close to bolt circle 


Option 6 = 3: 

Gasket close to inside of tube 


11 

2A6 

1-12 

Name of tube material 


2A6 

13-24 

Name of flange material 


2A6 

25-36 

Name of bolt material 


2A6 

37-48 

Name of gasket material (or 
seal, as applicable) 

12 

15 

1-5 

Number of loading phases 

13 

E10.4 

1-10 

Temperature differential between 
tube and flange 
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Units 


°F <°C) 
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